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SYNOPSIS 


The work presented in this thesis deals with heat transfer 
in buoyancy induced flows. The flow is a result of interaction 
between gravitational force and local density differences in 
the fluid. The density variations are due to temperature 
differences between a bounding surface and ambient fluid. 

Chapter ii of the thesis is devoted to analyse the flow 
and heat transfer over a surface which is itself curved or 
is parallel to the direction of gravity field/ while rest of 
the chapters are concerned with flows arising adjacent to 
surfaces that are either horizontal or inclined at a small ; 

angle to it. 

Chapter I is introductory and contains a brief outline 
of the background that have a direct bearing on our present 

t 

investigations, A review of the previous work of authors | 

mainly related to present study is also included in this , 

chapter, [ 

1 

I 

In chapter IIj, the laminar boundary layer description • 

f 

of planar/ free convection flow and heat transfer around two- | 

i 

dimensional bodies of arbitrary contour is examined for | 

I 

arbitrarily prescribed surface temperature distributions. The j 

[ 

Gbrtler-type transformation is developed to treat the general I 

wall thermal conditions. Series representations of the 



i± 


dependent boundary layer variables are developed on the lines 
of Gortler series method when the wall temperature variations 
are prescribed as 

S • 

k=0 ^ 

where x is a suitably normalized streamwise distance, a^ 
are constants with a / 0 and 

6 =0 for round- nosed cylinders 

and 

6 =1 for vertical plate. 

For round-nosed cylinders, the ordinary differential 
equations which determine the functional coefficients in the 
series expansion for the dependent variables are independent 
of the body contour and wall temperature data. Instead, some 
of the numerical coefficients in these equations are fixed 
by the class of body shape (namely, round-nosed here). Hence 
the solutions can be obtained for this class once and for all 
for fixed values of Prandtl number. The body contour and 
wall temperature data enter the problem only when the series 
are ass^nibled. Only two-terms of the series are found to be 
sufficient to represent velocity and tanaperature distributions 
around the lower-half of the circular cylinder when applied 
to specific wall-temperature distributions. The results for 



the heat- flux density at the surface of the non- isothermal 
circular cylinder compares well with previous analytical 
results. The effect of non- isothermal wall condition for 

V 

the circular cylinder influences significantly the heat 

transfer variations at the surface. When this method of 

solution is applied to isothermal circular cylinder^ the 

five- term sum of the series for heat- flux- density at the 

surface yields results accurate upto 87 percent of. the 

angular distance of 150® from the lower stagnation point of 

the cylinder and in the rest of the portion of it, the 

1 

results predicted are within an error of 3^ percent from the 
available exact numerical results. Also, heat transfer at 
the surface of an isothermal blunt elliptic cylinder of 
aspect ratio 2:1 is predicted quite accurately by five- term 
sum of the series for it. In the case of non- isothermal 
vertical plate problem, 24 functional coefficients associated 
with the first five-terms in the series expansion for 
dependent variables are calculated for Prandtl number 1 .0 
and 0.7. Using these functions which are independent of 
the prescribed wall temperature data, the boundary layer 
flow over non- isothermal vertical plate with sinusoidal wall 
temperature variation is studied. The results obtained for 
local heat transfer coefficient compares well with available 
finite difference data* 
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In chapter III, the boundary layer in free convection 
flow over the upper surface of a semi— infinite heated hori- 
zontal plate inclined at a small angle <X with the horizontal 
is studied analytically for flow and heat transfer by local 
non-similarity method of solution. The range of a considered 
here, as obtained from the analysis, is defined by the require- 
ment that A , the limit of tan a as Gr be 0(1) , 

where Gr is a suitably defined Grashof number. The conser- 
vation equations are transformed such that they can lend 
themselves to local non- similarity solutions. The forms of 
the local non— similarity method considered here are two and 
throe— equation models, obtained by retaining the transformed 
conservation equations as such and only selectively neglecting 
terms in the derived subsidisiy equations at the second and 
third level of the successive approximations. Numerical 
results for the local surface heat- transfer, wall shear-stress 
and velocity and temperature distributions are presented for 
Prandtl number 0,72 and 1,0 when A = 1,0, -1,0, While 
the two-equation model is found to give boundary layer 
characteristics sufficiently accurate for engineering purposes, 
the three-equation model over-corrects these solutions and 
the results thus obtained compare well with previous analysis. 
When the angle of inclination is negative, our method of 
solution predicts the separation of the flow and notably for 
the case of air, the two-equation model locates the separation 
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-I 

point within an error of percent from the previous exact 
value* It is to be noted that there is no evidence of 
singularity at the point of separation* 

In chapter IV, we derive in a systematic way all possible 
similarity solutions for both steady and unsteady free convection 
boundary layer flowsover a semi- infinite horizontal plate* The 
solutions presented are obtained not by assuming specific wall 
temperature distributions, but rather by deriving the solution 
forms from similarity conditions* This approach not only 
results in more general forms than hitherto obtained, but 
results in all other possible solution forms* Our investi- 
gations point out that there are in all four possible cases of 
similarity for free convection boundary layer flow over a 
heated semi-infinite horizontal plate* While this analysis 
has verified that the similarity possibilities for steady 
cases have essentially been covered in the literature, it also 
shows that somewhat more general wall thermal conditions could 
also be taken into account* The other two new possible cases 
deal with unsteady conditions* 

In the last chapter of the thesis, we initiate the 
study of free convection boundary layer flow over a semi- 
infinite, non- isothermal horizontal plate* A new Gortler 
type transformation is introduced to treat the problem of 
non- isothermal free convection with general wall temperature 



distribution. A new parameter characterizing the non- 
isothermal conditions of the wall is introduced* The 
transformed non-similar boundary layer equations are 
solved by the Gbrtler series expansion and the local non- 
similarity method* The series expansion is valid for an 
arbitrary surface temperature variation given by 

k=o ^ 

where n,p/G, are constants, G ^ O and x is suitably 

-K O 

normalized streamwise coordinate* 

Functional coefficients determining the terms in the 

series expansion for stream function and temperature depend 

only on n,p and Prandtl number of the fluid* The functions 

associated with first four-terms of the series expansion for 

dependent variables are determined for n = 0^ Pr = 0*72 and 

p =s 1,0, 0*5, 0.25* The application of the series method 

is illustrated for a class of wall temperature distributions, 

G (1 + e x^) where e is some constant* On the other 
o 

hand, the local non- similarity method at the second level of 
truncation is applied for the above class of surface 
temperature variations when Pr = 0*72, 5*0, 0*1 and 
p=l,0, 0,5, 0*25* A comparison of results for local heat 
transfer at the surface for air indicates that the local 
non-similarity method gives accurate results for a larger 
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region in the downstream flow than the four-term sum of 
the series expansion can give* Ihe effect of non- 
isothermal wall is significant on heat transfer at the 
surface and the variation- in heat transfer at the surface 
is greatly influenced by Prandtl number and the shape 
parameter, p, of wall temperature distribution* 



CHAPTER I 


INTRODUCTION 



X 


HEAT TRANSFER PROCESSES 


Heat transfer is defined as the transfer of energy 
across a system boundary caused solely by a temperature 
difference. The study of heat transfer is but the study 
of how energy is carried, distributed, and diffused in and 
by materials. Such a study is of significance in the 
energy-related applications,. For instance, the transfer 
of heat in power plants from the energy source, be it nuclear, 
solar or other, to the working fluid is one of the most 
basic processes in such systems. Similarly, the operation 
of refrigeration and aii>conditioning units depend on the 
effective transfer of heat in condensers and evaporators. 
Other applications of particular interest pertaining to 
environmental control include the minimization of building- 
heat losses by means of Improved insulating techniques and 
use of supplemental energy sources, such as solar radiation, 
heat pumps, and fire places. 

There are two basically different processes by which 
heat transfer is accomplished: conduction and radiation* 

Heat conduction is identified as the process of molecular 
transport of heat in bodies (or between thorn), due to 
temperature variation in the medium. Thermal radiation, 
on the other hand, depends only on the temperature and on 
the optical properties of the emitter, with its internal 
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energy being converted into radiation energy* 

The concept of convection embraces the process of heat 
transfer through the motion of a fluid* Heat is transported 
simultaneously during the process by convection and conduction 
By convection is meant the process of heat transport occurring 
through the movement of macroparticles of the liquid or 
gas in space from a region of one temperature to that of 
another* Convection is possible only in a fluid medium and 
the transport of heat is directly linked with the transport 
of the medium itself. Thus a consideration of problems in 
convective heat transfer involves the science of fluid motion* 
The practical importance of this mode of heat transfer needs 
no illustration* The use of moving fluids to transport or 
remove heat is well-known^ ranging from the circulation of 
coolant through a nuclear reactor to the mounting of a 
power transistor on a block with cooling fins. 

The convective heat transfer process is further divided 
into two basic processes* If the motion of the fluid arises 
due to a foixiing condition such as an externally- imposed 
flow of a fluid stream over a heated object, the process 
is termed as forced convection* If# the other hand, 
no such externally- induced flow is provided and the flow is 
driven simply by the interaction of a difference in density 
with a body force field, such as the gravitational field# 
then the process is termed natural convection or free 
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convection. The agencies which provide the density differences 
can also be very diverse* For instance^ they may be due to 
temperature an^/or composition variations* But in our present 
investigations, we deal only with temperature caused density 
differences* Also, our work deals only with free convection 
flows arising in a very extensive medium due to a given 
inhomogenity in toriperature between a bounding surface and 
ambient fluid* 

As a result of convective processes, local tCTiperature 
differences exist in the flow field resulting in higher 
temperature gradients* In most applications, it is of 
primary interest to know the quantity of heat exchanged 
between the surface and the stream* For a surface of area A, 
a heat transfer coefficient, h, is defined by the equation 

Q = hAe^, , (1*1) 

in which Q denotes the heat flux through the area A and 
is a characteristic temperature difference between the 
surface and the fluid* Assximing that the properties on whcih ' 
the heat transfer depends can be considered as constants 
with the exception of the density, the variation of which is 
expressed by a thermal expansion coefficient (Boussinesq 
assiimption), Nusselt [l] derived the following equation for : 
the heat flow# Q# per unit time 

, d^P^g 30 : 

Q * kd e / 5 — - ) • (1.2) 

w ^<^P m2 . 
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Here, g is the gravitational acceleration, 3 is the thermal 
expansion coefficient, d is a characteristic length, c 

XT 

is the specific heat of the fluid, k is the thermal conductivity 
P is the density and jn is the viscosity* Again, he extendec 
the thermal analysis to a gas in which the properties depend 
in the following way on the absolute temperature T 



k = 


o 


P 


= S® ) 

o 


(1.3) 


The index, 0, refers to a characteristic fluid temperature. 
The heat transfer coefficient h, then obtained by him is 



pr 



(1.4) 


when a » b = c. 

The names given to the dimensionless parameters appearing 
in the equations (1.2) and (1.4) are 


Nusselt niimber, Nu = ^ # 

fi c 

Prandtl ntomber, Pr = , 


and 


2 3 

gj3p"^d e 


w 


M' 


Grashof number, Gr 


(1.5) 
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The expression used to describe heat transfer in free 
convection can now be written as 

Nu = f(Gr, Pr) (1.6) 

corresponding to (1.2). When a^b,c are all distinct, 
the relation describing heat transfer is 

T 

Nu = f(Gr, Pr, ,a,b,c) . (1.7) 

o 

The form of (1.7) is an evidence to the fact that dimensionless 
relations describing heat transfer become more restricted, 
the more accurately one tries to describe the fluid properties. 
On the other hand, the relations of the form (1.6) hold for 
groups of fluids with the same Prandtl number. They are 
exact in the limit of small temperature differences. Apart , 
they describe with fair approximation, conditions at larger 
temperature differences when average fluid properties are 
introduced. Sparrow and Gregg [ 2 ] developed a reference 
temperature relation that permits good approximations of 
variable property results to be obtained quite simply from 
constant property results. 

2 . GOVEJRNING EQUATIONS FOR FREE CONVECTION 

The theoretical analysis of free convection requires 
that the fundamental laws of mass, momentum and energy and 

t 

the particular laws of viscous shear and conduction be 
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utilized in the development of mathematical theory for the 
free convection phenomena* A detailed derivation of the 
governing equations for the free convection phenomena# 
based on a Newtonian fluid with Fourier's law governing 
the molecular flux of heat# is available in [3] and therefore 
a familiar description will be selected as a starting point 
in our thesis* The difficulty of analysis in considering the 
full equations results from the coupling of momentum and 
energy equations and where important, from the variations 
in the properties of the fluid P # /i and Tc • An analysis 
including the full effects of variations in p, jj. and k is 
so complicated that some approximations become essential* 
Consequently the equations governing free convection flow are 
commonly used in a form obtained with the help of the 
Boussinesq approximations [3] • 

In the Boussinesq approximation, variations of density P 
is ignored except in so far as they give rise to gravitational 
force and the dependence of P on temperature T is given by 

p -p=P3(T-T^) (1.8) 

where suffix 'a' stands for the values of the quantities in 
the ambient medium* Under this assumption it has been 
formally justified by Ostrach [4] and Heliums and Churchill 
[S] that variations of all fluid properties can be ignored 
ccOTpletely* Ostrach has also shown that, within the framework 
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of this approximation^ viscous dissipation of energy may be 
neglected. For very low speeds found in free convection the 
effects of compression in the energy equation can also be 
neglected [ 6 ] . 


Under these assumptions, the governing equations for 
the free convection phenomena without internal heat genera- 
tion are 

V . V = 0 , (1.9) 


P ^ = - ^p + ju. V +pg , (1.10) 

^ + V. VT = a v^T (1.11) 


where 


a 



a is known as the thermal 


diffusivity, p the local pressure, T the local temperature, 
V = (u,v,w) is the local velocity vector and t is the time. 


In the force-momentum equation (I.IO) the local pressure 
p may be broken down to two terms, one due to the motion 
of the fluid, P, known as the viscous pressure and the other 
due to the hydrostatic pressure p^^. That is, we have 


P = P + Pa . 


It can be seen that the hydrostatic pressure with the body 
force acting on the fluid constitutes the driving mechanism 
for the free convection flow. 

pg - vp 


The vector 
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in the momentum equation (1,10) can be eiqjressed as 
,0 g - vp = p g - ( P^g + V P) , 

a, 

where we have used 

vPa = pj 

in the gravitational field. 

Using (1,8), finally we get 

p g - Vp = (P-P^) g - vP , 

= - ppdVT ) g-vP, (1,12) 

cl 

Hence the momentum equation (1,10) can now be written as 

~=-ivP + 3/V^V-g 3(T-t^) . (1.13) 

Here -g 3(T~T )# to be denoted by B, is known as the 
buoyancy force. Equations (1,9), (1,13) and (1,11) 
constitute the basic equations of free convection in the 
Boussinesq approximations, 

3, BOUNDARY LAYER APPROXIMATION 

An understanding of the detailed mechanism and 
interaction of flow and heat transfer processes was made 
possible through the application of boundary layer concept, 
introduced for forced convection flows by Prandtl in 1904. 
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According to this concept, there exists a thin region near 
solid objects in which certain quantities describing the 
flow change very rapidly. The idea is to consider the 
relative size of certain expressions in the differential 
equations governing the flow and discard some, due to their 
relative smallness* This leads to an impressive simplifi- 
cation of the governing equations for the flow* It is to 
be noted that in forced — flow problems, the validity of 
this approximation depends on the size of the Reynolds 
number for the flow and the theory becomes a better approxi- 
mation as the Reynolds number increases. 

However, this important development was brought into 
free convection phenomena only in 1930, The experiments 
by E, Schmidt and Beckmann [?] on an isothermal vertical 
surface in air indicated that boundary layer concepts of 
Prandtl were also applicable to buoyancy induced flows. 

However, there exists no unique prescribed characteristic 
velocity for natural convection flows, Reynolds niamber, 
therefore, loses its significance. So a new characteristic 
velocity must be defined which, as was first explicitly 
pointed out by Ostrach involves the three prime physical 

quantities of this phenomenon, namely, 0,g,0^* Unfortunately, 
the magnitude of the velocities in natural convection can 
not be determined from the boundary conditions, so that a 
characterization of the reference velocity must come from 

* 
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the differential equations* The explicit form of this 
velocity depends on the specific problem. 

For instance, in order to obtain non-dimensional 
equations for natural convection over a vertical plate , 

U = (g^e d)^/^ (1.14) 

^ w 

is taken as reference velocity scale and it was shown that 

Grashof number Gr, plays the role of Reynolds number in 

. 2 

this free convection flow and in fact Gr is simply Re 
based on U as given by (1.14). Also we note that the 
pressure in the region beyond the boundary layer in natural 
convection flows is hydrostatic unlike in the case of forced 
flow where it is being imposed by an external flow. Several 
other experiments have also upheld the validity of this 

a 

boundary layer approximation in free convection flows. 
Boundary layer theory has now been applied with success to a 
very large number of transport phenomena in buoyancy induced 
flows • 

4. THE BACKGROUND OF PRESENT INVESTIGATIONS 

In many interesting and important areas of natural 
convection, the flow is generated as a consequence of the 
thermal input from a surface x>rhich is itself curved with 
respect to the direction of gravity field. Surface curvature 
requires that the boundary layer must continually adjust 
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to the changing strength of the surface-parallel component 
of the buoyancy force. Thus an analysis of free convection 
flow including the effect of surface curvature is more 
complicated than that for a vertical plate. 

One of the earliest theoretical investigations to 
appear on free convection flow over a curved surface is that 
of Hermann [s] for cylindrical configuration. By means 
of a suitable transformation/ he was able to solve the 
boundary layer equations for an isothermal circular cylinder. 
The agreement between the calculated temperature field and 
the experimental results of Jodlbauer [9] was good. However/ 
his approximate method of solution is limited to isothermal 
case and cannot be extended to more general cases. Also/ 
it is not clear how other cylindrical objects are to be 
treated. Schuh [lo] found the natural convection solutions 
in similarity form for certain curved surface geometries. 

He considered in detail only the flow over two-dimensional 
body having everywhere finite curvature and a boundary 
layer of constant thickness. Later, Merk and Prins [ll] 
made a detailed study of both symmetric two-dimensional and 
axi-symmetric flows. They showed that the cone is one 
axially symmetric body having a similarity solution. It 
was pointed out that neither the horizontal cylinder nor 
the sphere has a similar solution. They obtained a 
similarty solution valid near the stagnation point and later 
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presented an integral solution for horizontal cylinders and 
spheres, Braun/ Ostrach and Heighway [12] generalised the 
previous two analyses and investigated such flows for 
similarty for both constant and variable surface tanperatures» 
A group of closed bottom shapes was found for each kind of 
symmetry. Herring and Grosh [13] and Herring [ 14] found 
similar solutions for the vertical cone for power-law 
surface temperature dependence. 

Although horizontal cylinders and spheres do not have 
similarity solutions for natural-convection flow, their 
importance has resulted in a nimber of approximate analyses 
also based on boundary layer equations, Chiang and Kaye 
[is] used a Blasius series expansion method to predict the 
heat transfer coefficients for laminar natural convection 
frem horizontal cylinders and spheres. For forced flow 
situations. Van Dyke [l6] evaluated the Blasius series for 
a parabolic cylinder of nose radius R and found that the 
resulting series diverges for 3</R > 0,62 where x denotes 
the usual coordinate measured from the nose along the arc. 

Also Schlichting [6] cautions that even the series for a 
circular cylinder may become divergent from a certain 

value of x/R onwards since the coefficients of the series 

9 

appear to increase after the term of 0 (x/R) • 

Recognising the need for a procedure to account for the 
effects of the body shape in a general manner, Seville and 
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Churchill [l7] used a G'^rtler type series for analysing the 
two-dimensional, planar and axL- symmetric boundary layer flows 
over smooth isothermal bodies of fairly arbitrary contours* 
When applied to horizontal circular cylinders or spheres, the 
Gortler series converges faster than the corresponding 
Blasius series [l5] . The ccxnputed temperature and velocity 
field at 90® and 150° from the stagnation point showed 
reasonable agreement with experimental data even if only 
one term of the series were used. Due to lack of information 
of the universal functions associated with the higher order 
terms of the series expansion for flow variables, no definite 
conclusion can be drawn on the ability of this method to 
predict heat transfer accurately for other cylindrical 
configurations. Again, the same problem was treated by 
Lin and Chao [is] using a series expansion of the form 
proposed by Merk [l9] and corrected by Chao and Fagbenle [20]. 
Their results agreed well with the previous analyses [l5, 17] 
for horizontal circular cylinders and spheres. The main 
difficulty in this method of solution is to evaluate the 
gradients with respect to the configuration parameter of 
the problem, which necessitates accurate evaluations of 
Zero-order solutions in the neighbourhood of each value of 
that parameter and hence calculations are to be performed 
on a large scale. Finite difference solutions have been 
obtained by Merkin [21, 22] for horizontal circular and 
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elliptic cylinders with uniform surface temperature or 
uniform heat flux. Lin and Chao [23] demonstrated that 
their series solutions for elliptical configuration agreed 
satisfactorily with the finite difference data and also 
pointed out that the various series results are largely 
dependent on body shape and much less on Prandtl niimber and 
coordinate transformations for the variables. Recently/ 
Muntasser and Mulligan [24] applied local non-similarity 
method of solution [25 ] to study analytically the free 
convection over an isothermal horizontal circular cylinder; 

This method of solution gave results comparable in accuracy 
with the finite difference data. If the boundary layer 
characteristics are required only at specific points/ the 
local non— similarity method [25] represents a substantial 
saving in computer time over full-scale digital methods. 

Recently interest is shown in solving the governing 
equations for free convection over isothermal/ horizontal 
circular cylinder without making boundary layer approximations. 
Kuehn and Goldstein [26] solved the governing Navier-Stokes 
and energy equations and obtained solutions over a wide range 
of Rayleigh number/ that is# 10® £ Ra- < lo”^/ where neither 
asymptotic matching techniques nor boundary layer assximptions 
are accurate. Their results also include the development 
of the buoyant plume (the plume begins to form at distance ot 
from the stagnation point/ given by a > 150®) which can not 
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be obtained using boundary layer methods . Pustovalov and 
Matveev [ 27 ] and Farouk and GUceri [28] have also solved 
these equations numerically* 

For free convection along non- isothermal surfaces, 
analytical solutions have been found for the situations 
which yield similarity solutions. However, for laminar 
free convection exterior to heated surfaces, there are 
a very few analyses known for arbitrary non-uniform surface 
temperatures* The free convection over a vertical plate 
with arbitrarily prescribed non-uniform surface toriperature 
was analysed by Sparrow [29] using the approximate integral 
method and later by Scherberg [30] , utilizing the classical 
polynomial profiles of Squire [3l] • Sparrow has treated a 
class of special surface- temperature variations [29] for the 
entire Prandtl number range, while Scherberg considers the 
non- isothermal surface problem by means of a quasi— uniform- 
surface temperature scheme, in which the effect of non- zero 
leading edge boundary-layer thickness is also taken into 
account* Kuiken [32] obtained an expansion for the case 
where wall temperature variation is expressed as a power 
series in terms of the streamwise coordinate. A Gdrtler- 
type series expansion has been tried by Kellher and Yang 
[ 33 ] to accommodate more general wall temperature variations. 
They demonstrated that the integral procedure of Sparrow [29] 
over-estimated values for heat transfer. Recently, 


gave 
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Kao/Domoto and Elrod [34] have developed another approximate 
method of solution for analysing free convection over non- 
isothermal vertical plate with arbitrary ( though smooth) 
temperature or heat- flux distributions, while Na [35] obtained 
finite difference solutions for several wall temperature 
distributions and he reported that the series solutions of 
Kao#Domoto and Elrod [34] are in good agreement with his 
data only in the initial portion of the boundary layer and 
further downstream they deviate from his computed results. 

The free convection over non— isothermal curved surfaces 
have not received much attention. Koh and Price [36 ] solved 
the non— isothermal free convection problem over circular 
cylinder as formulated by Chiang and Kaye [is] . They 
eliminated the dependence of polynomial coefficients in the 
surface condition from the governing equation of the flow 
by using special transformations. But still the number of 
terms needed to represent velocity and temperature profiles 
is directly dependent on the number of terms retained in the 
wall temperature data. Recently, Lin and Chem [37] 
extended Merk's type series method of solution to non— 
isothermal surfaces and presented numerical solutions to 
horizontal c ir cular cylinder for a specific wall— temperature 
distribution. Now, this method of solution involves 
algebraic manipulations also, besides the need to evaluate 
accurately the gradients with respect to the redefined 
configuration parameter. 
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In chapter II of the present thesis, our aim is to study 
the laminar boundary layer description of planar free conve- 
ction around two-dimensional isothermal and non- isothermal 
bodies of arbitrary shapes • The flow is a result of 
temperature difference between the surface of the object and 
the fluid. The Gortler-type transformation already proposed 
for isothermal surface condition [l7 ] is modified to deal 
with arbitrarily prescribed non- isothermal wall conditions. 
The resulting transformed boundary layer equations contain 
two streamwise dependent parameters which characterize both 
the non- isothermal conditions and the shape of the object. 
Series representation of the dependent boundary layer 
variables are developed on the lines of Gdrtler series method 
[ 38 ] , when the wall temperature variations are prescribed as 

( i) S a-.x'^ for round— nosed horizontal cylinders, with 

k=0 ^ 

a^ O and x, the distance from stagnation point 

OO - 

and ( ii) x S a^x^ for vertical plate, with a /= 0 and 
k=0 ^ 

X, the distance from the leading edge. 

round— nosed cylinders, the ordinary differential 
equations which determine the functional coefficients in the 
series expansion for the dependent variables are independent 
of the body contour and wall temperature data. Instead, 
some of the numerical coefficients in these equations are 
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fixed by the class of body shape, namely roxind-nosed here* 
Hence solutions can be obtained for this class once and for 
all the fixed values of the Prandtl n-umber. In the 
isothermal case, the first five -terms of the series are 
split into 24 universal functions and the governing 
equations for them are integrated for Prandtl number 0,7 
and 1.0 on a Dec-10 computer with reasonable accuracy. 

Again, for the non— isothermal case, the first three-terms 
of the series are split into 16 universal functions and as 
before the governing equations for those universal fxmctions 
which account for non-isothermality conditions are integrated 
for Prandtl numbers 0,7 and 1.0. The body contour and the 
wall temperature data enter the problem only when the series 
for flow variables are assembled. 

For an isothermal circular cylinder, it is found that 

a 5 -term s\am of the series representation for shape 

parameter accurately represents it over the lower-half of 

the cylinder while the results obtained for the heat-flux 

density at the surface with the same number of terms from 

the corresponding series for it are accurate upto 87 percent 

of angular distance of 150° from the lower stagnation point 

of the cylinder and in the rest of the portion the results 

1 

are predicted within an error of 3~ percent in comparison 
with the available exact results, thereby showing that this 
series method is valid as far as 150°. On the other hand. 
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for blunt 'elliptic cylinder with aspect ratio 2 : 1 , the shape 
parameter is well represented by a 5 -tenn Siam of its series upto an 
eccentric angle 45 ° measured from lower stagnation point 

and the results for heat-flux density at the surface obtained 
from sum of the first 5-tgrra^ of., its series is sufficiently 

accurate upto 115° on comparison with finite difference 
data* For non—isothermal circular cylinder, the method of 
solution is applied to two specific wall ’temperature 
distributions* For practical purposes, two terms are found 
to be just sufficient to represent accurately the velocity 
and temperature distributions in the vicinity of the lower half 
of the circular cylinder* Also, it is observed that non- 
isothermal conditions at the surface has a significant effect 
on the rate of heat transfer from the surface to the fluid. 

For vertical plate problem, we calculated 24 universal 

functions associated with the' firstf-terms in the series^-n^' 

f 

expansion for flow variables. These universal functions 
depend only on Prandtl number and are independent of wall 
temperature data. These functions are calculated for 
Prandtl number 0,7 and 1.0 and the botmdary layer flow over 
non -iso thermal vertical plate with sinusoidal wall tCTiperature 
variation is studied,. The results obtained for local heat 
transfer coefficient cempares well with available finite 
difference data [35]y 
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In the study of natural convection phenomena^ another 
area of importance is the flow over flat surfaces which are 
inclined to the direction of the gravity field. Sixrface 
inclination results in a pressure gradient across the 
boundary layer and leads to an analysis more complicated 
than that for vertical surface- For small angles of 
inclination with the vertical/ the pressxire gradient term 
in the momentum equation corresponding to the direction 
along the_ surface is neglected by an order of magnitude 
arguments. But retaining of this pressure gradient tem 
becomes very important for boundary layer over surfaces ■ 
inclined at large angles to the vertical. For surfaces 
which are horizontal or nearly horizont§l/ the buoyancy force 
tangential to the surface is very small. In particular, 
there is no component of the buoyancy force along a 
horizontal surface and there is therefore no direct drive 
of any tangential flow that may result. Rotem and 
Claassen [39] carried out a flow visualization experiment 
using a semifocussing Schlieien system which clearly 
indicated the existence of a boundary layer near the 
leading edge on the upper side of a heated horizontal s^lrface- 
Recently, Yousef ^ Tarasuk and Mckeen [40] also studied the 
laminar free convection boundary layer over a heated 
isothermal surface using Mach-Zehnder interferometer. 
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Stewartson [4l] was the first to give a theoretical 
description of how a botmdary layer is formed over a 
horizontal surface lander the action of buoyancy force. 

His analysis contains a sign erjzor in the pressure gradient 
term of the horizontal momentum equation and thus condition 
for existence of such a boundary layer is that the heated 
platen is to face downwards * Later, this error in the flow 
analysis was corrected by Gill, Zeh and del—casal [^42 ] 
and one finds a horizontal boundary layer above the heated 
plate.- driven by an induced pressure, which falls in the 
direction of flow; This is called '"indirect drive"" and 
is due to the pressure deficiency in the thermal transport 
region^next to the s'urface with respect to the ambient 
medium* It is _ indirect in that the buoyancy force does not 
drive the flow. Rather, it produces a pressure field which 
induces such a flow. 

Free convection from nearly horizontal surfaces placed 
in an extensive body of flxiid overlaying it above is vital 
to meterological and industrial applications. Pera and 
Gebhart [43 ] conducted ej<periments to determine laminar 
natural convection flow above horizontal and slightly inclined 
stirfaces. For horizontal surfaces,iheir experimental findings 
were in close agreement with those of Rotem and Claassen [39]. 
They observed that foxr nearly horizontal surfaces with an 
angle of inclination a to the horizontal, the boxindary layer- 
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region is shorter for a <4° while, at higher inclinations 
the flow persists further. Also, they studied analytically 
the effect of a small surface inclination on flow and heat 
transport as a perturbation of flow over a horizontal 
surface. They reported that the effect of surface incli- 
nation is reflected more strongly in the velocity rather than 
in the temperature field. For a heated plate inclined at 
a small angle a to the horizontal, the boundary layer flow 
over the surface is driven not only by the indirect mechanism 
but also by a component of the buoyancy force along the 

surfacei Jones [44] studied this flow by considering that 
except in the neighbourhood of the leading edge, the direct- 
drive due to dominates induced pressure effects, which 

can therefore be neglected; The resulting flow far down- 
stream is then, in his analysis, what is descitibed by the 
classical free convection solution for a heated vertical 
plate. He developed two series to obtain the solution, one 
valid near the leading edge and the other at large distances 
from it. The region where neither series is adequate, a 
step-by-step numerical procedure was used to obtain solutions 
of the boundary -layer equations. His numerical method of 
solution to obtain boundary -layer characteristics of the flow, 
is the selected points method of lanczos and hence, necessarily 
he had to linearize the transformed conservation equations at 
any downstream location before commencing his numerical 
procedure. In the initial portion of this patching region. 
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he solved the set of equations obtained by using variables 
which account for departures from similarity solution fbrr 
flow over horizontal surface^ while at some downstream 
location onwards^ he switched onto another set of equations 
using transformations suggested by anology with the classical 
case of free convection from a heated vertical plate. 

In chapter the boundary layer in free convection 

flow over the upper surface of a heated semi -inf inite 

horizontal plate inclined at small angle c£ with the horizontal 

is analysed by the local non -similarity method of solution [25j, 

The range of a considered here, as obtained from the analysis, 

is defined by the requirement that A , the limit of 
1/5 , . 

'"'r tan a as_Gr-*“,be OCl), where Gr is a suitably defined 
Grashof number; The forms of the local non-similarity method 
considered here are the two and three-equation models obtained 
by retaining the transformed conservation equations as such 
and only selectively neglecting terms in the derived 
subsidiary equations at the second and third level of the 
successive approximations in this method of solution; The 
corresponding two-point non-linear boundary value^’ problem for 
these two models at a streamX'^rise location are solved by the 
initial value technique- as modified by Nachtsheim and 
Swigert £45]. Particularly, the application of this 
numerical procedure for the highly non-linear three-equation 
model was coinmenced at some downstream location with good 
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initial estimates of missing wall conditions generated by 
quasilinearization technique [46] • Nxamerical results for 
the local surface heat- transfer, wall shear^stress and 
velocity and temperature distributions are presented for 
Prandtl number 0,72 and 1 .0 when A =1/ -1# 

While the two-equation model is found to give boundary 
layer characteristics sufficiently accurate for engineering 
purposes, the three- equation model ovea?- corrects these 
solutions and the results thus obtained, compare well with that 
of previous analysis. For small positive angle of 
inclination of the surface with the horizontal, unlike the 
previous analysis [44] where a composite non— similar solution was 
developed, our method of analysis shows solutions for the 
entire region of flow can be obtained with a single set of 
transformed boundary layer equations which makes no assumption 
on the neglect of induced pressure from some distance downstream 
away from the leading edge of the plate* 

When the inclination of the plate to the horizontal is 

negative, although the pressure gradient associated with the 

indirect process remains favourable, separation of the 

boundary layer flow eventually occurs, since now opposes 

the motion* Our method of solution predicts the point of 

separation of the flow and notably for the case of air, the 

two— equation model locates this point within an error of 

li per cent from the previous exact value [44 ] • Also, there 
2 
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is no evidence of singularity at the point of separation of 
the flow, as also pointed out in [ 44 ] • 

In the study of free convection boundary layer flow 
over nearly horizontal isothermal surface, we have seen 
that the nature of the flow is non-similar, the departure 
from similarity arising due to the inclination of the 
surface at an angle a with the horizontal. However, when 
a = 0, the flow is similar and various analyses are available 
for this situation. Stewaptson [41] published solutions for 
the flow for Pr = 0,72, while Gill, Zeh, and del-Casal [42 ] 
extended this to Pr = 1,0 and 10.0. All these solutions 
refer to the case in which the temperature of the bounding 
wall is constant, while Gill [42 ] also gave the solution for 
the case in which the boundary layer temperature varies as a 
power of the distance from the leading edge. Rotem [47] gave 
similar solution for sami- infinite horizontal plates for a 
power- law temperature difference between the plate and a 
uniform ambient fluid. Rotem and Claassen [39] ntimerically 
integrated the governing equations of two dimensional steady 
flow over an isothermal surface and gave the velocity and 
temperature profiles for various Prandtl numbers. Also they 
formulated the s±milarity solution for axi- symmetric flows on 
infinite surface [48 ] but no numerical results were given* 
Blanc and Gebhart [49] obtained some exact solutions to 
natural convection flows adjacent to horizontal disks. 
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including the energy effects of pressure and, viscous 
dissipation. 

Pera and Gebhart [43] considered both the surface 
conditions of uniform temperature and uniform heat flux 
for free convection boundary layer flow over sani- 
infinite horizontal plate and computed nxomerical solutions 
for a wide range of Prandtl numbers, Clarke and Riley [50 ] 
considered the prototype problem of Stewartson but allowed 
for density changes and temperature-dependent viscosity 
and thermal conductivity. Also they allowed the heated 
surface to be permeable and considered the effects of 
blowing upon the flow and thus their analysis can be taken 
to model the evaporative nature of the fuel surface. 

Again# they considered the free convective flow induced 
when a semi- infinite horizontal fuel surface bums in a 
quiescent oxidizing atmosphere [51], Recently Mahajan and 
Geihart [52 ] analysed the higher order effects for natural 
convection flow adjacent to a horizontal surface for a 
power“law surface temperature variation with leading term 
in this perturbation problem as the well known Zeroth— order 
boundary layer similarity solution treated earlier by 
Rotem and others [39# 43] . 

Very recently a few analytical studies [53# 54# 55 ] 
dealing with unsteady free convection flow over semi- 
infinite horizontal plate are made and they all deal with 
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the cases when the plate temperature oscillates in time 
about a constant non-zero mean tonperature# While 
specific cases of similarity solution are dealt for steady 
free convection boundary layer flow over horizontal 
surfaces with different and additional effects# no 
similarity solution was considered so far for unsteady cases* 

It is well known that most approximate methods depend to 
some extent on knowledge of the '■'similar*^ solutions to 
the boundary layer equations* As we shall see in the final 
chapter# it is possible to formulate the boundary layer problem 
in such a way that the similar solutions form the first# and 
by far the most significant# term in an exact series solution 
to the partial differential equations, so that they can be 
used with considerable success even for non-similar flows. 

The first study of a similar boundary layer appears 
to have been that for flow parallel to a flat surface by 
Blasius [ 55 ] not long after Prandtl [57] first propounded 
the concept of a boundary layer* The pressure— gradient in 
the stream direction in the solution by Blasius is zero* 

The range of similar solutions was then extended by Falkner 
and Skan [58 ] who included non-zero values of the pressure - 
gradient. A major contribution to our understanding of the 
concept of similarity was made by Mangier [59] who generalized 
the similarity conditions to include a much wider range of 
pressure gradients than earlier authors* Yang made a 
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corresponding study for free convection over a vertical 
plate and cylinder and gave all possible forms of similarity 
for both steady and unsteady cases* 

In chapter IV we derive in a systematic way all possible 
similarity solutions for both steady and unsteady free 
convection boundary layer flow over a semi- infinite horizontal 
plate* The method consists of defining a new independent 
variable and three dependent similarity variables in terms 
of the generalised transformation variables. When these 
similarity variables are substituted into the boundary layer 
equations, they are reduced to a set of ordinary differential 
equations only if the coefficients of the similarity variables 
and their derivatives can be made constant* These coeffi- 
cients when made constant are referred to as similarity 
conditions/ they usually contain derivatives of the trans- 
formation variables and are thus differential equations 
themselves* Any solution to the similarity conditions which 
will result in explicit forms for the transformation variables, 
will transform the original set of partial differential 
equations into a set of ordinary differential equations and 
these are called similarity equations. 

We have concluded from our present study that there 
are in all four possible cases of similarity for free 
convection over a horizontal surface* This analysis verifies 
that the similarity possibilities for steady cases have 
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essentially been covered in the literature, though somewhat 
more general wall thermal conditions could also be taken 
into account* The other two new possible cases deal with 
unsteady flow problems and these situations have not 
received any attention so far* 

In the previous chapter for free convection over a semi- 
infinite horizontal plate, we have been concerned exclusively 
with similar boundary layer flows* But the majority of 
practical engineering problems involve flows for which condi- 
tions for similarity are not satisfied* Each problem 
involving laminar boundary layers has its own particular 
distribution of surface temperature and to solve the partial 
differential equations as opposed to the ordinary differential 
equations derived using similarity arguments, even for one such 
problem is a major computational task. Only recently, finite 
difference solutions [35] for non- isothermal vertical plate 
free convection problem have been obtained and that too, 
for some specific wall thermal conditions only* The other 
methods which yield solutions to arbitrary wall temperature 
variations belong to the series expansion techniques* 

Unlike the vertical plate case where the surface 
characteristics to be determined relate to shear stress and 
heat flux, the horizontal plate problan includes also of 
determining the induced pressure at the wall* For vertical 
plate problem^ methods ranging from Karman-Pohlhausen 
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technique with relatively crude profiles to full scale 
numerical computation procedures have been developed* On 
the other hand^ to our knowledge, no single study either 
analytical or experimental is attempted so far for the 
corresponding case of a horizontal plate* 

In the last chapter of the thesis, we initiate the study 
of free convection boundary layer flow over a semi- inf inite 
non- isothermal horizontal plate for an arbitrarily prescribed 
surface conditions* A new Gortler-type transformation is 
introduced to reduce the governing boundary layer equations 
into non-similar form containing a streamwise dependent 
parameter which characterizes the non- isothermal conditions 
at the wall* The resulting equations are solved by the 
GOrtler series expansion and the local non-similarity methods* 

Series expansion is valid for an arbitrary surface 
temperature variations of the form 

S G. 
k=0 ^ 

where n,p,G^ are constants and 0 and x is a 

suitably normalized streamwise coordinate* Universal 
functions associated with the first 4-terms of the series 
expansion for flow variables are determined for air when 
p = 1*0, 0*5, 0*25* Application of the series solution 
method is illustrated for a class of wall temperature 
distributions given by 
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G (1 + s yP) for some e. 
o 

Also local non- similarity method at the second level of 
truncation, namely the two-equation model, is applied for 
this special case to predict local heat transfer coefficient 
at the wall for Pr = 0.1, 0.72, 5.0 when p = 1.0, 0.5, 0.25. 
A comparison of results for local heat transfer at the surface 
for air indicates that the local non— similarity method gives 
accurate results for a larger range of values of s(= 1+ s yp) 
covered here than the 4- term sum of the series ej^pansion. 

The effect of Prandtl number on local heat transfer coeffi- 
cient at the surface for various values of p is studied 
by the local non— similarity method. It is found that for 
fixed values of s and p, the percentage increase of non- 
iso thermal heat transfer at the surface over that for an 
isothermal surface decreases as Prandtl number increases. 



32 


REFERENCES 


[1] Eckert, E.R.G. , ASME, j. Heat Transfer, 103, 409-414 

(1981). 


[ 2 ] 

[3] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 
[ 9 ] 

[Ife] 

[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 
[17] 


Sparrow, E.M., and Gregg, J.L. , Recent advances in 

Heat and Mass Transfer (J.P. Hartnett, editor), 
Mc-Graw-Hill Book Co., N.Y. (1961). 

Jaluria, Y., Natural Convection Heat and Mass Transfer, 
Pergamon Press, Oxford (1980). 

Ostrach, S., NACA TN 2635 (1953). 

Hellxims, J.D. and Churchill, S.w., Ch^. Eng. Progr. 
Symposium Ser., No. 32, 57, 75 (1961). 

Schlichting, H,, Boundary layer Theory, McGraw-Hill 
Book Company, N.Y. (1960). 

Schmidt, E., and Bechman, W. , Tech. Mech. Thermodyn., 

1, 341-391 (1930). 

Hermann, R., NACA Tech. Mem. 1366 (1954), 

Jodlbauer, K., Forsch. Ing.-Wes. 4, 157-172 (1933). 

Schuh, H., Gt. Brit. Ministry of Air Production 

Volkenrode Repts. and Transls. 1007, Apr. (1948). 

Merk, H.J. and Prins, J.A,, Appl, Sci. Res. A4, 11-24, 
195-206, 207-222 (1953-54). 

Braun, W.H., Ostrach, S, and Heighway, J.E., Int, J. 

Heat Mass Transfer, 2, 121-135 (1961). 

Herring, R.G., and Grosh, R.J., ibid 5,. 1059 (1962). 

Herring, R.C., ibid, 8, 1333-1337 (1965). 

Chiang, T., and Kaye, J., Proc. Fourth National Cong. 
Appl. Mech., University of California, 

Berkeley, 1213-1219 (1962). 

Van Dyke, M., Proc. Ninth International Cong. Appl. 

Mech., Brussels, 3, 318 (1957). 

Seville, D.A. and Churchill, S.W., J. Fluid Mech., 29, 
391-399 (1967). 





33 


[^18] Lin^ P*N# and Chao, B. T. , ASME^ J, Heat Transfer, 96, 

435-442 (1974)* 

[19] Merk, H*J,, J. Fluid Mech., 5, 460-480 (1959). 

[2 0] Chao, B.T., and Fagbenle, R.O., Int. J. Heat Mass 

Transfer, 17, 223-240 (1974). 

[ 21 ] Merkin, J.H., ASME-AICHE Heat Transfer conference, 

St. Loiuis, Mo., August (1976). 

[ 22 ] , ASME^ j. Heat Transfer, 99, 453-457 (1977). 

[ 23 J Lin, F.N,, and Chao, B.T., ibid, 100, 160-163 (1978). 

[ 24 ] Muntasser, M.A., and Mulligan, J.C,, ibid, 165-167, 

(1978). 

[ 25 ] Sparrow, E.M., Quack, H, and Boemer, C.J., AIAA Journal, 

8, 1936-1942 (1970). 

[26] Kuehn, T.H. and Goldstein, R.J*, Int. J. Heat Mass 

Transfer, 23, 971-979 (1980). 

[ 27 ] Pustovalov, V.N, and Matveev, yu. ya., Inzhenemo- 

Fizicheskii Zhumal, 43, 2 H982). 

[2 8] Faxouk, B. and Gtlceri, S.I., ASME, j. Heat Transfer, 

103, 522-527 (1981). 

[2 9] Sparrow, E.m., NACA TN 3508 (1955). 

[ 30 ] Scherberg, M.G., Int. J. Heat Mass Transfer, 8, 1319-1331 

(1965). 

[ 31 ] Goldstein, S. Ed. Modem Developments in Fluid Dynamics 

Clarendon Press, Oxford (1938). 

[ 32 ] Kuiken, H.K., Appl. Sci. Res., 20, 205-215 (1969). 

[ 33 ] Kelleher, M., and Yang, K.T,, Quart, J. of Mech. and 

Appl. Math., 25, 445-457 (1972). 

[ 34 ] Kao, T.T., Domoto, G.A. and Elrod, H.G., ASME J, 

Heat Transfer, 99, 72-79 (1977). 

[ 35 ] Na, T.y,, Appl. Sci. Res., 33, 519-543 (1978). 

[3 6] Koh, J.C.Y. and Price, J.F., ASME, J. Heat Transfer, 

87, 237-242 (1965). 



34 


[ 37 ] 

[38] 

I3'9] 

[40j 

[ 41 ] 

[42] 

[43] 

[44] 

[45] 

[46] 

[47] 

[48] 

[49] 

[50] 

[51] 

[52'] 

[ 53 ’] 


Lin, and Chern, S.Yi, ibid| 103# 819-821 (1981). 

Gortler, H., J. Math, Mech* 6# 1-66 (1957). 

Roten, Z,, and Claassen, L,, J* Fluid Mechi# 39, 

173-192 (1969). 

Yousef, W.W, , Tarasuk, J.D,, and Mckeen# W.J.# ASMS, 

J, Heat Transfer, 104, 493-500 (1982). 

Stewartson, K, , Z.A,M.P., 9a, 276-281 (1958). 

Gill, W.N,, Zeh, D.w,, and del-Casel, E,, ibid# 16, 
539-541 (1965). 

Pera, L,, and Gebhart, B., Int. J, Heat Mass Transfer, 
16, 1147-1163 (1973). 

Jones, D.R., Quart, J. of Mech. and Appl, Math., 26, 
77-98 (1973). 

Nachtsheira, P.R. and Swigert, P., Development of 
Mechanics, 1, 361-367 (1965). 

Roberts, S,M, and Shipnann, J.S,, Two-point Boundary 
Value Problems: Shooting methods American 
Elsevier NY (1972). 

Rotom, Z., Proc. First Canadian Nat'l Congr. Appi, 
Mech., 2, 309-310 (1967). 

Rotem, Z., and Claassen, L;, Canadian J. Chem, Eng., 

47, 460-467 (1969) i 

Blanc, P. and Gebhart, B,, 5th International Heat ' 
Transfer Conference, Tolcyo, 3, 20-24 (1974). 

Clarke, J.F, and Riley, N., Quart. J. of Mech. and 
Appl, Math,, 28, 373-396 (1975). 

, J, Fluid Mech,, 74, 415-431 

Mahajan, R.L. and Gebhart, B., ASME, J, Heat Transfer, 
102, 368-371 (1980). 

Verma , R.L. and Singh, P., Aust. J. Phys,, 30, 

335-345 (1977). 



35 


[54’] Singh, P., Sharma, V.P, and Mishra, U»N., Acta 
Mechanica, 30, 111-128 (1978), 

[ 55 ] Verma, A.R,, Int. J* Engng» Sci*, Vol. 21, 35—43, ( 1983 ) • 

[56] Blasius, H,, Z, Math, und Phys, 56, 1-37 (1908). 

[ 57 ] Prandtl, L., International en Mathematiker-Kongresses, 

Heidelberg, 481-491 (1904). 

[S'S] Falkner, V.M. and Skan, S.w., British Aero. Res* 

Council, Reports and Memoranda, 1884 (1937). 

[ 5 - 9 ] Mangier, W., Z.A.M.M., 28, 97-103 (1948). 

fso] Yang, K.T., Trans. J. Applied Mechanics, 82, 230-236 
^ (1960). 



CHAPTER II 


LAMINAR FREE CONVECTION OVER TWO-DIMENSIONAL 
ISOTHERMAL x\ND NON-ISOTHEIMAL BODIES OF 
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1 INTRODUCTION 

This chapter is primarily concerned with prediction of 
heat transfer in planar free convection boundary layer flow 
over smooth bodies of arbitrary shape with uniform and 
non— uniform surface thermal conditions# placed in an unbounded 
Boussinesq fluid. There is# at present# adequate information 
about laminar free convection boundary layers over isothermal 
bodies of simple shapes such as a vertical plate or cone# 
for which the flow is similar [l ] • However# similarity does 
not arise in many practical situations in the study of free 
convection phenomena. Notably# both cylindrical and 
spherical configurations do not give similarity. Hence# 
several other methods have been employed for obtaining a 
solution to the governing equations of such flows. 

Hermann [ 2 ] has treated analytically the laminar steady 
state heat transfer frem an isothermal cylinder by a modified 
Pohlhausen's similarity solution for the vertical plate at 
Prandtl number 0.733. The boundary layer thickness at 
different angles around the cylinder was obtained by multi- 
plying the flat plate boundary layer thickness by a parameter 
that is a function only of angle from the stagnation point. 
Merk and Prins [3 ] presented analytical results of the heat 
transfer around horizontal cylinders and spheres using the 
Integral method first introduced by Squire for a flat plate. 
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Chiang and Kaye [4] used a different scheme for obtaining a 
solution to the boundary layer equations for free convection 
over a horizontal cylinder with prescribed surface temperature 
or surface heat-flux^ expressed as polynomials of distances 
from the stagnation point* The approach was similar to the 
Blasius series for incompressible viscous flow over arbitrary 
cylinders* Again^ for free convection flow around non- 
iso thermal circular cylinder, Kbh and Price [ 5 ] devised 
special transformations to treat the differential equations 
that arise out of applying Blasius series technique so as to 
remove their dependence on the polynomial coefficients in the 
prescribed surface thermal conditions* Later, Chiang, Ossin 
and Tien [6] also employed the Blasius series method of 
solution to determine the velocity and temperature profiles 
for angular positions around a sphere. 

Recognizing the need for a procedure to account for the 
effects of body shape in a general way/Saville and Churchill 
[7] introduced an analysis closely patterned on Gdrtler^s 
series for forced flow [8] so as to describe laminar free 
convection near isothermal horizontal cylinders with fairly 
arbitrary body contours# When applied to horizontal circular 
cylinders and spheres, the series converges rapidly than the 
corresponding Blasius series* Again, Wilks [9] and Lin [lo] 
used the Gbrtler's series method of solution in studying 
two-dimensional boundary layers over bodies of uniform 


1 
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heat— flux* Also this series method has been tried by 

Kelleher and Yang [ll ] for treating free convection over 
a non- isothermal vertical plate with more general wall 
temperature variations • Kuiken [ 1 2 ] has developed a 
particular series solution for the free convection over a 
non- isothermal vertical plate* Lin and Chao [l3]used a 
Merk type series to obtain solutions for various two- 
dimensional geometries with horizontal circular cylinder 
as a special case* Later Lin and Chem [14] extended this 
method of solution for free convection over non— isothermal 
horizontal circular cylinders* Also, finite difference 
solutions have been obtained by Merkin [l5#l6] for 
horizontal circular and elliptic cylinders of uniform 
temperature or uniform heat-flux* Na [l7 ] presented finite 
difference results for heat transfer in laminar free 
convection over non- isothermal > Plate , 

In the present study, we have developed a Gortler typo 
transformation to study the laminar free convection boundary 
layers over non— isothermal bodies of arbitrary shape* The 
resxilting .transformed boundary layer equations contain two 
parameters dependent on streamwise locations on the surtace*- 
The special data of the problems under our consideration 
appear only in these parameters* Series representations 
for the dependent boundary layer variables are developed 
on the lines of Gbrtler series method of solution for forced 
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flows# This method of solution is applied to flow over 
non— isothemal vertical plate and long horizontal round- 
nosed non— isothermal cylinders of arbitrary body shapes. 

Wall temperature variations in these two cases are 
prescribed as follows: 

OO 

S a,x for round-nosed cylinders, with a_ ^ 0, a, 
}c=0 ^ ox 

constants and x, a normalized distance from stagnation point 

and 

(ii) X S a, X for vertical plate, with a ^ O and x, 
kisO ^ ^ 

a normalized distance from the leading edge# 

For round— nosed cylinders, the ordinary differential 
equations which determine the functional coefficients in the 
series expansion for the dependent variables are independent 
of the body contour and wall temperature data. Instead, 
some of the mmerlcal coefficients in these equations are 
fixed by the class of body shape, namely round— nosed here. 
Hence, solutions can be obtained for this class once and for 
all for fixed values of the Prandtl niunber. In the isothermal 
case, the firs“t five -terms of the series are split into 24 
universal functions and the governing equations for them are 
integrated for Prandtl number 0.7 and 1.0 on a Dec-10 
computer with reasonable accuracy# Again, for the non— 
isothermal case, the first . three -terms of the series are 
split into 16 universal functions and as before the governing 
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equations for those universal functions which account for 
non- isothermal ity conditions are integrated for Prandtl 
number 0*7 and The body contour and the wall 

temperature data enter the problem only when the series 
for flow variables are assembled. 

For an isothermal circular cylinder, it is found that 

a five- term sum of the series representation for shape 

parameter accurately represents it over the lower— half of 

the cylinder while the results obtained for the heat— flux 

density at the surface with the same member of terms from 

the corresponding series for it are accurate upto 87 per cent 

of the angular distance 150® from the lower stagnation point 

of the cylinder and in the rest of the portion the results 

1 

are predicted within an error of 3^. percent in comparison 
with the available exact results [15 3 /.thereby showing 
that this series method is valid upto 150®. On the other 
hand/ for blunt elliptic cylinder with aspect ratio 2:1# the 
shape parameter is well represented upto an eccentric angle 
of 45® measured from the stagnation point and the results 
for heat-flux density at the surface is sufficiently accurate 
upto an eccentric angle of 115® on comparison with finite 
difference data [16 ] «. For non- isothermal circular cylinder/ 
the method is illustrated for two sp€K:ific wall temperature 
distributions* It is found that non- iso thermal conditions 
at the surface have a significant effect on heat transfer*. 
Also the effect of Prandtl nximberis studied*. 
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For vertical plate problem# proceeding as before we 
calculated 24 universal functions associated with the first- 
five terms in the series expansion for flow variables for 
Prandtl number 0.7 and l.O, They are dependent only on 
Prandtl number and are independent of wall temperature data. 
Using these functions#, we investigated the boundary layer 
flow over non— isothermal vertical plate with sinusoidal wall 
temperature variation* The result obtained for local heat 
transfer coefficient compares well with available finite 
difference data [l?] • 

2. BASIC EQUATIONS 

We consider steady laminar flow over two-dimensional 
bodies of non-uniform surface temperature t^ situated in 
an ambient fluid of temperature t^^# where t^ > t^ • For 
the steady state free convection phenomena# equations (1*9), 
(1*13) and (1*11) reduce tos 

Continuity equation 

V * V = O # (2*1) 

conservation of momentum 

V. vv = -^vp+ V - ^(t-t^)g # (2 #2) 

consejTvation of energy 

V. vt = a t # (2.3) 
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where t represents the temperature. Using the identity 
V*W = i t'(v . V) - vx vxv 

in (2.2) and then applying curl operator on both sidos> 
the resulting momentum equation is 

V .vw = - 3 vt X g + W (2*4) 

where W = vxv. 


The momentum equation (2*4) is now free from pressure 
variable though/ of course/ we have increased the order 
of the system. 

The boundary conditions of the problem are: 


and 


V = 0/ t =s t^^ on the surface 

V s= 0, t = t^ at infinity • 


(2.5) 


3. BOUN^RY LAYER MODEL OF FREE CONVECTION IN 
CURVILINEAR COORDINATES 

In order to formulate the governing equations (2.1)/ 

(2.3) and (2.4) for a curvilinear system of coordinates/ 
we use the special system of coordinates which is intimately 
connected with the geometry of the surface namely a horizontal 
cylinder for our considerations here and also symmetry properties 
of the problon. For the curves x = constant we take the 
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normals to the wall| and for the curves y « constant we talce 
the curves parallel to the wall* The coordinate system is 
then an orthogonal one based on distances x and y where 
X is measured along the surface from an intersection of the 
plane of symmetry parallel to the body/ which for flow near 
horizontal cylinder/ is taken as the forward stagnation point 
while y is normal distance from the wall (see figure 1). 
Tollmien [is] derived expressions for the complete Navier- 
Stokes equations in this system of coordinates* We shall 
now derive the boundary layer model of the free convection 
phenomena on the lines of Goldstein *^3 work [19] for forced 
flow over a cylinder* 


Let <P be the angle measured in the anticlockwise 

direction between a horizontal plane and the line tangent 

to the surface at P(x/y)* The components of PQ along 

X and y directions are respectively h^dx and h 2 dy 

where the scale factors h^ and h^# are given by h^= I + k y 

and h^=: 1/ k being the local curvature of the wall so 

that K = — , a function of x • The general e 5 <pressions 
dx 

for the operators div, curl/ grad, etc*/ in this system of 
coordinates are 


V = 


JL 

1+ Ky 


-f 

1 



V .9 = £ ^ +■ % (1+k5) V } , 

1+ Ky dx ay 
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V X V = k ( 2 l| « ^ (l+ic^u J . 

l+tcy ax ay 

v2 = I i. (-i__ -L) + -Lj (2,^6) 

1+Ky ax 1+K^ax ay ay 

where u and v are respectively the velocities in the x 

and y directions, i, j, k are the unit vectors along 

X, y and perpendicular to the plane of the flow respectively# 

We now introduce the following non-dimensional variables 


X 




u 



V 


K* = icL, e 


t— t 




(2#7) 


where 


t.,- = (t -t ) G(|^), t being some reference tcmpe- 

rature# Here L is a characteristic length and G(x) is 
a prescribed surface temperature variation* 

Using (2*6) and (2.7), the governing equations (2.1)/ 
(2i4) and (2.3) become 

+ ly 

r ,, ^ + V fix « 9- (l + ic'y)u1 

^ ( i+K^y) ax ^ ayJTi+i?^ ^-ax ay ^ j 


I 
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g0(T^-.T^)L"*G(x) 

v2 


f cos <f> de , ^ „ ae i 
'•r+ic-y + sin <0 — J 




ay j (l + K^y) 


X (i+x'y)^] (2.9) 


u 30 

XT+PyT 3 ^ 


+ V 


ay 


U0 1 dG 

+ TThTTT ^ ^ 


a r 9- 
57n+Py) L alJ 




30^ 

3x^ 


. ^ a..‘y, If] 


(2*10) 


Here <P is the angle between the surface normal and the 

body force. The momentum equation (2.9) contains a dimension— 

gP(T -T^)l3 

less parameter ^r~~ # called the Grashof number and 

2/^ mm 

to be denoted* hereafter* as Gr. We shall write cr =-^ and cr 
is a dimensionless parameter of the fluid and this occurs in 
energy equation (2.10). 

We now introduce asymptotically modified variables in 
equations (2.8) to (2.10) 

y = y Gr ' r u = u Gr ' * V « v Gr ' 


Ihe introduction of these variables is equivalent to the 
application of boundary layer approximations. The conser- 
vation equations (2.8) to (2.10) now become 


^ + i. 


V = O * 


(2^1) 
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[ S |L + ^ 1_^[ 1 g-i-Ci4.«*fGrV*)a] 

(1 + ic'yGr ^ By (l + K^yGr Gr^'"" 3y 


■g(x)[-t-,S 2S-L__,^ |a + sin<l>2|] _ 


H -h ' 


u 


ay 


3 0 . * a£ + u e 1 


Cl+x'jGr""''*) 3y (1+t'yGr"^'^^) ® 


^ (l + K^^r 


1 r 1 . 1 ^ 

-^r-^/4) S?72 3x (i^.K'5Gr-^/^) 3=4 


(2»12) 


+ 4 {l+«'^r“^/^) ^ ] 

3y ay 


(2^13) 


As Gr becomes large and if the curvature of the body , it , 
is small in order to insure that s'Gr y remains small# 
the above system of equations will yield the inner solution 
valid near the boundary surface and known as the boundary 
layer solution to the free convection over the surface [l9] • 
Thus# by neglecting the terms which vanish supposedly as 
Gr -* °°s we arrive at the system 

al + ai =: o # (2.14) 

3x ay 

^ " 3* 

9 + ^ i-)(- 3H) = -G(x)sincP 21 - # (2.15) 

9^ ay ay 3y 3y 
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u 


3x 


+ V 


* 

3y 


^ ^ _ 1 
G dx “ 0 





with boundary conditions t 


Y 

Y 


= 0 : 

OO • 


u=0=v ,9=1# 

U-*0 ,0-*-O» 


(2^16) 


(2.17) 


The momentum equation (2.15) is integrated once with 
respect to y and then, using boundary conditions, we arrive 
at the final form of the boundary layer model for the free 
convection phenomenas 


3x ^ ~ ^ * 

^ BY 

u ^ +• V ^ = G(x) sin <P + , 

3y 


*30 . * 30 ^ 1 dG 1 8^0 

oY oY 


with the boundary conditions (2.17). 


( 2 .18 ) 


(2.19) 

( 2 . 20 ) 


With the introduction of the stream function S such 


that 
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3x ' 


( 2 . 21 ) 


the boundary layer model for two-dimensional surfaces can 


be written as 
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( 2 ♦ 22 ) 

(2 *. 23 ) 


y = 0: = 0,0=1, 

Y 

Y '* 0 ° i O, e-*0, (2 #24) 

y 

where the sine of the angle ^ between the body force 
vector and normal to the surface of the object is S(x)« 

It can be seoi that for <P = 7i/2, the above equations with 
the corresponding conditions are the boundairy layer equations 
to the laminar free convection over a non- isothermal 
vertical plate* 

We now define the important physical characteristics 
for the boundary layer flow and heat transfer* Ihe Nusselt 
number, Nu, for the dimensionless heat— flux density at the 
surface is defined as 

wn - (heat-flux density) (characteristic length) ^ 

(thermal conductivity) (temperature difference) 

= -Gr^/'^ (^)* (2*25) 

dY y=o 


f xy 


♦ - 


G(x) S(x) e + , 

yy yyy 


-w 0 +1^® Q-ifi 

* X *X ^ G dx ^ “ a 

y y y yy 


with the boundary conditions: 
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(2.26) 

(2.27) 


(2.28) 


We now introduce the following transformation to reduce 
the laminar boundary layer equations to a system of differential 
equations containing two streamwise dependent parameters that 
account for the body shape and non-uniformity in surface 
themal conditions: 



§ = r [g( 

O 

z) S(z)]^'^^ dz , 

(2.29) 


7) = (|)^/^ 

A 

[s(x)G(x)1 . 

(2.30) 

f (x/y) 

(4)3/4 g 

F(g/h) / 

(2.31) 

T(g,n) 

= 0(x/y) • 


(2'. 32) 


The Nusselt number averaged over a length/ L/ is 


Nu = i- r Nu(x')dx' . 
o 

The dimensionless shear stress at the surface is 


T (ai,_ , 


f V fl 


ay y=0 


and the average shear stress at the wall is 


f / T (x')dx' . 
f ^ o ^ 


4. TRANSFORMATION OP BOUNDARY LAYER EQUATIONS 
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It is to be noted that for isothermal bodies, the above 
transformations reduce to that of Saville and Churchill [?] 
in their treatment of boundary layer free convection over 
isothermal bodies of arbitrary contour* 

The transformed boundary layer equations are now 


^ = 1 - Vr>n'> - 1 

K(g ) P^ , 

(2.33) 

§ =1 

F77T , 

(2.34) 

where 



K(e) - i + U_ dCSG) 

2 ^ 3SG d§ 


( 2 . 35 ) 

and 



u( £)= ^ . 

3G dg 


(2.36) 


The boundary conditions are now 

r? = 0 : P(g,r]) = 0, t(§,T)) = l , 

on 

71 ^ oo . - 0, T(g,n) -» 0. (2.37) 

an 

It should be noted that the data of any problem appear 
in the transformed equations only through functions K(S) 
and where it is understood that x is a function of § 

defined uniquely by (2.29). These functions play an 
important role in the analysis to follow and they will be 
called the Principal Functions of the problem. 
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5. GORTLER SERIES METHOD OF SOLUTION 


For the purpose of generating series solutions for the 
transformed boundary layer equations (2*33) and (2 #34) the 
forms of K(g) and w(_g) should be obtained in terms of 
power series in g and hence to be determined frcsm equation 
(2*29). Now, we develop the method for long horizontal 
cylinders of rounded-nosed cross-section and a semi- inf inite 
vertical plate for an arbitrarily prescribed surface 
temperature distribution. For either of these configurations’ 
S(x) is known and once G(x) is fixed, we can express K(g) 
and in a series form. 

Yang and Kelleher [11] considered the general wall 
temperature variations for a non- isothermal vertical plate 
in the form: 


G(x) = x" S a. 

k=0 ^ 

where p is any constant and n is the surface temperature 
variation parameter and x is the distance measured from the 
leading edge. The governing boundary layer equations for a 
non- iso thermal vezrtical plate with the above prescribed 
surface temperature variation were solved for CJ = 0.7 with 
n =0, 1/5 and p = 1, 1/2. For our consideration of free 
convection over a non— Isothermal vertical plate, we shall 
assume 


G(x) 


V. 2k 

X 2 a, X 

k=0 


with a^ / O. 




17517 
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^ 3 ^ = 0 for k > 1, it is seen that G(x) varies 
linearly with distance from the leading edge and this 
distritsution gives rise to a similarity solution as shown 
by Sparrow and Gregg [20 ] « This surface condition is of 
special significance as it yields similarity solution for 
the full Navier-Stokes equations governing the flow over 
non— isothermal vertical plate as demonstrated by Reeves and 
Kipponann [ 21 ] • 

Chiang and Kaye X ^ ] considered laminar free convection 
boundary layer flow over a horizontal circular cylinder with 

2 4 

G(x) = 1 + a^x + a.^x • 

Later Koh and Price [ 5 ] solved the above problem for a pair 
of values of a^^ and a 2 • For our studies on non— isothermal 
cylinder and a non— isothermal vertical plate/ we shall take 

G(x) = x-S r a x^^ / a 0 (2.38) 

k=0 ^ 

where 6=1 for vertical plate/ 

and 6 = 0 for round-nosed cylinder. 

Using (2.38) and expressions for S(x) in (2.29)/ (2.35) 
and (2.36)/ we can write £//3(S) and 03(|) as a power 
series in x and then applying the Lagrange^ s formula for 
for reversion of power series [22] # the variable x is 
obtained as a power series in g (see Appendix A). This 
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series for x is used in the ^ and q expansions to 
yield a series for 3 and w in powers of g • The 
coefficients in 3 and W series will ultimately be 
functions of the coefficients in the expansion of G and S, 

To obtain the recjuired series expansions certain operations 
such as multiplication^ division, and inversion^ as well as 
differentiation and integration term~by— term will be performed 
on SG* Such operations require that the series for S 
and G be convergent. Also it should be noted that 
a^ 0 is essential for all the operations to be meaningful. 

Thus K(g) and W('g) for a non- isothermal horizontal 
cylinder are obtained as 

K(g) = K^+ (2,39) 

and 

(ji)(g ) = ( 2 ,40) 

v^cre K = 3/4 and K. and W., j > 1 are the coefficients 

O J J 

to be calculated using the procedure given in Appendix A, 

For non— isothermal vertical plate with G defined as 
in (2 ,.38), we have 

K(g) I + + K3g^/^+ K^|^+ ... (2-.41) 


and 


(0(g) = K(g) - 1/2, 


( 2 . 42 ) 
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where again coefficients J to be calculated 

using the procedure outlined in the Appendix A. 

To solve the equations (2*33) and (2*34)^ we. now 
assume the following series for P(5 yT) ) and T(g/‘>7): 


F(|y77) 

CO 

= S F . (h ) # 

j=o J 

(2.43) 

T(§y77) 

oo 

= S T.(r?) 
j=0 ^ 

(2.44) 

, 


where a = 3/2. 


Substituting (2.39) to (2*44) in equations (2.33) and 
(2.34) and comparing the coefficients of like powers of g y 
we get the following system of ordinary differential equations: 


(a) for vertical plate 

2 

p'" + F p" - p' + T =0 0 

o o o o o 


i- T" + PT'-F'T =0 *, (2.45) 

a o o o o o 


for 


j > 1 

Fr'+ vj« - |(aj + |) 

I TV + F^T'.-(l+ I aj) 


F'F'+(1+ 

P'T.+(1+ 

O j 


4 

3 


4 

3 


otj ) 


F" F 4 K.F' +T.= RP . ^ 
o J 3 j o 3 J 


aj)F^.T;-Fp^. 


.4K.F'T = 
J O o 


RT, 


(2.46 ) 


where 
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KF. = 2 [ i (ak + |) FiF* ,-(1+ T ak) f.fV 


4:' k"j-k 


k j-k 


+ I K 


5 ""j-k .!, ^k-i3 


2^ [(1+ I ak) T^F'__^-(1+ I ak) 


k j-k 


+, 4 K . - 2 F*. T, . ] . 

3-K , „ i JC—l j 


for round- nosed cylinders 


p«/ + F P" ^ K F-' + T 

o o o 3 o o o 


-t 

±, rnf * 4. •p 

a o o o 


= 0 ; 


(2.47) 


for j > 1 


p^'+ F^F'^ I (aj + 2 K^)pJf'. + (1+ I aj) Fq'Fj- | k^fJ +t^.=ff^.^ 


T" + F^T^.- I aj F'Tj+ (1+1 aj)FjTj-4«^jFjT^= RTj 


4 


(2.48) 


where 


* jIL 

Z [ |(K^+ (1+ I ak) P^P”^ 
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and 


RT. = 
3 


j-1 

S 

k=l 


ak (1+ I ak) 


k 

+ . , E 

^ ^ i=0 



The boundary conditions for both the cases are 

at 71 = O 

as 7} -*• oo» (2*49) 

for j > 1 


Fj = F' = T^. = 0 / 

at 71 = 0 


F' - 0/ T. -♦ 0 / 

3*3 ' 

as 71 -♦ oo. 

(2.50) 


The essential feature of this Gbrtler type series 
expansion is that the zeroth— order term satisfies all the 
boundary conditions of the problem and hence already 
represents a solution to the full problem. Succeeding terms 
in the series, only serve to Improve the accuracy of the 
solution inside the boundary layers. These advantages in 
this type of series solution method have been demonstrated 
already in various studies on laminar free convection 
[7ji9, 10/11] .. The zeroth-order equations for vertical plate 
(2.45) comprise the same system as investigated by Sparrow 


= F' = 0/ T = 1 

o o ' o 

O , T 0 

o ' o 
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and Gregg [ 20 ] for non— isothermal vertical plate with linearly 

varying surface temperature distribution* Again# the zeroth- 

order equations (2*47) for the cylindrical body under our 

consideration are the same as that studied by Chiang and 

Kaye [ 4 ] for free convection near the stagnation point on a 

horizontal circular cylinder* The system of equations for 

P.(^) and T.(7?) for 1 > 1 are dependent on the nature of 
J 3 “ 

the driving force through the appearance of the coefficients 
Kj and in the equations (2*46) and (2*48)* Employing 

the following special transformations# these coefficients 
are eliminated from direct consideration in the solutions of 
the resulting differential equations: 

( i) non— isothermal vertical plate and isothermal cylinder: 

Fi = . 

'^1 " ' 

■^2 = S '^22 ' 

^3 = + ‘* 1^32 + ' 

'^3 ~ ^‘'2^31 *'1^32 **■3^33 ' 

^4 = “^^41+ '^"'2^42-^ *T."^3^43+ ^^44+ V4S ' 

T 4 = *1^^42''' *i*'3^43''' *4^45 ’ ^2-51) 


and so on* 
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( ii) non— isothermal cylinder: 


h = +“1^12 ' 

’’i " ■'■“1^12 ’ 

^2 = * V2 2-^“l^2 3 +“2'^24 =^“ 1^25 ' 

’'2 = ^^21 + ‘S222+“^23 +“2^24 + ' V > 1^25 ' ‘2-52> 

and so on. 

The differential equations for ys and Zs together 
with their boundairy conditions are given in Appendix B, 

Now these functions are universal with respect to the special 
data of the problem like the body contour and wall temperature* 
For a specific a, they can be determined once for all and 
then applied to any special case of either of the two classes 
of problems simply by assembling the series for the stream 
and temperature functions* When the series is assembled/ 
then S,G/^ and the expansion for principal functions must be 
introduced for the specific problem on hand. 


6. METHOD OF INTEGRATION 

The non-linear partial differential equations (2*33) 
and (2*34) are essentially reduced by the application of 
series method of solution followed by special transformations 
to a non-linear together with a sequence of linear two-point 
asymptotic boundary value problems. Also# these linear 
Equations have their coefficients determined by functional 


by 

coefficients of the proceeding terms in the series* The 

non-linear eguations are numerically integrated using a 

fourth— order Runga— Kutta integration scheme with a fixed 

step-size* Direct niomerical integration requires that the 

two— point non-linear boundary value problem be converted to 

an initial value problem by estimating two of the five 

required initial conditions* The exact outer boundary 

conditions are satisfied to an arbitrary degree of accuracy 

by integrating to some large value of V and using a modified 

Nachtsheim and Swigert [23] iteration scheme to select new 

values of the initial conditions until convergence is 

achieved* The solutions to the unloiowns F"(0) and 

o 

T'(0) are obtained accurately upto six decimal places* 

The linear boundary value problems are solved by the 
method of complementary functions as given by Roberts and 
Ghipman [24] and we briefly outline this technique here by 
applying to first-order system for isothermal cylinder 
(see Appendix B). We rewrite these equations as follows; 
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Boundary conditions are: 

Yii(O) = YJ^(O) = Z^3_(0) = 0 

and 

Yj^(oo) = = 0 • (2.54) 

Choosing (0,0, 1,0,0)^ and (0,0,0,0,1)^ as initial 
vectors and integrating the homogeneous ecjuations obtained 
by dropping the non-homogeneous term in (2.53), we get two 
linearly independent solutions 

Y^ = (r3_(h), r^Cr?), 

and 

Yg = (s^(‘n), S 2 (^)/ 3^(7)), s^(v), S^ivV 
respectively. Let = (Pi(’l). 

obtained by integrating (2.53) with (Yj^^( 0 ),Y'^( 0 ) , 0 ,Zj_^( O) ,0) 
as initial vector , be the particular solution of (2.53). 

Then Y = Yp + Y^ + Y^ 

is a solution of the linear boundary value problem (2.56) and 
(2.54) where and C 2 are given by 


r n 

=1 

1 


"1 


1^2, 

^11 

1 

r^(V 

! 

P 4 '"e) J 
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Here ri^ is the suitably chosen boundary layer edge* A 
step-size of 1/64 and 77^ =15 or 13 according as the 
Prandtl nxmber is 0,7 or 1*0 for both classes of problems 
under our consideration were chosen to perform the integra- 
tion of the linear equations governing the universal 

functions Ys and Zs* The conditions that both lY^ (tl )i 

j r e 

and IZ. (h )l are less than 10 , were always satisfied 

J 3r © 

for any range of values of j and r considered here. 
Throughout the computation^ calculations were carried out 
upto eight decimal places in a Dec— 10 computer* 

As described earlier/ the solutions of these linear 
boundary value problems are universal with respect to the 
special data of the problon like the body contour and wall 
tcanperature* For cr = 0*7 and 1 *0/ these functions have been 
calculated upto and including the fourth order terms for round- 
nosed isothermal cylinders and non— isothermal vertical plat< 2 / 
for non- isothermal cylinders, they have been calculated upto 
and including second order terms* The missing wall conditions 

P''(0) and T^(0) for zero-order solution and Y"(0) and 
o o 

Z-'CO) for all universal functions are listed for a = 0*7 
and 1*0 in Table 1 for isothermal cylinder, in Table 2 for 
non— isothermal cylinder and in Table 3 for non— isotliermal 
vertical plate* These values are used in the calculations 
of surface heat transf-^^r characteristics of the problems* 

I 
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7 . RESULTS AND DISCUSSION 

In order to denonstrate the usefulness and limitations 
of the present series solution^ the method is applied to 
the following particular cases and the results are compared, 
where possible, with experimental and exact solutions: 

I. Isothermal Surface Condition 

(a) Horizontal circular cylinder: 

We begin by presenting results for isothermal horizontal 

circular cylinders for which experimental and analytical 

results are available for comparison. The experimental data 

obtained by Jodlbauer [ 25 ] is used for comparison with our 

results obtained for the boundary layer velocities and- 

5 

temperatures • For high Grashof values, such as 7 x 10 , we 
can expect the theoretical results agreeing well with 
experimental data* 

The Shape Parameter, k(5) 

First we calculate the shape parameter K(g) for 
isothermal circular cylinder, for which 

o(| ) =0, 

and taking diameter of the cylinder as reference length, 

and S(x) = sin = sin 2x 
where is the angular position of the location x from the 


x_ _ r 

" 2R " 2 
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lower stagnation point and R is the radius of the oylinder* 
Using S(x) in (2.35)^ we find K(0 at various locations on 
the body contour* Now we shall consider the ou-fstion of how 
well the shape parameter K(|) is represented by series 
(2*39) at downstream locations from the stagnation point* 

This consideration leads one to determine the radius of 
convergence of the series for K( ^ ) * S ince only a limited 

nxomber of terms can be calculated, we assess the quality of 
the series by studying what happens to the sum of the series 
when the number of terms included is varied. For this purpose 
more number of coefficients Kj of the series are to be 
evaluated than those are already known [?] • From the general 
computer routine to find the coefficients Kj, we calculate 
K^/ K^/ for this particular case of an isothermal 

horizontal circular cylinder* With diameter as reference 
length, we obtain 

= -0.3265986 , 

= -0.1263288 , = -0.08378184 * 

We observe that these coefficients show a monotonically 
increasing tendency and are bounded above by zero. Table 4 
depicts the 2-term, 3-term, 4-term partial sums of the series 
for , K( I) against actual values for K( ^ ) obtained from 
(2.35) at various locations on the surface beginning from 
the lower stagnation point. A two- term sxam approximates K( §) 




= -0.1955556, 
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very well within 1 percent error upto the downstream location 
^=5 0® from the lower stagnation point. Further downstream, 
it can be seen that a two— term representation for K( ^ ) 
involves increasingly large errors and so the results for 
flow variables and heat transfer calculations obtained by 
Seville and Churchill [?] beyond the location 'V'= 50° down- 
stream carry with them this inaccurate representation for 
K(^), For that matter, the five term expansion for K(g) 
represents it well upto 90° and beyond 100° the 

error in representing K(g ) is large though less compared to 
the two term result for it# The general ability of the 
series representation for K(| ) can be described in terms 
of our findings as follows: on the lower side of the circular 
cylinder, K(g) can be represented well by taking the first 
five terms of the series for it, while on the upper side so 
many additional terms of the series are needed to represent 
I<(§ ) satisfactorily that it may turn out to be Impracticable 
owing to the tedious algebraic process involved in the 
calculation of Kj • 

The Velocity and Temperature Profiles: 

The profiles for velocity and temperature distribution 
in the boundary layer are calculated using the expressions: 
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uL 


V Gr 


T72 


v- i 

3 


g S 


1/3 


P. 




+ S^'’^(f^f^Y'j^(’n)+KjY32(’l)+K3'y33(>l5) 
+ |«(l4Y'^(>l)4-I^K2Y42(n)+K^K3Yj3rt) 

+ K*Y'^(51) + K^YJ^CT?)) +... ] , 


+ g ( K^ K^ (?? ) +K^Z 32 ,(‘n ) +K 3 Z 33 (r? ) ) 

+ g^(K^Z^^(Tl)+i^K2Z^2^^>+^K3243(^) 

+ i4z44(’’)+K^Z45(’^)) +““* 


The velocity and temperature profiles at the location 
=150® are shovm in figures 3 and 4 respectively# It is 
seen that at this location the present results for local 
temperature distribution agree satisfactorily with experimental 
data# In the case of velocity profile^ figure 3 shows a 
regular deviation in the sense that the measured velocities 
are greater than theoretically computed ones# But one 
significant aspect is that the position in the flow where 
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maximum velocity occurs do agree with the theoretical 
prediction as seen from figure 3*. Deviations occuring at 
large distances for velocity measurements from the computed 
result may be due to the measurement of small velocities 
there. It is to be noted that a five-term result for 
velocity profile improves significantly over the correspon- 
ding two-term result around the region of maximum velocity 
occurrence at 'y = 150°, 

Heat transfer: 

In technological applications, the rate of heat transfer 
from the surface to the fluid is of primary importance and it 
can be written in terms of the local Nusselt number as 

[T;(0)+S=/\zf^(0)H-53(Kfz|,(0)+K2Z'2(0)) 

Gr s 

+ |^/2(Kj^K2Z'^(0)+K^Z'2(0)+K3Z^2^°^ ^ 

+ 5 ® ( 0) z'2 ( 0) +K^K3/:' 3 ( 0) 

+ K^Zj^(0)+K^Z'g(0)) +••• ] . 

Taking reference length as the diameter of the cylinder, the 
results for heat transfer at the surface in terms of local 
Nusselt number are given in Table 5 for G = 0.7 • For 
comparison purpose, the finite difference data is also 
Included in this table [is] • From the table we see that the 
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fivG-term series expansion predicts heat— transfer accurately 
upto the location 130® and underestimates it by about 
2 per cent and percent at 140® and 150® respectively. 

Thus the large inaccuracy in representing heat transfer at 
these locations by a two— term expansion [ 7] is greatly 
reduced and thereby, it is shown l±iat G&rtler series can be 
used as far as the location 150®, It should be noted 
that the extensive calculations show though convergence is rapid 
in calculating heat transfer on the lower side of the cylinder, 
more terms are needed to reduce the inaccuracy beyond 
U = 130®, 

Now we compare the mean heat transfer coefficient as 
obtained by the present method with that of e 3 q>erimental 
result. The average Nusselt number (averaged over the 
perimeter of the circle) is 

1 ^ 

Nu = i" Nu(x)dx • 
o 

Taking radius as reference length and using (x,^) transfor- 
mation (2,29) to transform the integral with respect to S 
with the limits of integration given by 0 and ^ (7i) = 2,59, 


we find 
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VJhile a two “term siJin of this series gives for a value 

Gv ' 

0,3027 [ 7 ] , the five-tem sum for it is 0,3069, Thus we 
see Jodlbauer data gives for average heat transfer coefficient 
a value 9 percent higher than a two-term sum and 7 percent 
higher than the five-term sum of the series, 

(b) Horizontal Elliptic Cylinder: 

Now that we have calculated considerable number of 
universal functions associated with the present series method 
of solution/ an attempt is made here to study the applicability 
of this series to elliptic cylinders also. We consider the 
variation of heat transfer with distance from the stagn- 
ation point for a 2;1 elliptic cylinder with major 
axis perpendicular to the direction of gravity (see figure 2)« 
If 2a and 2b are the lengths of major and minor axes 

respectively and e is the eccentricity of the elliptic 

2 b^ 

section of the cylinder then e = 1 — -y. Also x and 
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S(x) = sin (P are given parametrically in terms of.'Ot the 
eccentric angle, by 

X = “ = ! 1-e^ sin^^ dX , 

Cl 

O 


and 


sin 



V 


O 

sin v' 

I 2.2 

l~e sin 


<V' 

L/ 


We then transform the elliptic section into a circle 
and the angular displacement § from the stagnation point 
around the circle in the transformed plane is chosen as the 
independent variable and now it is given by 

I = / 1 1 sin a (1-e^sin^ da. 

o 

The coefficients (j =1,2, 3, 4) for this problcan were 

calculated using our general computer routine for calculation 
of coefficients of the principal function and they are found 
to be 

= -0.8981463, = -0.8955555 ^ 

= -0.6441484 and = -0.08554964. 

It is important to note that in this case a five-terra sum 
for K(g) could represent it well only upto '^= 45°. 

Since these KjS are larger in magnitude in comparison with 
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chose for circular cylinder, we have to do additional work 
in order to extend the validity of the series representation 
-Po-r K(g) beyond 45®* 

The heat transfer coefficient at the surface of elliptic 
cylinder in terms of local Nusselt nximber is given by 


Nu 


Gr 


174 




'/■ 1 2 . 2 ,9- 

1-e sin V 


[t'(0) + 

+ g K2 ( 0) tK^ Z32 ( 0) +K3Z' 3(0)) 

+ 0®(K^Z|^(O)+kJk2ZJ2(O)+K^K3Z'3(O) 

+ kJz|^(0)+K4Z'5(0)) +.•> ] . 

The results of the calculations for heat transfer coefficient, 
taking the length of semi-major axis as reference length/ 
are given in Table 6 for a = 1.0 with the finite difference 
data of r-’crkin [15 ] • It is seen from the table that wjiile 
a two-term sum agrees with finite difference data upto about 
= 68 deg*, the five-term expansion gives sufficiently 
accurate results upto about 115® beyond which the 

agreement between the two data rapidly deteriorates* For 
this blunt orientation the heat transfer coefficient first 
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increases with the distance from the lower stagnation point 
and reaches a maximum at about 92® and then decreases. 

On the other hand, we note that there is reduction in heat 
transfer for circular cylinder as we go downstream locations 
away from the stagnation point [ see figure 5] • 

II Non— isothermal surface conditions 

In great many technological applications the heat 
transferred from the surface is non— isothermal. To treat 
this case, we start with G(x) as given by (2*38), To 
illustrate the application of the series, we have to fix the 
geometry of the body and also the coefficients in wall 
temperature distribution. First we shall deal with the 
simple geometry, namely a vertical plate and then take up 
the horizontal circular cylinder, 

a. Vortical Plate: 

To deal with this case, the universal functions associated 
with the series method are worked out and the missing wall 
values for these functions which are important in the 
calculation of rate of heat transfer at the surface are 
tabulated in Table 1 , For this geometry, 

S(x) = 1 

and 


W(g) = K(g) - 1/2, 



We now consider the wall temperature variation of the form; 


G(x) = sin X# 

Then comparing this with (2*38)/ we get 

% = ^ ^2 = ■ n '■■■" 

With these values of the coefficients in wall temperature 
data, we obtain Kj( for j=l/2/3/4) as follows: 

= -0.1154701 , = -0.02444445/ 

= -0.0055830 and K. = -0.001309096. 

The range of validity of the series representation for 
which describes non- iso thermal condition at the plate/ can be 
discussed on the same lines as we have done for isothermal 
circular Cylinder. The velocity and temperature profiles 
are shown in figures 6 and 7 respectively for x = 0.8/ 1.2 
and 1.6/ x measured from the leading edge of the plate. 

To facilitate comparison of our result for heat transfer 
at the surface with that of the finite difference data [17] / 
we use the following expression for local Nusselt number: 



is the heat transfer coefficient defined by 


where h 
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Now intiroducing the local Grashof number 

3 

g3(t (x) - t J X 

Gr = 21- —22 

X 2^2 

Nu^ may be written as 

X 

1/3 

^ ^74 [Tj(0) + g^/2K^Z'j^(0)+|®(Kjz|j^C0)+K22;'2(0)) 

+ 5®/2(j^j^2^(g)^^2'2(0)+K3Z'3(O)) 

+ |«(!^Z'3(0)+k2k3Z'2(0)+KiK3Z'3(0) 

+ l4^i4(0) + K^Z'5(0)>+...] . 

The results from the series solution and the finite 
difference data are depicted in figure 8 for G = 0,7 and 1 ,0 
in terms of the combined parameter Nu as a function 

X X 

of the dimensionless distance x from the leading edge of 
the plate. It is clear fron the figure that the results 
of our analysis agree well with finite difference data up to 
a downstream distance x = 1,6 from the leading edge of the 
plate though/ of course/ the ability of the series in terms 
of nxamber of terms to be evaluated vary according to the 
location under our consideration. It is found that the two- 
term Siam can give accurate results for the local Nusselt 


Nu__ 

X 



X : 
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number upto x = 0*8 while a three, four and five-term 
partial stum yields results accurately upto x = 1.2/ 1.4, 

1*6 respectively. Also it is clear from the figure 8 
that the range of validity of the series is independent of 
the Prandtl number values considered here. To increase the 
range of validity of the series solution, it is necessary 
to evaluate some more additional terns • But for the range 
of value of x considered here, the series converges very 
well and the results can be considered accurate. Also, from 
the figure it is observed that the effect of the Prandtl 
number increase is to increase the heat transfer at the 
surface, 

b. Horizontal Circular Cylinder 

To illustrate the application of the series to this case, 
we consider some specific values for the coefficients a^ 
in G(x). We choose two distributions of G(x) by taking 

(i) - 1/ a^ = -0.1 and a^ = O for j > 2 

and 

(ii) a^ = 1, a 3 _ = 0,2 and a^. = 0 for J > 2' 

The radius of the cylinder is taken as reference length and 
the temperature at the stagnation point, T^^/is taken as 
reference temperature the above two cases of surface 

temperature distribution, G(x), the coefficients K^. and 0^ 


are 
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(i) 


when 



G(x) = 1-0.1 

-0,1847521, = -0.05442963, 
-0,07698004,^2= -0,02449383? 


(ii) when G(x) = 1+0,2 x^. 


Kj_ = 0,02309401, = -0,05727407, 

= 0,1539601 f ^2 ~ -0-05056790. 


The principal function K(|) is represented by a three- 
term sum of the series for K(§) accurately upto an angular 
distance of 60° in both the cases and is approximated 

well to a distance 50° from the lower stagnation point by the 
first two— term sum of the series for it. 

The velocity and temperature profiles are calculated 
from the esqjressions 

(SG)^/^ + 

V Gir^/ ^ ^ 


and 


^ = T AV) + 53/2(K^z^j_(t7)+ai^Z^2(t)))+g3(K2z2j^(n) 

w °° 
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for a = 0*7 at an angular downstream distance 90® from the 
stagnation point when G(x) = 1— 0*lx^ and are shown in 
figures 9 and 10 respectively* We observe from the figures 
that two- term expansion of the series represents adequately 
well the local velocity and temperature distributions as 
far as 90° downstream from the stagnation point* It is 
only to be anticipated in view of the fact that the dependent 
and independent variables are so defined that as. much infoi>- 
mation as possible is concentrated into the zeroth-order 
solution in the expansion of the series for flow variables# 
while the succeeding terms in the series only serve to 
improve the accuracy of the solution inside the boundary 
layers* 

Also# we give tho velocity and temperature profiles in 
figures 11 and 12 respectively for a = 1 *0 at the location 
60° for both cases G(x) = 1+0, 2x^ and G(x) = l-0*lx^ 
including the isothermal case* An increase in wall tempe- 
rature corresponds to an increase in maximum velocity while 
a decrease in it tends to suppress this value. Also this 
property holds good for the absolute value of the temperature 
slope near the wall* 

The dimensionless heat transfer is computed in terms of local 
Nusselt number from the e 3 g)ression 



77 


Nu 


-( 3 ) 1/4 


(SG)^/^G^/3 

el /4 


[t' ( 0) +g^/2 (Kj_Zj[^ ( 0) +6)^ Z '2 ( O) ) 
+f^K2z'^(0)+KjZ'2(0)+«2z'3(0) 
+0)2 Z'^(0)+Kj_63^Z'g(0) )+..,] . 


The r'esul'ts for the heat transfer are given in figure 13 
for wall temperature distribution G(x) = 1+0. 2x^ for g = 1.0 
and 0.7 and these are compared with those obtained in [s] for 
O’ = 1,0. It is observed that the heat- transfer increases as 
d increases. In figure 14, the heat transfer is compared 
with that of Blasius series solution [sj for the surface 
temperature distribution G(x) = 1-0.1 x^ when d = 1.0. 

We now define a quantity q’^ as follows: 


•J5- 

q = 


[q/(yt„)Gr^''^] 


non— isothermal wall condition 


isothermal wall condition 


# 


where q is the surface heat- flux. This simplifies to 



(Nxi/Gr ^ non— isothermal 


_ 1 (Ku/Gr ^ non— isothermal ^ 

can then be used to give the effect of the non-isothermal 
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wall on heat -transfer. For a = 0,7, fig* 15 shows the 
effect of non— isothermal wall at various locations for the 
two prescribed surface temperatures. For G(x) = l-r0,2x^/ 
the heat transfer at the location 60° from the stagnation 
point is 19 percent higher than that for an isothermal 
cylinder. When G(x) = l-0,lx^/ the heat transfer at 
that location is 13 percent less than that for an isothermal 
cylinder. In other wordS/ the heat transfer is very much 
influenced by the surface temperature variations. 
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Table 1 


Missing wall values of universal functions 
for Isothermal cylinder 



a = 0.7 

(7 s= 1.0 

F^'(O) 

o 

0.859345 

0.817010 

T'(0) 

-0.37023f 

-0.421431 

hqj(O) 

-0.091465 

-0.076768 

l'Zfl(O) 

0.032265 

0.032783: 

Y‘£.(0) 

0.01844? 

0^14393- 

z'i(o) : 

-0.004^15 

i . ^ 

-0^04477 


-0.076577 

—0.064994 


0.O27?63 

0.027981 

|y“(o) 

0.028544 

0. 022518 

( 


-0.007031 

-0.006957 


-0.00340 

-0.002510 


0.000525 

o.ogo5i7 

:v55<o, 

-0.066790 

-0.057186 

|Z'3(0) 

0.023532 

0.024093 

;y”(o) 

0.011210 

0.008936 


-0.002728 

-0.002676 


-0.007388 

-0.005540 
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Table 2" 

Missing wall values of universe functions 
for Non— isothermal cylinder 



a = 0*7 

0 = 1 

Fj'(O) 

j w 

0.859345 

0.817010 

T'(0) 

-0.370234 

w 

-0.421431 

YJJ(O) 

-0.091465 

.J 

-0.076768 

Z'^(O) 

0.032265 

0.032783 

Y-(0) 

-0.183803 

-0.187802 

Z'j(O) 

-0.S80182 

-0.637876 

Y''(0) 

0.018443 ' 

■J \ 

0.014393 

Z|i(0) 

-0.004615 

J 

-0.004477 


-0.076577 

-0.064994 


0.J027363 

0.027881 

Y''(0) 

0.069578 

0.072933 

Z|3(0) 

0.297854 

0.333228 

Y-ICO) 

-0.11670 

-0.120081 

- iJt 

Z'^CO) 

-0.4683187 

-0.5143137 

Y-(0) 

0.029079 

0.026188 

a 


0.020705 

0.018563 
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Table 3 

Missing wall values of universal functions 
for Non— isothermal vertical plate 



or = 0,7 

CT = 1 ,0 

F'(0) 

0,782428 

0.793502 

O 


V 4 

T^(0) 

-0,52807 

-0.595092 

o 


4 

Y|'(0) 

-0,196956 

J 

-0.185822 

Z'^(O) 

-0,356712 

-0.392551 

YJ'(O) 

0,071950 

0.068073 

Z'l(O) 

0,163080 

0.178188 

Y''(0) 

-0,145781 

-0.137417 

Z'j(O) 

-0,307786 

-0.338492 

1 

0,096403 

0,091061 

25 ^( 0 ) 

0.246494 

0,269208 

Y"(0) 

-0.021996 

i 

-0.020957 

4 

Z^jCO) 

-0,069845 

-0.076452 
. 1/ 

Y'|(0) 

i 

-0,117329 

-0.110569 

4 

Z'3(0) 

-0,274678 

f 

-0.302051 

Y'«(0) 

0,033158 

0,031281 

z'i(o) 

0.094^85 

t 0.103264 ' 

1 . j . 


Y''CO) 


0,041588 


0,039532 




a = 0,7 


a = 1,0 
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Table 4 

Shape parameter K(^ ) for isothermal 
horizontal circular cylinder 


0 

10 * 

20 * 

30* 

40* 

50* 

60* 

10 * 

80* 

90* 

100 * 

110 * 

120 * 

130* 

140* 

150* 


0 

0.03655 

0.09192 

0.15729 

0.22972 

0.30741 

0.38902 

0.47343 

0.55965 

0.64677 

0.73388 

0.82012 

0.90453 

0.98614 

1.06383 

1.13626 


K(|) 


0,75 

0.7477 

0.7407 

0.7288 

0.7115 

0.6879 

0,65 7T. 

0.6173 

0.5661 

0.5 

0.4133 

0.2968 

0.1347 

-0.1029 

-0.4794 

-1.1531 


two • 
term 


0.75 

0.7477 

0.7409 

0.7296 

0,7140 

0,6943 

0,6708 

0.6436 

0.6133 

0.5801 

0.5447 

0.5074 

0.4690 

0.4307 

0.3916 

0.3544 


three 

term 

four 

term 

five 

term 

0.75 

0.75 

0.75 

0.7477 

0.7477 

0.7477 

0.7407 

0.7407 

0.7407 

0.7289 

0.7288 

0.7288 

0.7117 

0.7115 

0.7115 

0.6887 

0,6880 

0.6880 

0.6592 

0,6574 

0.6571 

0.6229 

0.6185 

0.6176 

0,5790 

0.5697 

0.5671 

0.5272~ 

0.5094 

0.5033. 

0.4674 

0.4360 

0.4229 

0.3996 

0.3478 

0.3223 

0.3243 

0.2439 

0 ,1980 

0.2426 

0.1240 

-0.0469 

0.1562 

-0,0107 

-0.1321 

0.0675 

i -0.1569 

-0.3372' 
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Table 5 

1/4 

Local heat transfer parameter/ Nu/Gr ' / for 
isothermal horizontal circular cylinder (cr=0.7) 


■O' 

1 

1 Nu/Gr^/'^ __ . 

two 

term 

three 

term 

four 

term 

five finite diff 
term erence [15] 

0* 

0 

0.4402 

0,4402 

0,4402 

0,4402 

0.4402 

10* 

0,03655 

0.4397 

0.4397 

0.4397 

0.4397 

0.4395 

20* 

0.09192 

0,4379 

0.4380 

0,4380 

0,4380 

0,4377 

30* 

0.15729 

0,4350 

0.4350 

0.4350 

0,4350 

0.4348 

40* 

0.22972 

0.4308 

0,4309 

0.4309 

0.4309 

0.4307 

50* 

0,30741 

0.4254 

0,4256 

0,4256 

0.4256 

0,4255 

60* 

0.38902 

0,4187 

0,4191 

0.4192 

0,4192 

0.4190 

70* 

0.47343 

0.4106 

0,4113 

0.4114 

0,4114 

0,4113 

• 

O 

OD 

0,55965 

0.4010 

[ 

' 0.4021 

0,4024 

0,4025 

0,4024 

90* 

1 

' 0.64677 

0.3899 

0.3915 

0.3920 

0,3922~ 

0.3922 

100* 

0.73338 

0,3770 

0,3793 

0.3801 

1 0,3805 

0,3806 

110* 

0.82012 

0,3621 

0,3652 

0,3666 

0.3672' 

0.3677 

120* 

0.90453 

0.3450 

0,3489 

0.3509 

0.3520 

; 0.3532' 

130* 

0,98614 

0.3252 

0.3300 

0.3327 

0.33-44 

0,3370 

140* 

1 ,06383 

0.3019 

0.3075 

0,3111 

I 0.3135 

0.3190 

150* 

1 ,13626 

i. 

0,2740 

0.2801 

0.2845 

; 0.2877 

0.2986 


Notes The diameter of the cylinder is used as reference length. 
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Table 6 

1/4 

Local heat transfer parameter/ Nu/Gr ' s for 
isothermal/ horizontal blunt 2:1 elliptic 
cyl inder ( a =1 . 0 ) 


/ 

g 

K(|) 

Nu/Gr^/^ _ 

two 

term 

five 

term 


0.0 

0 

0.75 

0.3542 

3542 

0.3542 

0.2 

0.069284 

0 ,7333 

0.3555 

0.3555 

0.3555 

0,4 

0.172049 

0,6811 

0.3589 

0.3590 

0.3593 

0.6 

0.288335 

0.5871 

0.3648 

0,3652 

0,3657 

0.8 

0.409145 

0.4436 

0,3733 

0.3745 

0.3747 

1.0 

0.528110 

0.2544 

0.3841 

0.3872 

0,3861 

1 .2 

0.641007 

0,0829 

0,3954 

0.4021 

0,3984 

1 .4 

0.746501 

0,1340 

0.4016 

0.4142 

0.4081 

1.6 

0.847175 

0.5821 

0.39 37 1 

1 

0.4146 

0.4088 

- 

1 .8 

0.948753 

1.0573 

0,3680 ^ 

0.3997 

0.3958 

2.0 

1.056329 

1,1553 

0.3318 

0.3778 

0,3713 

2.2 

1,171286 

0.9443 

0,2934 

0,3585 

0.3407 

2.4 

1.291330 

0.5313 

0.2559 

0.3454 

0.3081 

2.6 

1,411567 

-0.1646 ! 

j 

0,2190 

0.3365 

0,2752 

2.8 

1*525005 

-1,6801 . 

0.1801 

0.3212 

0,2418 

3.0 

1.621315 

••3 » 44:1. 4: 

0.1301 

0.2687 

0.2056 





Notes The length of semi major-axis of the elliptic 
cylinder is used as reference length. 
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APPENDIX - A 


Eixpressions for K(g) an d in terms of § 

Recall the expression (2.29) for i / that iS/ 


I = /^ (SG)^/^dx . 

o 

We assume the following power series for SG: 


SG = S^x + S^x^ + SgX^ + S^x"^ +••• (Al) 

with / 0, This is true for either vertical plate or 
cylindrical body with necessary changes in the value of the 
coefficients . 

We write P = SG 

By Lagrange^ s foimula for reversion of power series ^22 3 / 
we get 

X = ^ (A2) 

}c=l ^ 

with coefficients <Xy, given by 

= CA3) 


vrhere 


H^(t) 


hht) ' 


(A4) 


with 

h(t) = s^+ S^t^+... (A5) 



88 


I can now be written as 

§ = / p ^/3 Ip (^p 

as P is analytic and P'(0) O. 
Using (A2) in (A6), we get 


(A6) 


/ pl/3 ( 2 k OL P^'"^)dP 

o k=l 


On evaluating the above integral , g is obtained in the 
following series form 

= I ^ otgP^ + H 

On cubing both sides and simplifying/ it becomes 


3 S 7 9 

- = a. P + a P + a-P + a.P +•.. 

p vj,p ^2 4 


Further simplification leads to 

|V4 ^ 1/4 p (1 g(p))l/4 . 


(AS) 


whore 


• • • 


g(p) = a 2 P^ + a 3 P^+, 

E 3 (panding (A8)/ we arrive at the form 

= bj_P + b^P^ + b2P^+*.« . (A9) 

|p £-l - bj^ {AO), we apply Lagrange's 


Using the fact that 



89 


formula to equation (A9) to obtain 




“h# • « 


(AlO) 


where are to be determined in the same way as 

Substituting for P in (A2) and on rearranging the coefficients 
of like powers of g , we get 




(All) 


We recall the expressions (2*35) and (2 .36) ’which can be 
written as 


=1 + ;;l73S ' 


(A12) 


and 


Q(|) = 


dG 


3P' 


J73 35c • 


(A13) 


Using (AlO) in (A7) and (All) in the derivative of (Al), 


3P' 


173 


and 


dP 

35E 


can be written as power series in § 




3P 


173 


• • • / 


(A14) 


and 


~ = f. + f„ + f.,§^+... 

dx 12=“ 3^ 


(A15) 


With (A14) and (A15), (AI 2 ) gives K(g ) in the following form 
K(t) = K^+ 


for both vertical plate and round-nosed cylinders. 
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Again^ expressing in terms of S (using (All)) 

(3.G 

and substituting in the expression for # we obtain a power 
series in | for Using this expression for ^ in 

powers of ^ , (AlO) and (A14) in (A13)/ o(|) can be written 
as 

(0(g) = (0j^|V2 + +... . (A17) 

For non— isothermal vertical plate/ we have 

0)(§) = K(g ) - 1/2' (see 2 #42) 

which with the help of (A16) becomes the following series 
in ^ 

0(§) = (K^-1/2') + • 


(A18) 



Non- Isothermal Vertical Plate* 


Using the transformations (2*51) in equations (2*46) 
and scaling out coefficients from the resulting equations^ 

we get the following set of differential equations for 


universal functions: 


Li(Yii,Zii) - ^ » 




^1^^22'^22^ “ 3 ^o * 


^2 ^^22 '^22^ 


4 p' ; 

o o 


^2 ^^21 '^21^ 


= 3 3Y,^Y- + I P'YJ, . 




A ^ 

I F' » 

3 o 


^2^^33'^33^ 
^1^^32'^32^ 
^2 ^^32 '^32^ 


= 4FJTo I 

= 3Zj3_Y'i-3Y^^Z'i+ SYj^Zji-SZfiYjj^ I 




= 8^U^22-3^iU22-5’'ii''22'" I ' 

+ 4[f'(Z^^+Z22)+Tq^‘^11'^'^22^ ^ 


^2 ^^31 '^31^ 



92 


Li(y^ 2 ,z^g) - I Fq / 






4 p/T , 

o o ' 


lOYj <-^'21-'^2l'^n > 


+ I ^^o'^22*ni^n'> ' 


+7 ( YJ 1 Z33-ZJ jY 33 )+4 ( Pjzls+Yj ^ Zj^ 4Z^ ^ Y'^ +T„Y^; 
i(Y^ 3 ,Z 43 ) = 10Y£4Y'3-3Y4jY''+ |p' CYJ^+Y' 3)-7Y” Y 33 . 


'^ 2 '\ 3 ' 243 > = 2Cz^iY'3-Y4^Z'3)+7(YfiZ33-Zi^Y33) 


h'^ 42 '^ 42 ' 


S*^ 42 '^ 42 ^ 


h(Y 4 l'^ 41 > 


+^<^i 233 +FjZ 4 ^+T^Y' 3 +T^Yi^) ; 


lOYjiYji-SY.^Y^i+lOY'jYl^-TYf^Yji 

-5 <^ 21^2 2 « 2 l ^22 > + I^^P'Y'.+Y^J ); 

+ l<^ix’" 22 +^;’' 3 l’'' 

^11^31 “%1^31^''^^^21^22'‘‘^22^21~^22^21“^?L^ 
+ 7 C^{ "> *^'-^o^ 21'*'^11^11 '*''^o ^21 


■•■^ 0^31 ■'■^1 1 ^2 2 1^22 ^ ® 

^^^22~’^22^2p ■'' 3 *' 0^22 ' 


^2 ^^41 ^^ 41 ^ 


5 ( 222 ^ 22 -'' 22 ^ 22 >-^‘^i 222 +VM>' ' 
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where 


Ys' and Zs* 




L^CYj 


Z 

J 


satisfy 

= Y'"+P^Y''- |(aj+ |)p'Y'^+(l+ I aj)p''Yj^+z.. 

>=H- 

+(1+ I aj)T-Yj_ -T^Y'._ • 


The boundary conditions are 


Y. (O) = Y' (0) = o = Z. (0) , 

J- J~ J- 

Y'_(<») = 0 = Zj_(oo) . 


Iso-Thermal Cylinder: 


In this case/ 6)^ = 0 for all j in (2*48)* Using the 
transformations (2.51) in equations (2.48) and scaling out the 
coefficients Kj s from the resulting equations, we get the 
following sets of differential equations for universal functions: 


L, CYi.,Z,,) 


4 -r.#' 




^ 2 ^^ 21 ^^ 21 ^ 
^ 1 ^^ 22 '^ 22 ^ 


0 


= I (K^+a)Y'j-(l+ |a)Yj^jYq+ I F^Y'j_ 




4 /^ 
3 o 


^ 2 ^^ 22 '^ 22 ^ = ° ’ 
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|(3a+2K_^)y^'^Y'2-(l+ |a)Y^^Y''-(H- |cc)Y''Yj2 


^2 ^^31 '^31^ 


l''i<^n'^ ^|2> • 

l“^11^22-‘"-^ I “"^11^-22- 


-(1+ lajYjjZfi } 


^l(^32'^32> 


^2‘'^32'^3P = 


L^(Y33,Z33) 


L2(Y33.Z33) 




’^2‘'^41'^41> = 
'^^^42'^42^ ” 


|(3a+2K^>Y'jY'i-(l+ I a)Y^j_Y'{-(l4 fcDY^Y^i 

+ I ^o^21''' 3^11.' 

I aZj^j^Y'j^-(l+ |a)Yj^jZj3+ 

-(1+ |a) zJ^Yji , 

3'^0 * 


0 t 


^2 ^^42 '^42'^ 


|(2a+K^)Y'2-Cl+ ^'^'^22^22'^ 3^0^22'* 

l^^22’^22“^^‘^ §^^^22^22 » 

I ( 4a+2 K^ ) Y^ - ( 1 + |a ) Yj^ ^Y^l - ( 1 +4a ) Y^ J Y 33 _ 

+ |(2F'Y'i+Y'i+2F;y^3^+2Y|^Y'^): , 

I aZ3^3^Y'^+ 4aZ33_Y{i~(i+ |a)YiiZ^i-'(i+4a)z£j^ 

+ |^^^22^2l-^^2l222'^‘"^^"‘ l^^^^2l^22-*-^2Y^2l^ ? 
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^'^43'243> = l'4a+2K^)Y{^Y'3-(l+ |a) YjjY''-(l+4a)Y''y33 

^2‘^43'243> = I ®2u^ 33+ “’fil233-<l+ l“>^U^33 

-(l+4a)Y23Z|ji_ ; 

~ §‘^^11^32'" 

+ |(2a+K^)Y|2-(l+ |a>Y22^2^l‘''o^32+5^l'l^22>. 
I'2(^44'244> = I“2u^32-^?^11^32-'1-^ I“>^1^32 

-(l+4a)2'^Y33+ faZ^jY-j-dt §Ct)Yj22|2 , 

4<^45'^45’ ' 

4<^45'245) = ° 

where 

'^j J = T.-^V”.- I‘“:!+2 V^o^L 

+ (1+ |aj) F”Yj_+ Zj_ ^ 

I.2(Yj_/Z^_) = ? Zj'.+Vj.- I • 

The boundary conditions are 

Y. (O) = Y' (O) = Z. (0) = 0/ 

3 * J - Jr 

Y'^ («) = 


0 =S Zj^(o°) 
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Non- Isothermal Cylinder: 

Using the transformations (2*52) in equations (2.48) and 
scaling out coefficients s and s from the resulting 

equations^ we get the following sets of differential equations 
for universal functions: 






= 0 i 


q(yj2-Zi2> 

'^2 '^21 '^21^ 
’^1^^22'^22^ 


= 4F' » 

= |(K^+a)Y£j-(l+ |a)Y^iYjf+ I PjYji / 

= |aZllYii-(l+ > 

/> 2 

= 1^; ' 


'^2^''22'^22^ 


= 0 t 




^2^^23'^23^ = |^2i2^12“^^'^ |a)Yi2Zi2+^^^0^12'*'Vl2 ^ ’ 


^i^^24'^24^ - 0/ 


^2'^^24*^24^ 


4P'T t 
o o 


= |■(a+K )Y/iY/2“^1+ ^^11^12'*’^11'^12^'‘' 3^0^12 


^o'ni"i2 
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where 

'Zj - ) = ^5" +Vjl - 

+ (1+ |a3)F^'Y. +2. . f 

•3 O J •• J «| 

h'h-'h-’ 'H'--" Vj.- I 


The boundary conditions are 


Yj (0) 

Y' («) 
3 *“ 


Y' ^ (0) =- Z^ ^ (0) = 0 * 

Zj ^ (oo) = 0« 
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Fig. 1 Physical model and coordinate system 
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1-5 



► T} 


Fig. 4 Temperature distribution around a heated 
isotnernrKil circular cylinder;'2>=150'’ 
and cT'rO'?. 



Fig. 5 Heat transfer over isothermal circular 
and elliptical cylinders; cr=1 -0 . 








-ig. 9 u-vetocity distribution around a 
non -isothermal circular cylinder: 
G (x) = 1 -0-1 and o~ = 0-7. 





j{x) = 1-0-1 
s other mat 


Fig. 11 u-velocity distribution around a 

non-isothermal circular cylinder 
2> = 60° and cr = 1-0. 



G(x) = U0-2x2 
G(x)^l-0-1 x^ 
isothermal case 


Fig. 12 Temperature distribution around a 
non - isothermal circular cylinder', 
=60° and cr=1’0. 



— Three term sum(G6rtler) 

O Three term sum (Blaslus) L5J 

o 



Heat transfer variation over a non 
isothermal circular cylinder; 

G(x) =1-0 + 0-2 x2. 



3-t<2rm Sijrn ( Gort ler series) 

O 3 -term =^um ( Blast us scries) [5] 
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CHAPTER III 


LOCAL NON-SIM ILARITY SOLUTION OF FREE CONVECTION 
FLOW AND HEAT TRANSFER FROM AN ISOTHERMAL PLATE 
INCLINED AT A SMALL ANGLE TO THE HORIZONTAL 
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1* INTRODUCTION 

In this chapter^ we study analytically the laminar free 
convection boundary layer flow over a flat surface that is 
inclined at a small angle to the horizontal. Surf ace inclination 
results in a pressure gradient across the boundary layer 
(being induced by the surface normal component of the buoyancy 
force) and leads to a theoretical analysis more complicated 
than that for vertical surfaces* In the case of heated vertical 
surfaces maintained at a temperature greater than the surrounding 
fluid, a local region of lower density produces an up\> 7 ard 
buoyancy force and a boundary layer is formed over the body 
which is driven directly by this force* But if the surface 
is horizontal or nearly so, the buoyancy force along the 
surface is very small and a very different kind of flow drive 
arises* In the particular case of a horizontal surface, 
there is no component of the buoyancy force along the surface 
and there is no direct drive of any tangential flow that may 
result. Rotem and Claassen [i] carried out a flow visualiza- 
tion experiment using a semifocussing S chi ierei system which 
clearly indicated the existence of a boundary layer near the 
leading edge on the upperside of the heated horizontal 
surface though its horizontal extent is limited* Stewartson 
[ 2 ] was the first to give a theoretical description of a 
boundary layer over horizontal surfaces under the action of 
buoyancy force. Since there is no component of buoyancy 
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along the surface^ the accelerating flow must be driven 
indirectly by a buoyancy - induced pressure gradient. 
Stewartson's analysis contains a sign error in the pressure 
gradient term of the horizontal momentum equation and thus 
condition for existence of such a boundary layer is that the 
heated plate is to face downwards. Later/ this error was 
corrected by Gill/ Zeh and del-casal [ 3 ] and one finds a 
horizontal boundary layer above the heated plate driven by a 
negative pressure gradient. This is called "indirect — 
drive ' * and is due to the pressure deficiency in the thermal 
transport region next to the surface with respect to the 
ambient medium. It is indirect in that the buoyancy force 
does not drive the flow. Rather/ it produces a pressure 
field which induces such a flow. 

For a heated plate inclined at a small positive angle a 
to the horizontal/ the boundary layer flow is driven not only 
by the above-mentioned indirect mechanism but also by a 
component of the buoyancy force along the surface. 

Jones [ 4 ] made an analysis of this type of flow by considering 
that except in the neighbourhood of the leading edge/ the 
direct drive due to dominates the induced pressure 

effects/ which can therefore be neglected. The resulting 
flow far downstream in his analysis is then vAiat is described 
by the classical free convection solution along a vertical 
plate. Pera and Gebhart [s] considered the free convection 
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flow over a horizontal and slightly inclined surfaces. The 
effect of inclination on the flow and transport is analysed 
by perturbing flow over horizontal surface. Their study 
includes experimental data for the tanperature distributions 
and heat transfer parameter* When the inclination of the 
plate is negative, the pressure gradient associated with 
the indirect- drive is favourable for motion while the 
component of the buoyancy force along the surface is opposing 
the motion. Jones [4 ] dononstrated that in this case the 
flow eventually separates though, of course, separation is 
not accompained by the development of a singularity in the 
solution of the governing equations* 

In this chapter, the boundary layer in free convection 
flow over the upper surface of a heated semi— infinite 
horizontal plate inclined at a small angle a with the 
horizontal is studied analytically by the local non-siinilarity 
method of solution [^6 ] * The range of ot considered here, as 
obtained from the analysis, is defined by the requirement 
that A , the limit of Gr^^^ tan a as Gr - «>, be 0(1), 
where Gr is a suitably defined Grashof number* The conser- 
vation equations are transformed such that they can lend 
themselves to local non— similarity solutions. In this 
method of solution, all the terms in the transformed 
conservation equations are retained and terms are selectively 
neglected only in the subsidiary equations. This solution 
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method, in general, offers to non— similar convective heat 
transfer problems many of the advantages of similarity 
analysis while at the same time, providing an assessment of 
accuracy of the results obtained within the method itself. 

The forms of the local non— similarity method considered 
here are the two and tnree -equation models obtained by 
retaining the transformed conservation equations as such and 
only selectively neglecting terms in the derived subsidiary 
equations at the second and third level of the successive 
approximations in this method of solution. The corresponding 
two-point non-linear boundary value problems for these two 
models at a streamwise location are solved by the initial 
value technique as modified by Nachtsheim and Sv^igert [6] • 
Particularly, the application of this numerical procedure 
for the highly non-linear three-equation model was ccmmenced 
at some downstream location with good initial estimates of 
missing wall conditions generated by quasilinearization 
technique [7] . Numerical results for the local surface heat- 
transfer, wall shear-stress and velocity and temperature 
distributions are piiesented for Prandtl number O.TT and 1.0 
when A = 1 and -1 • 

For small positive angle of inclination, our study makes 
no assumption on the neglect of induced pressure effects frcxn 
scrnie distance downstream away from the leading edge of the 
plate and the method of analysis adopted here facilitates 
solving for the entire region of flow. 


On the other hand. 
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J ones [ 4 ] had covered the whole region by working with two 
sets of transformed boundary layer equations in order to 
give non— similar corrections which bridge the gap betx-jeen 
his series solutions, one valid near the leading edge and 
the other far downstream of the flow* It is shown that our 
results for transport rates at the surface and also the 
velocity and temperature distributions compares well with 
that of Jones in the entire flow region* When the angle 
of inclination is negative, our method of solution predicts 
the point of separation of the flow and notably, for the 
Case of air, the two— equation model locates this point 

-I 

within an error of percent from the exact value as 
reported by Jones [4 ] • This stands in sharp contrast to 
its behaviour in the case of forced flow over a circular 
cylinder, where this technique yields large deviation in the 
boundary layer characteristics of the flow while approaching 
the region near separation point* Perhaps, this behaviour 
may be due to the fact that there is no evidence of. singula- 
rity at the separation point, unlike in forced flow, where 
a classical singular behaviour is to be expected in skin- 
friction factor near the point of separation* 

2* BASIC EQUATIONS 

We consider the steady two-dimensional free convection 
flow over a semi-infinite flat plate inclined at an angle a 
to the horizontal placed in an infinite expanse of fluid. 
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The plate is maintained at a tanperature which is 

different from the ambient temperature of the fluid. 

Ihe Boussinesq form of the governing equations for the flow 
as derived from equations ( 1 .9 ) / (1.13 ) and (l.ll) in 
cartesian coordinates x and y, along and normal to the 
plate respectively, with origin at the leading edge may be 
written as t 


^ = 0, (3.1) 

3x ay 


u 25 + v 2 i = -i 2 l + y (a! 5 _ + + gg sin a (T-T„) . 

3 x 3y ^ dx ax"^ dy 

(3.2) 

+ = + ^ cos a (T-Tj 

ax 3y ^ ay ay 

(3.3) 


- aT . 
u — + 

V = a 

9 2 

,a^T . a^T. 
-2 ^ -2r 

ax 

ay 

ax ay 


(3.4) 


Here quantities u and v are the components of velocity 
in the directions along and normal to the plate respectively, 
P and T are respectively the pressure and temperature 
fields induced by motion, v and a are the appropriate 
transport coefficients, P the density and a the angle of 
inclination of the plate with the horizontal. 


The variables are now non-dimensional ized and stretched 


in the following ways 
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X 


p = 


2 . V - 

y.„ 

* u 

_ \ 

L ' y “ 

LGr*’^'^^ 

% 

— 2 

PL^ 
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(3.5 ) 


V 7 here L is a characteristic length along the surface and 
Gr is the Grashof niimber given by 


Gr = 


g3(T^- Too) cos a. 


(3.6) 


Using(3.5), equations (3.1) to (3.4) are reduced to 


+ az ^ 0 , 

3x 9 Y 


(3.7) 


u fi + V 
9x 9y 


u + V az 
ax ^ ay 


u + V ^ 

9x ay 


_ 2E + £“ + if" Gr~^^^+ Gr^^^ta.n a e , C3.8) 

3^ 37= 3x2 

_ az Gr2/5+ eGr2/S+ a!z + Gr'^/S ^ , t3.9) 


9Y 


a y 


dyt 


1 [ 2!| + Gr-2/S a!|] 

^ ay^ ay*^ 


(3.10) 


where 0 is the Prandtl nvimber of the fluid. 

The boundary conditions in terms of non-dimensional 
variables are 


u = 0 = V and 0=1 on y = O for x > 0^ 
u,e,p ■* 0 as y -» oo (3.11) 
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and 

u=p=0=o on x = 0 for y>0* (3.12) 

By ignoring terms which are 0(Gr"*^'^^) relative to 
tnose retained in the limit Gr -♦ we get the boundary 
layer equations 


3u 

ax 



(3.13) 


u 


ax 


ay ax ^ gy2 ' 


(3.14) 


0 = - aE + e 


(3.15) 


u 


as + =i 2-2 

3X ^ iY a gy2 


(3.16) 


where the range of angle of Inclination a considered here 
is defined by the requirement that 

A = Lt Gr^^^ tan a 
Gr* ^ oo 

be 0(1) so that the buoyancy force term is formally comparable 
with the induced pressure gradient along the plate* For A = O^ 
the flow is over a horizontal surface and for A -* «> we will 
get the classical free convection problem for vdiich the 
scalings used here became inappropriate. 

We now choose the stream function which satisfies 
the continuity equation (3.13) and now the boundary layer 
equations become 
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^ W ~ w w 

y xy X * yy 


-P^ + f 


yyy 


+ A e 


# 



^y X X y a °yy ' 


with boundary conditions j 

^ = 0 = 0=1 on y = 0 for x > 0 , 

-*0 as y ^ oo 


and 

Wy = p= 0= O on x = 0 for y > O , 

3. TRANSFORMATION OF BOUNDARY LAYER EQUATIONS 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


(3.21) 


We now introduce the following transformations 


T? = y X 


'2/5 


(3.22) 


S = F(x/r)) , 

p = x^'^^ G(x^n) * 

0 = HCx^-n) (3.23) 

to transform the boundary layer ecjuations (3.17) to (3.19) 
into non-similar form: 

5F'''+3PP'* -F* + 2r?G^-2G+5Ax^/^H=5x(|| + |~ -F" || ) 

(3.24) 

G' = H (3.25) 
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SH'-' + 3aFH' = 5x(F' H - H' |-|) • (3.26) 

In the foregoing equations^ the primes denote partial 
differentiation with respect to n .* Wie boundary- 
conditions are^ now: 

F(x,0) = F'(x,0) = 0, H(x^O) = 1 , 

F'(x,oo)= G(x*~) = H(x,«) = 0 . (3.27) 

Here T) would be a similarity variable if A s= 0/ as 
in that case the boundary layer over horizontal isothermal 
semi-infinite plate is similar. Otherwise^ the -transformed 
momen-tum and energy equations remain a system of partial 
differential equations wi-th — terms on -the right hand 
side providing -the major obs-tacle to -the solution. We say 
this set of equations is in non— similar form, -the non— similarity 
arising as a consequence to A 0, 

local non-similarity model 

To correct the drawbacks of the local similarity model/ 
Sparrow# Quack and Boemer [8 ] # described a local non- 
similarity me-thod for ob-taining solutions of the non-similar 
boundary layer equations. This preserves -the most attractive 
features of the local similari-ty model (-that is, quasi odes 
and locally Independent solutions)# while retaining all of 
the terms in the conservation eq^uations# The 
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departure from local similarity is accounted for through a 
sequence of successive approximations as explained in [ 8 ] . 
Only terms in the equations subsidiary to the conservation 
equations are selectively neglected* This technique is 
particularly useful for solving problems in vdiich many 
unknowns appear and for which forward marching technique 
needs considerable conputation procedure in knowing the 
solution at downstream locations* Moreover/ in this method, 
manipulation to introduce a finite difference expression is 
not necessary and the derived ordinary differential equations 
can be directly integrated n'umerically* The formulation of 
the system of equations for the local non— similarity models 
will proceed as follows* Taking 

Q = A x^^^ (3*28) 

and 

f = 21 and g = li , (3.29) 

3x 

the boundary layer equations and boundary conditions 
can be written as 


1 FF"- J- F'^+( Q+ i “G = x(g+F'f'-l’'"f) » 

5 5 3 ^ 

(3.-30) 

G*"+ J' OFG'* = ax(P'g' — G'^f) / (3.31) 

5 

F(x/0) = F'(x/0) = 0 G'(x/0) = 1 


and 


F'(X/<»)= G(x,~) = G'(x,°°) = 0 


(3*32) 
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Differentiating the above set of equations and boundary 
conditions with respect to x and taking = A we get 

£-»+ I M-. 7 I F-.£ _ I g+(|„+Q)g-+ 3 

= X 1 ^ (g+F'f'-F'^f) , (3.33) 

g'"+ I aFg^' + I aG"f - aF'g' = fJx (F'g'-G'"f) , (3.34) 

f(x^O) = f'Cx/O) = g'(x,0) =0 

and 

f'(x,oo) = g(x,oo) = g'(x/<») =0 (3.35) 

P,f,G and g appear in the equations (3,30) to (3.35). If 

not for the presence of the term x ^(g+F'f'-F'’'f ) and 
3 

dx -^(F^g' — G*'f)/ these equations could be treated as 
9X 

ordinary differential equations and solved locally. To 
provide the desired local autonomy to equations (3.30) ^ (3.31), 
(3*33) and (3.34), it is now assumed that the terms 
X -- (g+F'f'-.F"f) and cJx ■— (F'g'-G^'f) can be dropped 
from these equations. This reduction is readily justifiable 
for x?“values that are close to 0, whereas when x is not 
small, we can Invoke the postulate that the dropped ones 
are small. 

The philosophy behind the reduction is similar to that 
employed in deriving the local similarity model, but with a 
fundamental difference in the outceme* 


In the case of the 
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local similarity/ it is postulated that for x not small/ 
g fF"+ p'f' and - G"f are small so that right 

hand side of (3.30) and (3,31) can bo approximated by zero* 
Thus the approximation introduced leaves the boundary layer 
equations truncated. Whereas in the non-similarity model/ 
the eguations (3,30) and (3,31) are left intact and the 
approximation is introduced only in the subsidiary equations 
(3,33) and (3,34), Thus we can anticipate that local non- 
similarity method will yield accurate results for boundary 
layer characteristics than the local similarity model can give. 

Thus the governing equations for local non— similarity 
model are now given by 

(i) Equations (3,30)/ (3,31) and boundary conditions (3,32) 
and 

( ii) approximating equations 

I Ff'i- I P'f'+ Ip^'^f- |g +(|77+Q)g'+- | Qj_g' 0 , 

I cjpg" + I aG''f-aF'g^ = o 

and boundary conditions ( 3 • 35 ) • 

For a given ^ Q and would be known at any 
streamwise position x# Therefore^ at a fixed streanawise 
positioner x#Qand itiay be regarded as known constants 
and equations (3i» 30 (3# 31 (3 #33^ and (3«344 can be 

treated as a system of ordinary differential equations p 
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Solutions of such equations are straightforward by techniques 
regularly employed for similarity boundary layer equations. 
Equations (3.30)^ (3.31),(3.33fi)and (3.34&.)with boundary conditior 
(3.32) and (3.35) constitute what is known as the two- 
equation local non-similarity model for our free convection 
problem. It represents the first stage in a succession of 
locally applicable multi- equation systems which are expected 
to provide increasingly accurate results. 

We now derive the Three— equation model for the problem. 

Put ^ “ lx ^ ~ lx (3.36) 

so that 

-/ 2 

X 9- (g + F^f^-F-^'f) = x(g - fF"+f F'-ff"+f" ) 

3x 

and 

qx |- iF'g'-fG'*) = Ox if'g*- fg"+ F'g'-fG") 


Then differentiating (3.33)/ (3.34) and (3.35) with respect 

—7 /5 

to X and taking Q 2 = A ^ z we obtain 

f"+ |2. p-'?'+ 11 F"f - i + (|i +Q)g'+ I 

. |_ o^G'^ iSff”- f''= X . 

i (3.37) 

I aFg''+ afG''-20(F'g"+f^gO + afg" 

= ax |-(f'g'-fg''+F'g'-fG'0 


(3.38) 
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f(x,0) = f'(x/0) = g'(x/0) = O# 

f'(x,oo)= g(x,oo) = g,(x,oo) =0. (3.39) 

To make the system locally autonomous^ we drop right hand 
sides of (3.36) and (3.37). On the other hand, all the 
terms in the transformed boundary layer equations (3.30) 
and (3.31) and in the first subsidiary equations (3.33) and 
(3.34) are retained. Since approximation is made in the 
secondarily subsidiary equations (3.37) and (3.38)/ the 
boundary layer characteristics should be even more accurate 
than those of the two- equation model. 

Thus the three- equation model comprises of the following 
three sets of equations and their boundary conditions: 

(i) equations (3.30) and (3.31) with boundary conditions (3.32) 

(ii) equations (3.33) and (3.34) with boundary conditions (3.35) 

and ( iii) the approximating equations 

1 '"+ P'f'+ ~ P'' f - g + (|n +Q)g»+ I 

2 

. |_ Q G' + ~ ff** - P f' =0 (3.40) 

25 2 5 5 

and 

g'"^+ I gpg" + crfG" - 2CTP^g' -2af'g'+ afg" = o 

(3.41) 

with boundary conditions (3.39). 



115 


5 • NUMERICAL SOLUTIONS 

We have already seen that the application of non-similarity 
raethod of solution to this problem^ the two and three -equation 
models give rise to four and six coupled equations respectively 
to be solved in each case as non-linear simultaneous equations. 
To initiate forward integration of the sets of equations in 
these two models requires N nvimorical values for the functions 
and their derivatives at h = 0 where N is the total order 
of the system of equations in each case. Thus, for example, 
for the three-equation model, a total of eighteen function 
and derivative values are needed at T? = 0 but only nine of 
than are available from equations (3.32), (3.35) and (3.39). 
Therefore it was necessary to determine simultaneously the 
nine missing starting values. We are not aware of previously 
published boundary value problems where it has been solved 
for as many as nine simultaneously unknown starting values. 

Each of the systems of equations for the above forms of the 
local non- similarity method was solved by employing the 
Runga-Kutta integration scheme in conjunction with Newton- 
Raphson shooting method as modified by Nachtsheim and 
Swigert [6] so as to fulfill the-lDoundary conditions at the 
edge of the boundary layer . By solving directly the 
original equations in these models through the initial value 
technique, we retain the main feature of local non— similarity 
method that no approximation is made in the conservation 
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equations. This is in sharp contrast to the numerical method 
of solution^ namely the selected points method of LanczoS/ 
used by Jones where it is necessary to go in for linearized 
form of the conservation equations at each location of the 
downstream flow. Also/ unlike his numerical method of 
solution where the downstream region is to be covered in 
steps^ these local non-similarity forms offer scope to solve 
the problem at any downstream point independently of the 
other when the estimates of the missing wall conditions 
there/ are available* But/ because of the highly non-linear 
character of the sets of equations in the two forms of the 
local non-similarity method/ it was essential to have good 
initial estimates of the missing wall values to enable us 
to integrate the sets of equations in these two models. 

We obtained good estimates of missing wall conditions for 
the two-equation model by considering the reduced linear 
boundary value problem obtained after dropping off non-linear 
terms in this model. The solution of the linearized problem 
thus obtained is 


^(”5 = 30oV ‘2 - I 




G(ri) = i 

^ 2 % 
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f(^) - Xo ^ 1^00 (2r] - 3 ^ T}2 + 1 r)3) , 

oo 

g(v) = 0 , 

where rj^ is an arbitrarily guessed value for the edge of 
the boundary layer r)^* At the downstream location x = 0,1^ 
when A = Ij the initial guesses for = 1 are 

F'Mx^O) = 0.160706 ^ G(x^O) = - 0.5^ G'Mx^O) = -1.0^ 

f'Mx,0) = 0.50238 , g(x,0) = g-'Mx^O) = 0. 

We used these initial values to integrate the non-linear 
equations of the Two-equation model across the boundary layer 
region in steps by successively increasing the range of 
integration until the correct missing wall conditions and 
the edge of the boundary layer are found by some convergence 
criterion. The solution is taken to be convergent if the 
boundary conditions and the edge of the boundary layer are 
found by some convergence criterion. The solution is taken 
to be convergent if the boundary conditions 

I (x,«>) I t I gCx,**) I , I G' (x/oo) I < 5 X 10 

and 

lf'^(x/«>)l, ig(x^~)l/ Ig'^Cx,®®)! < 5 x 10 ^ 

are all satisfied simultaneously. Once this is done at 
X = 0.1, the resulting solution provides good initial 
estimates of wall values for other streamwise locations 
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and the initial value techniques can be used to solve the 
corresponding non-linear boundary value problem. 

In solving the non-linear boundary value problem for the 
three- equation model at any streamwise location it was necessary 
to obtain very good initial estimates for the nine missing 
wall conditions which allow integration of the equations in 
the model over some interval [O/'^oo] across the boundary layer 
region. We generated these estimates by applying quasilinear- 
ization process [?] to the non-linear boundary value problem 
in this model. Before proceeding to give an outline of this 
numerical technique, we shall first write 


f = Cf, 




Ihe essential idea in applying this technique is that given 
one approximate solution to the system, we approximate 

the differential equation for f near by the linear 

differential equation in • These are 

^i+1'*’ ^i'^ ^ ^r^i+l“ t ^i+1 

-f (|71+ Q)Gi+i+ ^r^i+1” ^1+1" ^^i+1 


= cl F. + xf,.) » 


■5 i 


I G- t cl xf.)G-^,t I G-F.^,-xgpi^l+ -<^rfi+l-^i5i+l 


O ^i+1" '5 1 


I F^GC' + xGV^i * ^i • 
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(I F.+ I f.+ C? f; + I F'' 


-i+1^ '5 S "r ^"i""i+l "5 i 


5 %*^l+i'*‘^5 2xf f j|_f 


7 


- f ^i+l-^ cf ^+Q)gi^l+ ^r^i+l- ^l^/+l“ ^^i+l 


_ ^ P^^T? j.® V* ' -f 4- VP f"- — f^F^- X (f^5^ 

“5 i i^5 i i^ i 5 i i ^ ^ 


+ xFj_' fj- xF^f^ » 

I 9iVl + I 9i'^i+l- ‘9l ^ £.+ xf^) G”i+(fG» +xgi')£j^i 


- X + ^fp^i+l+^S ^i+1'*’ ^i+1 


- - n5i+i= I + f ^r 


' 4. i aV 


+ Xf^g^' + - xF< g^ 5 


s?;!- 1 n "w- *1 r ^1 nil- K «2 =1^1* r 


16 


*'i^l+i-^ r ^i^lii-" t ^ ''i'^i+i- r ^1 ^1+1 


+ 1 - r 5i+i+ 'I”* “’51+1 


5 "i"i+l 5 


f^P'- “ F.fV 
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1 # 3 •• 

— rf^ .4- — n 






■"2figi+ g^''f-.+ G^'f .-2F'5' +|p.gV (3.42) 

XID J-1 XI D XI* 


The linearized ecjuations comprise a linear two— point boundary 
value problem for along with boundary conditions 

^i+l^x^O) = Fl^i(x,0) = 0 , {x,0) = 1 


i+1 


i+1 


^i+l^^*'°^ “ f<^i(x,0) = g'_^^(x,0) = O 


'i+l 


i+1 




and 




^4 Li ^ (X,0o) 

i+1 i+1 


= 0 


* 1+1 = 9l+l'=^"> = 9l+l‘^'~> = ° 


^1+1^^'"^ = 54^.1 (x,») =5;4., (x,») =0 


'i+1 


i+1 


(3*43) 


As usual with Newton's method^ the hope is that the sequence 
of functions 

^o' ^1 ' * * * ^i^ * * * converge to solving the 

original problem. Thus, we have reduced the original non- 
linear boundary value problon to the transformed problem: solve 



121 


a sequence of linear boundary value problom for 5 #••• 

converging to From the equations (3 *42) it is clear that 
out of 18 boundary conditions required only 9 are prescribed 
at V = 0, As described in Chapter II/ we proceed by assuming 
the solution into 10 linear parts, 9 of which are complemen- 
tary solutions and one is a particular solution# Complementary 
solutions are obtained by dropping non— homogeneous terms from 
the equations and assuming one of the unknowns as unity and 
the rest zero at a time. The particular solution is obtained 
by retaining the non-homogeneous terms and taking all the 
unknown conditions as zero# Then the solution Y is given by 


■i+1 


9 

= S 
r=l 


^h. 


+ Y„ 


P 


where Y, , r = 1 to 9, are complementary solutions and Y 
r ^ 

the particular solution. The constants C^(r = 1 to 9) can 
be obtained by satisfying 9 boundary conditions at TJ^# To 
start the quasilinearization process at the downstream 
location x = 1.0 when A = l.O, we supplied the initial 
profile by combining the solution for F/G, f,g as given by 
the two-equation model at that location together with 
f = g = 0 over the interval of integration [O/l] . If Idie 
maximum of the absolute difference between the two trial 
solutions for the missing initial conditions is greater than 
e = 10“®/ the process is repeated till the desired accuracy 
is achieved. From Tablei(a)/we note that the iterative process 
converged in only 3 trials. We then form the initial profile 
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Yq for an application of the technique for the extended 
interval [0,2] by combining the stored final trial solution 
Yp over [O^l] with Y^(7l) = Yp(l) for 1 < 7? < 2 . We 
observe again from Tablel(b)that the quasilinearization 
process converged for the same e in 5 trials. With the 
final trial values for missing wall conditions thus obtained 
in quasilinearization process, we now switch over to the initial 
value technique and solve the non-linear boundary value problem 
for the three-equation model at x = 1.0 when A = 1.0. The 
initial value technique is taken to converge if 

I (x/«>) I , I G(x/«>) I / 1 O' (x/<») I <5 X 10 ’“'^ 

|f''(x,‘»)|, 1 g(x,«»)| / 1 g' (x,*’) 1 < 5 x 

and 

|f'"(x,«’)l/ I gCx^oo) I ^ I g' (x,<») 1 < 5 X 10 ^ 

are all satisfied simultaneously. The solution thus obtained 
at X = 1.0 was used to solve by initial value technique at 
other location and the process is thus extended to all other 
streamwise locations. Our nxomerical experimentation shows 
that the computer time taken at any location in the streamwise 
direction to solve the three— equation model is very large in 
comparison with that for the two-equation model. Moreover, 
as will be seen in the next section, the two-eqpation model 
gives result for temperature distribution and heat transfer 
rate very accurately while the skin friction and velocity 
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distribution are of accuracy somev/hat enough for the purposes 
of engineering applications. Hence we found solutions to the 
two-eguation model only for other cases A = -1 and A = 1 
when 0 = 0,72/ 1,0 and a = 1,0 respectively, 

6, RESULTS AND DISCUSSION 

Numerical results were obtained for Prandtl numbers 0.72 
and 1,0 when A = 1,0 and -1,0, The ntamerical values of 
missing wall conditions for the two and three— equation models 
are listed in Tables 2 ^3, 4, 5/6. The tabulated information 

has been used for the preparation of graphs which depict the 
trend for skin— friction and heat— transfer variations at the 
surface in the streamwise distance x from the leading edge 
of the plate. 

First/ we shall consider the case when A = 1 . Figures 1 
and 2 show skinfriction and heat -transfer at the plate as a 
function of x. The results from the three-equation model 
are seen to be in close agresnent with those from the 
Selected Points method of Lanczos, The two-equation model, 
on the other hand, is observed to provide results that are 
enough for practical purposes in engineering applications. 

We note here that near the leading edge, the results from the 
two and three -equation models agree well among themselves and 
with the series solution of Jones, The skin-friction trends 
as predicted by the local non-similarity method of solution 
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for o = 0*72 and 1»0 are consistent with the expected effects 
of buoyancy force after a short distance x in downstream 
of the flow# This component of buoyancy force modifies the 
shear stress at the wall at a short distance x in the 
sense that the decreasing values for skin— friction near the 
leading edge of the inclined plate starts increasing as seen 
from the figure 1 * An inspection of the same figure also 
reveals that at a fixed point x the buoyancy force is seen 
to have a larger influence on the friction factor for a = 0#72 
than for a = l.O* This is only to be expected as air is 
comparatively of lesser density and hence greatly sensitive 
to buoyancy force effect# thereby causing large changes in 
velocity gradients at the •wall which is thus reflected in the 
wall shear results# Also shown in figure 1(b) are results 
for large downstream distances and the series solution results 
at -those locations as obtained by Jones are also included for 
comparison# It is seen that -there is a good agreement between 
results from our me-thod of analysis and the series ej^ansion 
in this region also# 

Within the scale of the figure# the results frcxn the 
two and three-equation models for heat transfer coefficient 
at various streamwise locations are indistinguishable and 
therefore they are plotted as a single curve in figure 2# 

The agreement between results from local non-similarity and 
selected points method as seen frcsn the figure 2 indicates 
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the high accuracy of our method of solution# Also# at large 
downstream distances# the series result coincides with that 
given by the local non- similarity method, (see fig. 2b) 

The close agreement of results for skin-friction and heat — 
transfer between our method of solution and that used by 
Jones is somewhat remarkable and quite surprising* It is 
thus shown that the non-similar equations formulated on the 
basis of a horizontal scale#- v^en analysed by the local non- 
similarity method of solution predicts with reasonable 
accuracy the boundary layer characteristics of the flow in 
the entire region as covered here. Again# an inspection 
of the figure indicates that the heat transfer rate at the 
surface for a =1,0 is larger than that for air. This is 
because as the Prandtl number increases the thermal layer 
thickness decreases and the surface temperature gradient 
increases which results in a high rate of heat transfer 
from the surface. 

Representative velocity and temperature profiles for 
several values of x for O = 0.72 and 1 .0 are shown in 
figures (3;i2). The velocity gradient at the wall increases 
in downstream direction of the flow. This is accompanied "by 
increase in velocity near the surface and decrease in 
thickness of the velocity layer. The results of the present 
analysis are compared with those of Jones in figures 11 
and 12# the coordinates of the figure being those interpreted 
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in terms of that author. In terms of the present analysis^ 
the ordinate and abscissa groupings may be written as 

V = f F = F 

fj n 

The velocity profile calculated in the two-equation model 
at X = 4.5 is seen to agree satisfactorily with that given 
by using the selected points method of Lanczos for the non- 
similar formulation of the problem^ wherein the neglect of 
induced pressure effects has been assumed. In our coordinate 
system^ we find the agreement between the three-equation 
model and the selected points method of Lanczos for the 
velocity profile at x = 10.0 is highly satisfactory as seen 
from figure 10. Also seen in the same figure is that the 
three-equation model somewhat over corrects the t'wo— equation 
model and provide results of significantly greater accuracy. 

With regard to temperature profiles^ the temperature 
gradient at the wall is seen to increase as x increases 
with an accompanying decrease in the thermal layer thiclcness. 
It is noted from figure 9 that the temperature gradient at 
the wall for a = 1.0 is somewhat steeper than that for 
a = 0.72 as is to be expected. Since the two and three ~ 
equation models predict temperature distributions v^ich are 
almost identical, we plotted this as a single curve at 
X = 4.5 in terms of Jones variables (see figure 12). The 
comparison of results from the local non- similarity with the 
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Jones results shows very good agreement and thus the capability 
of the local non— similarity method to predict boundary layer 
characteristics is again supported. 

We shall now examine the free convection phenonena 
when the angle of inclination of the plate is negative and 
given by A s= Por this case, the buoyancy force component 

along the surface now opposes the motion and hence both the 
Nusselt number and skin-friction factor decrease in down- 
stream direction of the flow as seen from figure 13 and 14 
and ultimately the flow separates. The two-equation model 
predicts this point of separation at x = 3.65 for a = 0.72 
and X = 3.15 for a = 1.0. It is noted from the figure 13 
that the separation of the flow is not accompained by any 
singularity. But the numerical method converges near 
separation point only after a large number of iterations 
than that is required in the region near the leading edge of 
the plate. Analogous to the characteristic behaviour near 
the separation point — the skin friction and heat transfer 
coefficient were found to behave in a regular manner - was 
also encountered in other boundary layer studies such as in 
the study of interaction between a laminar boundary layer and 
a moderately strong shock [^9] and also in the study based 
upon an incompressible boundary layer calculation where it is 
insisted that the displacement thickness associated with 
the boundary layer behaves regularly at separation '[lO ] . 
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Contrary to this findings for flow under consideration, most 

boundary layer calculations including those with heat transfer 

exhibit singular behaviour at the point of zero skin— friction. 

This was the case for the mixed convection flow over an 

isothermal vertical plate placed in a uniform stream parallel 

to it with buoyancy force opposing the motion* Merkin [ll ] 

has shown that the skin -friction behaves like x -x near 

s 

the separation point x . This is the classical singular 

o 

behaviour to be expected when there is no heat transfer. 

Also Buckmaster, [l2] , has proved analytically that when 
there is heat transfer at a plate cooler than the fluid 
outside the boundary layer and the pressure gradient is 
prescribed, then a singularity of the type '^x.-x log (x -x) 
may be expected* In the light of these discussions, it is 
no wonder that the local non-similarity method as applied 
to this problem can go up to separation point unlike in the 
case of forced flow over a circular cylinder [sj • While 
the selected point method of Lanczos predicts the point of 
separation of the flow for air at x = 3*704, our method of 
solution puts it at x = 3*65 and thus deviating from exact 
result just by per cent only. 

The velocity and temperature profiles in this case, 

A = -1, are shown in figures 15-22 for several values of x 
and a = 0.72 and 1*0. The effects of opposing buoyancy 
force on the fluid motion is to reduce the velocities and 
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the velocity gradients at the wall in the downstream 
direction* The velocity profiles corresponding to the 
locations x = 3.65 and x = 3*15 for a = 0.72 and 1.0 
respectively are typical of that which is obtained at the 
point of separation. With regard to temperature profiles, 
the temperature gradient at the wall now decreases with 
accompanying increase in the thermal boundary layer thickness 
for increasing downstream distances* 
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Table 1 


Convergence history of guasilinearization 
process applied to the three— equation 
model at x = 1.0 and a = 0.72 

1(a) Interval of integration [o,l] : 


Missing 
wall values 
at rj = 0 

Trial Number 

1 

2 

3 

F' ■'( x,0) 

0.582515 

6.582444 

0.582443 

G(x,0) 

-0..506109 

-0.506795 

-0.506796 

G'^(x,0) 

-1,011215 

-1.009053 

-1.009052 

f"( x,0) 

0.147377 

0.148055 

0.148055 

g(x^o) 

0.001001 

0.000674 

! 0.000674 

g* *{x,0) 

-0.003494 

-0.002439 

-0.002439 

lE"(x/0) 

-0.030106 

-0.030003 

-0.030003 

g(x,o) 

-0.000183 

-0,000305 

-0.000305 

g"(o) 

' 0.000659 

0.001095 

0.001095 
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Kb) Interval of integration [o,2] 


Missing 
wall values 
at T) = 0 


Trial Number 


F"(x,0) 

1.133427 

G(x,0) 

-0.904286 

Q 

O 

-0.596535 

f"( x,0) 

0.220240 

9 (x/O) 

0.019520 

g"(x,o) 

-0.018744 


-0.038858 

g (x,o) 

-0.009545 

g'^'C x,o) 

0.008501 


t 


1K72933 

•0*948659 

-0.557730 

0,232218 

0.011932 

-0.012028 

•0.044931 

•0.005668 

0.005210 


3 

4 

5 

1.173112 

1U73112 

1*173112 

1 

-0.948649 

-0.948649 

-0.948648: 

-0.557663 

-0,557664 

-0.557664 

0.231951 

0,231953 

0.231954 

0.011946 

0.011945 

0.011945 

-0.012005 

-0.012005 

-0.012005 

-0.04474$ 

-0.04474-? 

-0.044746 

-0.005615 

-0,005615 

— 0.005 ol 6 

0*005161 

0,005161 

0.005162 
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Table 2 

Missing wall values for the two-egaation model at 
various streamwise locations when A = 1.0 and O = 0.72 


X 

F"(XyO) 

G(x/0) 

G' (x.O) 

f"{x,0) 

g( x,0) 

g' (x,0) 

0.01 

1.012956 

-1 .719960 

-0.360899 

1.724501 

0 . 75305.1 

-0.186420 

0.05 

1.069364 

-1.696415 

-0.366502 

0.907459 

0.378576 

-0.095723 

0.10 

1.116645 

-1.677632 

' 

-0.371126 

0.688496 

0.276625 

-0.071170 

0.35 

1.272207 

-1.620259 

-0.385624 

0.454658 

0.166785 

-0.044162 

0.50 

1.342495 

-1.596513 

-0.391889 

0.361759 

0.123178 

-0.034248 

0.95 

1.513424 

-1.543596 

-0.406458 

0.278862 

0.084736 

-0.024862 

1.0 

1.529997 

-1.5 38736 

-0.407811 

0.273060 

0.082115 

-0.024214 

1.5 

1.680533 

-1.497327 

-0.419839 

0,230984 

0.063333 

-0.019518 

2.0 

1.810740 

-1.464564 : 

-0.429767 

0.204758 

0.052066 

-0.016637 

2.5 

1.927391 

-1.437320 

-0.438322 

0.186252 

0.044410 ■ 

-0.014637 

3.5 

2.133092 

-1.393490 

-0.452706 

0.161074 

F 

0.034513 

-0.011979 

4.5 

1 2.313399 

-1.358878 

-0.464654 

0.144211 

0,028308 

-0.010253 

i 

5.0 

2.396558 

-1.343969 

-0.469977 

1 0.137598 

0.025985 

-0.009591 

6.0 

2.551901 

1-1.317677 

-0.479260 

0.126757 

0.022328 

-0.008526 

8.0 

2.829239 

-1.275195 

-0.495994 

0.111112 

0.017421 

-0.007043 

10.0 

3.074487 

i -1.241 664 

-0.509654 

0.100127 

0.014266 

-0.006045 

12.0 

3.296487 

1-1.214052 

-0.521440 

0.091848 

0.012064 

-0.005321 

14.0 

3.500643 

-1.190646 

-0.531842 

0.085311 

0.010438 

-0.004768 

16.0 

3.690517 

i-1 .170367 

-0.541174 

0.079976 

0.009188 

j 

-0.004329 

18.0 

3.868646 

[ 

i-1 .152512 

-0.549652 

0.075512 

0.008198 i 

-0.003972 

20.0 

f-, 03689 2 

-1.136586 

-0.557431 

0.071704 

0.007394 

-0.003675 

25.0 

4.422832 

Ll. 103071 

-0.574515 

[ 

0.064192 

0.005922 

-0.003111 
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Table 3 


Missing wall values for the three- equation model at 
various streamwise locations when A =1.0 anda=0»72 


3(a) 


Missing 

wall 

values 

X 

0,01 

0.05 

0.35 

0.95 

1 .0 

F''(x,0) 

1.01198 6 

1.065471 

1.263427 

1.497710 

1.5139758 

G(x,0) 

-1,7202 37 

-1.715802 

-1.622007 

-1,545579 

-1.540737 

G''(x,0) 

-0,360887 

-0.362808 

-0.385 377 

-0.406105 

-0.407474 


1.869169 

0.978137 

0.454658 

0.304145 

0.297836 

1 

g(x,o) 

0.829650 

0.41979G 

0,166785 

0.095255 

0.092318 

g' ' ix,o) 

-0.201031 

-0.101215 

-0.044162 

-0.027186 

-0.026478 

' (x/O) 

-54.410390 

-5.647714 

-0.371913 

-0.093376 

-0.087040 

g(x,o) 

-23.711487 

-2.558342 

-0.178260 

-0.0449 36 

-0.041819 

g^'C x,0) 

6.154328 

0.635890 

0.044219 

0,011178 

1 

0.010414 






- 
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3(b) 


iM is s ing 

wall ] 

values 2*0 3.0 

F'^(x,0) 1,782448 2.006773 

G(x,0) -1.483585 -1.414805 

G'-'Cx^O) -0.424620 -0.445793 

f'"(x,0) 0.221720 0.187483 

g(x,0) 0.058825 0.043658 

g"(x/0) -0.017794 -0.014386 

l"(x/0) -0.033752 -0,019734 

g(x,0) -0.015539 -0.003841 

g"(x,0) 0.003876 0.001073 


X 

4.5 5.0 10.0 

2.281138 2.515875 3.031266 

-1.35881 6 -1.316811 -1.2 39290 

-0.464810 -0.480017 -0.510572 

0.1565577 0.137270 0.107781 

0.031723 0.024914 0,015746 

-0.011194 -0.009285 -0.006540 

-O.OllSOS' -0,007849 -0.003966 

-0.004538 -0.002862 -0.001217 

0.001252 0.000820 0.000379 






3(c) 
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Missing 

wall 

values 


F"' (x/0) 

G( X/O) 

G' " (x,0) 

f'Mx,0) 

g(x^o) 

g^'(x,o) 

?"(x,0) 

g(x,o) 

g' ' (x,0) 


12.0 


3.250478 

•1.211157 

•0.522571 

-0.098623 

0.013252 

-0.005741 

-0.003101 

0.000886 

0.000286 


14.0 


3.452147 

-1.187328 

-0,533159 

0.091401 

0.011417 

-0.005131 

■0.002516 

•0.000675 

0.000224 


16.0 

o 

# 

CD 

20.0 

3.639764 

3.815818 

3.982125 

-1.166708 

-1,148574 

-1.132414 ; 

-0.542656 

-0.551283 

-0.559196 

: 0.085517 

0.080601 

0.076415 

0,010012 

0.008902 

0,008004 

-0.004648 

-0.004256 

-0,003931 

-0.002097 

-0.001784 

-0.001542 

-0.000531 

-0.000429 

-0.000353 

0.000181 

0.000150 

0.000127 
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Table 4 

Missing wall values for the two-equation model at. 
various streamwise locations when A = 1,0 and o = 1.0 


X 


G(x,0) 

G'(x,0) 


mm 

g'(x,o) 

0.01 

0.897840 

-1,556704 

-0,394791 

1.662952 

0.762325 

-0.225220 

0.05 

0.952344 

-1.533030 

-0.401577 

0.875860 

0.381551 

-0.115461 ■ 

o 

• 

o 

0.997883 

-1,514165 

-0,407106 

0,664845 

0.277932 

-0.085758 

0.35 

1.148327 

-1.457100 

-0.424540 

0.403675 

OJ. 48450 

-0.048685 

0.75 

1.313259 

-1.402417 

-0,442347 

0.297050 

0.096073 

-0.033511 

0.9 

1.365020 

-1.386699 

-0.447686 

0.275809 

0.085873 

-0,030508 

1.0 

1.397551 

-1,3771 38 

-0.450985 

0.264179 

0,080349 

-0,028870 

2.0 

1.668790 

-1.305718 

-0.476958 

0.198017 

0,050219 

-0.019710 

in 

• 

1.781351 

-1.279797 

-0.487019 

0,180057 

0.042612 

-0.017299 

3.0 

1.884238 

-1.257697 

-0.495893 

0.166433 

0.037076 

-0.£)15504 

3.5 

1.979606 

! -1.238433 

-0.503863 

0.155602 

0.032845 

-0.014103 

4.5 

^ 2.153136 

-1,206048 

' -0.517789 

0.139214 

0.026771 

-0.012034 

5.0 

2.33085 

-1.192163 

-0,523975 

0.132786 

i 

0.024509 

-0.011242 

6.0 

2.382333 

-1.167790 

-0,535165 

I 

0.122248 

1 

0.020965 

-0.009972 

7.5 

2.585286 

-1.137501 

-0.549700 

0. 110309 

0.017215 

-0.008579 

10.0 

2.883433 

-1.097999 

-0.567800 

0,096376 

0.01 3232 

-0.007026 

12.0 

3.095916 

-1.072863 

-0.583328 

0.088341 

0.011144 

-0.006171 

15.0 

3.383422 

-1.042151 

-0,600712 

0.079294 

0.008986 

-0.005249 

17.5 

3.601155 

-1.021040 

-0.613251 

0.073523 

0.007722 

-0.004684 

20.0 

3.803352 

-1.002860 

-0.624460 

0.068824 

0.006760 

-0.004240 


1 

1 



— 

mm .j-iirti-.r-J 
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Table 5 


Missing wall values for the "two-equation model at 
various streamwise locations when A = -l.o and a = 0.72 


X 

P"'(x,0) 

G{x,0) 

G'Cx^O) 

f''(x,0) 

g(x,0) 

g'(x,0) 

0.01 

0.943830 

-1.750323 

-0.353826 

-1.719238 

-0.795719 

0.191783 

0.05 

0.887897 

-1.776237 

-0.347923 

-0.900073 

-0.437562 

0.103122 

0.1 

0.841584 

-1.798698 

-0.342903 

-0.679722 

-0.344704 

0.07969 3 

0.2 

0.772022 

-1.834286 

-0.335747 

-0.511752 

-0.277228 

0.062213 

0.55 

0.605718 

-1.929735 

-0.315390 

-0.333371 

-0.214694 

0,044624 

1.0 

0.456267 

-2.030749 

-0.295935 

-0.253277 

-0,195101 

0.037544 

2.0 

0.228145 

-2.222106 

-0.262401 

-0.173744 

-0.187952 

0.031515 

2.5 

0.144948 

-2.305768 

-0.248838 

-0.149781 

-0.189038 

0.029916 

3.0 

0.075563 

-2.381886 

-0.236951 

-0.130560 

-0,190945 

0.028689 

3.5 

! 0.17220 

-2.45035 

-0.226594 

-0.114546 

-0,193116 

0.027695 

3.6 

0.00668 

-2.46309 

-0.22469 

[ 

-0.11166 

-0.19355 

0.02752 

3.65 

0.00125 

-2.46923 

-0.22370 

-0.11034 

-0.19377 

0.02745 



Table 6 

Missing wall values for the 
various streamwise locations 

X P"(x,0) G(x,0) G'(x,0) 

0.01 0.831200 -1.587410 -0.386230 

0.1 0.732990 -1,636616 -0.373036 

0.2 0.666422 —1.672996 —0.363666 

0.55 0,508456 -1.771033 -0.339859 

1.0 0,368555 -1.874817 -0.316634 

1.5 0,251956 -1.976660 -0.295481 

2.0 0,160522 -2.068045 -0.277639 

2.5 0.086696 -2.149706 -0.262430 

2*75 0.054921 -2.186870 -0.255688 

3.0 0.026692 -2.219176 -0.250236 

3.15 0.009902 -2.241192 -0.245980 


;wo- equation model at 
when A = -1 iO and <^ = 1.0 


f"(x,0) 

-1.655614 

-0.652926 

-0.490496 

-0,317114 

-0.238451 

-0,191135 

-0.159536 

-0.135806 

-0,125898 

-0.117154 

-0,112070 


g(x^o> 

-0.8099 31 
-0.353640 
-0..285961 
-0.224295 
-0.205926 
-0.200758 
-0.200465 
-0.201903 
- 0.202881 
-0,203735 
-0.204562 


g'(x^o) 


0,232119* 

f 

0.096678 

0.075578 

0.054337 

0.045732 

0.04105>9 

0.038241 

0.036207 

0.035386 

0.034678 

0,034267 
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Fig. 9 Temperature profiles when A =10. 



Fig. 10 


Velocity profiles when A=1-0 

x = 10-0. 
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Pig. 11 Velocity profiles when A =1* 0, 0“= 0*72 
and X = 4*5. 
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Fig. 12 Temperature profiles 
and x = 4-b. 
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Fig. 19 Velocity profiles when A =-1-0 and 




Fig.20 Temperature profiles when A=-1-0 and <r=10. 








CHAPTER IV 


POSSIBLE S3MILARIIY SOLUTIONS FOR 
LAMINAR FREE CONVECTION ON 
HORIZONTAL PLATES 
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1. INTRODUCTION 


Natural convection flow adjacent to horizontal surfaces 


bounded by an extensive body of fluid above is of considerable 


importance in micrometerological and industrial applications. 
Most of the earlier theoretical studies [l/2/3, 4/5/6] are 
devoted to the problem of finding similarity solutions for 
steady free convection over horizontal surfaces j in all of 
these analysis/ boundary layer approximations are enployed 
and power form variations of wall tanperature distribution 
are assumed. Recently/ a few analytical studies dealing with 
unsteady free convection flow over semi- infinite horizontal 
plate have appeared [7/8/9] • But they all deal with the 
cases when the plate temperature oscillates in time about a 
constant non-zero mean tanperature. In view of what has been 
done so far, it leaves us Hi doubt as to whether these studies 
have eidiausted all possible similarity solutions or not for 
free conveatlon over semi-infinite horizontal plate. 


in this chapter, we derive in a systematic way all 

*1 ■fnr' "hoth st03.c3,y Slid uiist03idy 

possible similarity solutions for botn suea y 

T flow ovcr a semi— infinite 

free convection boundary y 

a. 1 1 +-ra The solUtlons presented here were obtained 

horizontal plate, soru 

/ o^-ific wall tanperature distributions/ but 

not by assuming specific wa f 

Finn forms freon similarity 

rather by deriving the solu 

This approach not only results in more general 


conditions , 
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foirns thsn hithsirto obtcimsd,^ but iTssults in 311 othsn 
possibl© solution forms » These solutions ate of impojrtance 
as they represent test cases for approximate and numerical 
approaches. Also/ they could as well be used for under- 
standing the effects of pressure gradient and heat transfer. 

The method employed here is similar to that used by 
Hansen [lo] / Yang [ll] and Nanda [ 12 j| for certain other 
laminar boundary layer problems. It consists of defining a 
new independent variable and three dependent similarity 
variables in terms of the generalized transformation variables. 
When these similarity variables are substituted into the 
boundary layer equations, they are reduced to a set of 
ordinary differential equations only if the coefficients of 
the similarity variables and their derivatives can be made 
constant. These coefficients when made constant are referred 
to as similarity conditions; they usually contain derivatives 
of the transformation variables and are thus differential 
equations themselves. Any solution to the similarity 
conditions which will result in explicit forms for the 
transformation variables, will transform the original set of 
partial differential equations into a set of ordinary diff- 
erential equations and these are called similarity equations. 

We conclude from our investigation that there are in 
,all .fo ur possible ca ses of similarity for free convection 
boundary layer flow over a heated horizontal plate. While 
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this analysis has verified that the similarity possibilities 
for steady cases have essentially been covered in the literature/^ 
it also opens a new possibility that somewhat more general 
wall thermal conditions could also be taken into account. The 
other two new possible cases deal with unsteady conditions. 


2 BASIC EQUATIONS 


Let us consider two-dimensional laminar free convection 
boundary layer flow over a semi- infinite horizontal flat plate# 
placed in an infinite expanse of Boussinesq fluid. The 
governing equations to the flow are then given by 


^ = 0 , (4.1) 

3y 

+ V ^ (4,2) 

ay 


(4.3) 


(4.4) 


Where f and \ are the temperatures of the boundary layer 
and free stream respectively, 5 and v are the velocity 
components along the 5 and y directions respectively, x 
being the coordinate along the plate and y normal to It, 
g is the acceleration due to gravity, 5 Is the pressure 


a|+uai + va|s-i ai 
ot 3x ay ^ 3^ 


i iii * g0 (T - Tj# 

ay 

af - af . - if - oc — 
at 3x ay 9y 
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inside the boundary layer, t is time, 0 is the coefficient 
of thermal expansion, a is the thermal diffusivlty, p 
1s the density and V is the kinematic viscosity. In our 
discussion the heated plate with tanperature T^(xit) is 
taken to be facing upwards and is higher than 

While the initial conditions generally depend on the 
particular problem, the boundary conditions are given as 

y=0» u=0=v,T = T„(x,t) or ~ t 

ay 

yssoo* u=0, p = 0 and T = • (4.5) 


Here q is the prescribed surface heat flux and k the 
thermal conductivity. 


We now introduce the following transformation 


t a , X « , y 

L 


XXL „ vL , 


g0L^(T - 5J 


G a 




-t 2 

n - 

/ P - 2 

pv 


(4.6) 


to render the equations (4.1) to (4.4) in dimensionless form: 


a2i 

ax 



0 , 


a t 


4- U 


au 

a X 


+ V 


ay 


- - 

“ ax 


9Ji 

3Y^ ' 


s G , 

ay 


ac 

at 


d G 


+ V 


aG _ 1 if® 
3y “ ® ay2 


(4.7) 

(4.8) 

(4.9) 


( 4 . 10 ) 
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The boundary conditions now are 

Y - Oi n = 0 = V, G = G (x,t) or ^ * -q(x^t) j 

W oj 

y>eoo| 11=0^ p3=0 cind G 5= 0 # 

We introduce the stream function / defined as 
d'i ^ 

u == II and V = - / 

in equations (4*7) to (4.11) to get 

*yt y xy x yy x yyy 

Py = G / 

°t ^y ®x “ \ ^y “ a '^yy ' 

with 

/ N 9G 

y.O* «*«y = 0t ® = 3y “ 

y»<»Jfy*p = G*0, 


3i SIMILARIIY CONDITIONS 

i_ ^ +-he similarity analysis is to 

The first step in the simu-a 

a similarity variable V In the 


n 



(4*11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

C{,(X/t) / 

(4.16) 

introduce 

(4.17) 


and the stream 
G are talcen as 


function * / pressure p 


and the temperature 
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f(7l) 

■ f 

(4.18) 

7T(T|) 

= ' 

(4.19) 

eiv) 

G 

" G (x/t) * 

w 

(4.20) 


The proble:n now essentially reduces to determining the functions 
9^(x,t), and G^(x.t) so that the partial 

differential equations (4.13) to (4.15) are reduced to 
ordinary differential equations for f(V)/ and 

with proper boundary conditions to be obtained from (4.16). 
Though we will be deriving possible forms of the similar 
solutions for the case of prescribed surface tenperature 
variation, the corresponding cases for prescribed surface 
heat flux distribution will also be indicated. 


We now pass on 


to the key step in the analysis by 


substituting the identities 


(4.21) 


V 



O Cp ^ a 1 

- [(-33;)f +<Pi ' ax' 


P = V ' 


(4.22) 

(4.23) 




and 


(4.24) 
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in eqruations (4*13) to (4*15) so as to reduce them into the 
form 


.2 

f"'-a2fE"'-(a2+a3)f"+ a^ff * '-(a^+a^) f ' -a^^Jr-agTin'* 0 > (4..25) 

n' = a^e k (4 *.26) 

e^'-CagTI-a^f) 0'- (ag + a^f^)e = 0 (4*27) 


where 


a< = 


3 % 
C-y^) 


"1 “ ' 
1^2 


d(p. 


1 

1- ^^2 
®3 = ^57- '"rt’ 
1 2 

^4 * ^ ' 


8<P, 


®5 * ^ 


(-^) f 


i- i- (!2w) , 

e To ?r at'^ ' 


<l>2 =w 


* 


<(>2 3 G , 


(__S!) 

«’l°« 


^3-._ (1!^) , 
a« = -T~ ^ 3 X 
8 


w 

Sg * (p^<p^ 


(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 

(4.33) 

(4.34) 

(4.35) 


(4.36) 
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The boundary conditions with prescribed temperature variation 
are 

f(o) s 0 = fMo), e(o) = 1 1 

f'(oo) = 7 r(oo) = e(oo) = 0. (4,37) 


It is clear now that the similarity solutions are possible 
if and only if a^ to ag are all constants. Thus the 

ecjuations (4*28) to (4*36) comprise of what is known as 
similarity conditions for our problem. Since the number 
of these conditions exceed the niimber of unknown functions^ 

determined/ not all of these constants 

can be taken arbitrarily. For instance, it can be verified 
that 


a^ ^^2 ^ 

By « 2ag + 3a^ » 

a^Sg** ^2a^+ ^5)^3 * 

a^agW 2aga^+ • 


(4.38) 

(4.39) 
(4*40) 

(4.41) 

(4.42) 


A general method of determining the unknown functions 

differential equations (4.29) to (4.32). Then equation 

(4.35) yields (Py The surface temperature variation for 

is oossible can be obtained frxxn 
which a similarity solution is posi.x 


(4.36). 
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p,,,-.TCT,E FOmS_ OT:. -TUTT.APTTf SOWIIOMS 

, Steady .ree convection wit. surface 

Case 1 • 'jT-t-au.jr ^ 

Power o£ a Linear Function of x. 

Varying with any 


Since we deal with steady oases. 


a =S a^r 0 • 

^3 6 


^ , .n (4 31 ) and C 4 . 32 ) 

4 « a aeneral solution to C4.3i; 

Vje derive a genercix 

combining them to yield 


d<P, 

>4 "l ilil . 
r = T "dV 

dx 


, -1 4 - i nn we have 

!§. and solving the equation 

Taking a - sr “ 


(4*43) 


SC a, 


^10 


'Wr\ XV/ — 

on substituting 

«+-ant of integration. 

a, „ being a constan ,oi ving the resulting 

Ill we obtain on solving 
(4.43) into (4.31)# 


equation for = / 1' 


<Pj = C ®11 * ®10 


4. (e-l)r-' 


(4.44) 


'f- ^ - X-v XV/ 

1. of integration 

in (4.43)/ we can write 
using this ib 




Ho ^ 


Ho 


(e-l)] 




(4.45) 


/a we obtain 

( 4 .. 35 ) and ( 4.3 ' 


From 
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G 


^ 8 ^ 9^10 


w 


[a 


11 


+ X 


"10 


(e-l)3 


(4,46) 


Here constants ag,ag, ^io'^4'^ll ® arbitrary. 

But still it is possible to assign specific values for some 
of these constants without actually restricting in any way 

the fortn of G in (4,46). Since characteristic length L 

w 

is included in x and a^^^ is entirely arbitrary in the 
sense it does not occur in any of the similarity conditions, 

a., a„ and a^ can all be assigned specific values leaving 

4 8 9 

e, and a^^j^ as arbitrary constants. 

We now introduce a constant n such that 


n = 


3+26 

•nr* 


Prom this we get 


e s 


n+3 

n-2 


(4,47) 


Let us now assume 


B.A ** 

“4 


n-2 _ 

, ag 


2^ and a^ 


= 1< 


Thus we have 


a SS 


from (4#47}# 
5 


-5 n frcxn (4»39)# 


= 


^n±2 (4 •40) 

5 


and 
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Now the forms of the functions and G are 

n~2 

cfi = + f” ^ ^ ^ 


and 


10 

n4-3 

^2 ^10^^11 ■' 


fa I X ^n 

a, ^ a,, 1 + -r- J • 


w " “10 '“11 ^ a. 


10 


Eejuations (4.25) to (4.27) can now be written as 

ft*' 4. ff'^- f* - |(2n+l)Tr + ^^rin* = 


n' = 0 / 


i. e'' + fe' - rif'e = o 

a 5 

with boundary conditions (4.37). 

The velocity components u^v and the temperature G 
the boundary layer may now be written as 


2n+l 

^ 2L.) 5 f/ , 

10 ^“11 ^ ^10 


uasa</.v(aii + a ^ 


n-2 


11 a^Q 


G = a|o(aii + * e 

and tha rate of heat transfer at the surface is 


(4.48) 

0 , 

(4.49) 

within 

( 4 . 50 ) 

( 4 . 51 ) 

( 4 . 52 ) 
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6n~2 

q= 0'(o) 5 ^ (4.53) 

10 

For n = O and a^^ =0, equations (4*49) together with 

boundary conditions (4.37) represent the steady free 

convection flow over an isothermal semi— infinite horizontal 

plate. Rotem and Claassen have dealt with this case 

and gave numerical solutions for several Prandtl numbers 

including large and small values of it. By putting ® 

and n = 1/3 in (4.52) and (4.53), one finds that the 

surface temperature varies with one third power of the 

distance away from the leading edge while the resulting 

heat-flux at the surface is uniform. Nimerical solutions 

for this case have been given by Pera and Gebhart [6] for 

1 

a range of values of Prandtl number Yq 100. Previous 
similarity solutions given by Rotem and Claassen can be 
obtained by taking a^j^ = 0 and n arbitrary. The present 
analysis shows that the numerical results discussed previously 
can be directly applied to cases where a^^^ / O and ^2.1^ 

The physical Implication of ^ ^ that the surface 

temperature at the leading edge of the plate is different 
from that of the ambient fluid, a feature which is certainly 
more realistic, since it is difficult physically to heat the 
plate without raising the temperature of its leading edge. 
However, it is to be noted that for such a case of non-zero 
a^^j^ the velocity component u 


is not identically zero along 
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X = 0, but is only zero at y = 0 and y = oo and thus we 
observe some plate length and temperature difference must 
precede x = 0 in order for the boundary layer thickness 
to build up. Thus the effect of a^^ is merely to allow 
a shift m the x coordinate. If we assume that x = 0 
corresponds to the leading edge of the plate/ we have 
a^i =0 and a simplified general form for G is 

given by 

n~2 n+3 

where = 5D^ and D^, are constants. 

For e / the boundary layer thickness is constant 
or increasing with x if n < 2 and that the same is true 
for n if n > -1/2. When n = -1/2/ the thermal energy 
that is convected at any location is constant and hence the 
surface is adiabatic and all the convected energy is released 
by a line source at the leading edge which means, physically, 
that there is no boundary layer flow. Thus the restriction 
on n, namely - | < n < 2, becomes necessary to retain the 
boundary layer approximations. 

In view of eguation (4.53), it can be easily seen that 
the present similarity solution is equally possible when the 
surface heat flux is prescribed. The condition is that the 
surface heat flux q must be of the form 
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q = (a 


Ho’ 


(4.54) 


where 0 is another constant parameter. If we let 


6n-2 


or n 


5 3+2 


■g— / these two problems then become 


identical, 


Case 2. Steady Free convection with Surface Temperature 
Varying with an Exponential Function of x. 

We deal again with steady laminar free convection 
flow, hence 


^2 = ^3 = ^6 = 


Now taking e = 1, equation (4,43) becones 


^2 = • 


The equation (4.31) yields: 


^1 “ ^12 


a, o e 


=^10 


which on substituting in (4.55) gives 

a. 


"10 


^2 * ^ 10^12 ® 


Prom (4.35) and (4.36), we obtain 

(p3 cpj 


w 


^ 8^9 (p3 (p2 
1 2 


4 

5 ^ X 
5 2 10 
4:^ ^i2^0 ^ 


(4.55) 


(4.56) 


(4.57) 


(4.58) 
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Since the integration constants are entirely arbitrary 
and a characteristic length is included in x# we assign 
=1, a^ =1 and Ug = We then have a^ = 5 from 
(4.39) and ^ from (4.40). Finally equations 

(4.25) to (4.27) are 

2 4 1 

f"'+ ff"- 2f' - i Tt - ih jr' =0, 

71' = e , 

I e"+ fe'- 5f'0 = 0 , (4.59) 


with the boundary conditions given by (4.37). 
5a , 


Putting m 
we have 


^4 5 


and using (4.56)/ (4.57) and (4.58) 



^10 ^10 



mx 


■“l 

= a ^2 e 5 , 

(4.60) 


mx 


•^2 

- 5 a e 5 

” S ^12 ® ' 

(4.61) 

°w 

5 _5 -mx 

=* ”2^12 ® • 

Wl 

(4.62) 


Also from (4,21)/ (4.22) and (4.24)/ we obtain 

2mx 

m “12 


5 2^5 

u = r a^o e , 


V ss^9t, 


12 


mx 

e ^ (f+neO / 


(4.63) 

(4.64) 


c 5 mx Q 

S = ^ ^2 ® ® 


m 


(4.65) 
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and the rate of heat transfer at the surface is given by 

6mx 

5 6 

q ^ a^ e e'(0). (4,66) 

m 

By assuming the distribution of wall temperature to be ejq^onen- 
tial# Jaluria [l3] has arrived at the same eqfuations as (4«59)» 
We note here that the parameter m is absent in the resulting 
similarity equations and this will be so, since the value 
of the parameter s is taken unity at the outset. According 
to equation (4.66), the present similarity solution is again 
seen to be identical to that of prescribed plate-surface 
heat flux 


q 


^ ,6 
5^2 ^2 


/ 


where 



Case 3. Unsteady Free Convection with Surface Temperature 
Vairying Inversely as Square Root of a Linear Combination 
of X and t: 

To our knowledge, no similarity solution for free 
convection boundary layer equations over horizontal plate 
under unsteady surface thermal conditions exists in the 
literature. But recently a few analytical studies on free 
convection over horizontal surfaces with fluctuating surface 
taciperature about a non-zero mean have been published [7,8,9] 
The similarity conditions already derived can now be used to 
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explore the forms of the possible similarity solutions under 
unsteady conditions* 

An integral of (4.29) is of the form 

= (A(x) - 2a2t)“^/^ (4.67) 

where A(x) is to be determined. Also from (4.29) and (4.30), 
we get 

(f >2 - B^x) (A(x) - 2a2t)“’^'^^ ^4.68) 

where co = ~ , a. / 0, 

^2 ^ 

Again, frcp’nxequations (4.31) and (4.32), we have 

<p2 = C(t) <P® (4.69) 

^5 

where e == — ~ , a. ;)t 0. 

Solving (4,31) using (4,69), we obtain 

a (e-1) ^ 

9^ g X + D(t)] (4.70) 


and hence from (4.69) we get ^ 

a, (6-1) 

^2 = C(t) [ 3C + D(t)] . (- 

On comparing the equations (4.67) to (4. 70), a suitable 
solution is obtainable if the following identities hold 



157 


e *= 0 ) = -1, B(x) = C(t) « a^ 3 , 


2a- 


A(x) = - T—^ X and D(t) = - 2a- 1; 

^13 2 


(4.72) 


From (4.35) and (4,36), is obtained as 

w 


G _ ^8^9 2 

°w-— ^3 


(-!!i X . 2 a,t)V2 

^13 2 


. (4.73) 


Though arbitrary values may be assigned to ag, a^, a^, 
we have to keep as parameter in the reduced similarity 
equations for this case. We now have 


3 2 5 4 

From (4.38) and (4.39) it follows that 

Eg = a2, a.^ = a^. 

Prom (4,40), we get = 0, Taking a^ = -1, Eg = -1 and 
ag = 1, the similarity equations for the flow may now be 
written as 

■£**t - + ff" =0, 

7T' = 0, 

i 0'' - (a2^-f) e'-(a2-f') 8= 0 (4,74. 

and the boundary conditions (4.37) are retained as such. 
Final forms of the unknown functions 
similarity transformation are: 
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(p 




G = 
w 


(|25 - 2 a , 

(4.75) 

a. - - 2 a t)^/2 ^ 

13 a^3 2 

(4.76) 

^13 (|x _ 2 ,1/2 • 

(4.77) 

ai3 2 



The boundary layer characteristics can now^ be written ast 


u “ ^1 3 ^ ^ 

V .(225 - 2 a ^ t )“^/2 
^13 2 


(i2 - 2a^t) 


172 


0 


(4,7S) 




The rate o£ heat transfer at the surface is 

•eMo) . 


q a - 




( 4 . 79 ) 


(22 . 2a, t) 


‘13 


It la evident that the similarity solution breaks down 
if the quantity (|f - 2a, t) is non-positive. Thus it is 

1 -3 

t- a rarefullv for numerical calculations 

important to choose ^2 

and also the ranqe of validity of the solution for x and t 

ibod beforehand. This restriction is necessary 
is to be prescribed betorenaxii. 1 . 

r: 4 -omThfsrature and surface heat flux within 

to Tceep the surface temperatur 

rv is to be discarded from 

bounds. Obviously, the leading 
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consideration in respect of this similarty solution. 
Otherwise/ at any finite x 0/ the surface tonperature 
decreases monotonically with time for a„ < 0 and 
increases with time for a„ > 0 and t < ) » As a 

and a ^2 arbitrary constants/ this time limit is not 

very restrictive. Since/ physically/ the parameter a^ 
is an indication of the unsteadiness of surface tanperature 
variation/ this similarity solution could prove useful in 
determining the validity of quasi-steady solutions. By 
taking prescribed surface heat flux as 

9 = ^3 < 1 ^ - 2 * 2 ^’ ' 


we arrive at the same similarity solution as given above. 

Case 4, Unsteady Free Convection with Surface Temperature 
Varying Directly as Square of a Linear Function of x and 
Inversely as ^ Power of a Linear Function of t» 

We shall consider here a^ = 0. Now from equation 
(4.29)/ we have 

where a^^^ is a constant of integration. 

From (4.30)/ we obtain 

= E(x) - 232 1) / / 


where E(x) is to be 


determined from (4.32) using (4.80) 
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The equation for E(x) is 


2a2t)-V2 ^ g 


With 6 )= 1, E(x) = a^x + ^2 becomes 


<P 


a^x + a^5 

2 u”4.2a2t)^/2 • 


(4.81) 


From (4.40)/ we have 33 .= 9 since a^ /£ 0 * If a^= 0 , then 
which in turn gives G as functions of time alone 

and thus a^ = 0 and a^ = 0. But then, it implies that 
the flow is not driven by the "indirect'"' mechanism. Thus 
a^ = O will lead tO/ physically# unacceptable solution for 
free convection flow over a semi- infinite horizontal plate. 
Hence / 0. Now we obtain by combining (4.36), 

(4.28) and (4.34), as 


^1^ 


(agX + aj^g ) 


G « 

« (3^4 - 2a^t) 


572 • 


(4.82) 


In this case, we have a^ = 3.^ and a^ = 5 a 2 * The parameters 
a 2 .a 5 ,a^,ag all have to remain in the similarity equations 
which are now 


f'" - 2a^f'+ ag 


ff" — a-f' — a-^ = 0 / 


Ti' = ag9. 


i 0"- (a2^-*a5f5®^ * ^^^2"^ 2a5f')e = 0. 


(4 ,83) 
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The boundary conditions are the same as given by (4»37), 
Also we have 


u 


^ 5 ^ ^15 

^14 “■ 




V = - 


a f 




Tn ' 


^ 1^9 

2a^ 


(a^x + a^^g) 


(ai4-2a2t) 




0 


(4,84) 


and the rate of heat rransfer at the plate is given by 


^ 1^9 ^^ 5 ^ ^ 1 5 ^ 

'iS'a'r 


(ai4 - 2a2t)^ 


0 ^0). 


(4.85) 


Similar to previous case/ constants a^^ and a2 should be 
so chosen that the quantity ^^14'" ^^2^^ remains always 
positive. It is to be noted that for a2 < 0/ the surface 
temperature ( 4 . 82 ) at any x changes directly with time 
while this trend is reversed for a2 > 0«- For this latter 
case, a limit on time t is not necessary since the arbitrary 
constant can be assigned large values. Similar to the 

previous case/ this similarity solution is also applicable 
to the problem with surface heat flux prescribed by 

’ “ ^ (a,, - • 
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In view of the generality of the present analysis, it 
may be concluded that the four cases discussed here represent 
possible similarity solutions for laminar: free convection 
boundary layer flow and heat transfer over a two dimensional 
heated semi- infinite horizontal plate* 



163 


REFERENCES 


[ 1 ] 

[ 2 ] 

[3] 

[^3 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[ 8 ], 

[ 10 ] 

[ 11 ] 

[ 3^3 

[13] 


Stewartson^ K., Z.A,M*P., 9a 276-282 (1958)4 

Gill, W.N*, Zeh, D.W. and del-Casal, E., ibid 16, 539-541 
(1965). 

Rotem, Z., First Canadian National Congress of Applied 
Mech-, Proceedings, 2b, 309-310 (1967). 

Rotem, Z,. and Claassen, L. , J, Fluid Mech, 39, 

173-192 (1969). 

Rotem, Z. and Claassen, L. , The Canadian Journal of 
Chemical Engineering, 47, 461-468 (1969). 

Pera, L, and Gebhart, B,, Int# J. Heat Mass Transfer, 

16, 1131-1146 (1973). 

Singh, P#, Ph.D, Thesis, I.I.T,, Kharagpur, India (1967). 

Verma. R.L. and Singh, P., Aust. J. Phys., 30, 335-345, 
(1977). 

Verma, A.R., Int. Jl. Engng. Sci. 2l, 35-43 (1983). 

Hansen, A.G., Trans. ASME 80, 1553-1559 (1958). 

Yang, K.T., Trans. ASME Jl. of Applied Mechanics, 27, 
230-236 (1960). 

Nanda, R.S., Proc. Symposium on Magnetohydrodynamics, 
I.I.Sc., Bangalore, India (1965). 

Jaluria, Y., Natural Convection Heat and Mass Transfer, 
Pergamon Press, Oxford, 79-84 (1980). 



CHAPTER V 


LAMINAR FREE CONVECTION OVER NON-ISOTHERMAL 
SEMI-INFINITE HORIZONTAL PLATE 
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I, INTRODUCTION 

In this chapter we study the laminar free convection 
£lov/ and heat transfer over a non- isothermal semi-infinite 
horizontal plate with an arbitrarily prescribed surface 
temperature distributions. In free convection over horizontal 
surfaces/ most of the previous studies deal with either 
uniform surface temperature distribution or the surface 
temperature varying according to power-law in streamwise 
distance from the leading edge* These cases are easily 
amenable mathematically because they yield to similarity 
solutions [ l,2/3/4/5/6/7] • Though these conditions are 
readily realizable in the laboratory^ however arbitrary 
non— isothermal surface conditions are expected to exist 
in most physically realistic situations* For the arbitrarily 
proscribed non-uniform surface temperature distribution for a 
semi -infinite vertical plate, the free convection flow 
and heat transfer has been studied by many authors [8,9,10, 

II, 12,13,14] • Unlike the vertical plate problem where the 
surface characteristics to be determined relate to shear- 
stress and heat-flux, the horizontal plate problem includes 
also of determining the induced pressure at the wall. For 
vertical problem, methods ranging from Karman-Pohlausen 
technique with relatively crude profiles to full scale 
numerical computation procedures have been developed. On 
the other hand, to our knowledge, no single study either 
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analytical or experimental has been attempted so far for the 
corresponding case of a horizontal plate. Therefore, wo 
initiate here this study and give some exact solutions of 
the governing boundary layer equations when the surface 
thermal condition is non*“unifonn, with the exception of the 
cases that allow similarity solutions. 

In what follows, a Gdrtler-type transformation is 
introduced to treat the problen of non- isothermal free 
convection with general wall temperature distribution. A 
new parameter characterizing the non- isothermal conditions 
of the wall is introduced. The resulting transformed non- 
similar form of the boundary layer equations are solved by 
the Gdrtlor series expansion [IS] and the local non- 
similarity method [16 ] ». The series expansion is valid 
for an arbitrary surface temperature variations given by 



whore n,p,Gj^ arc constants, / 0 x is a suitably 

normalized streamwise coordinate. Such wall temperature 
variations are generally encountered in experiments. In 
the expansion of flow variables as a series, the zeroth- 
order term in the expansion satisfies all the boundary 
conditions in the complete problem and hence already 
represents a solution to the full problem. Succeeding 
terms in the series only serve to Improve the accuracy 
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of the solution inside the boundary layers# Also/ since 
equations governing higher order terms are linear^ they are 
integrated independent of the coefficients of wall 

temperature data (5#1) for a particular problem by splitting 
the higher order terms into universal functions# Universal 
functions associated with the first four terms of the series 
expansion for flow variables are determined for air when 
p != 1#0, 0#5ji 0*25, With the help of these functions^ wo can 
investigate the boundary layer flows by simple algebraic 
process# An application of the series solution method is 
illustrated for a class of wall temperature distributions 
given by 

°w ^o ® 

for some constants G and e* On the other hand/ the 

o 

local non- similarity method at the second level of truncation, 
namely the two-equation model, is applied for this special 
case to predict local heat— transfer coefficient at the wall 
for Pr » 0.1, 0.72, 5.0 when p = 1/ 0.5, 0.25. A 
comparison of results for local heat-transfer at the surface 
for air indicates that the local non-similarity method gives 
accurate results for a larger region in the downstream flow 
covered here than the four-term s\jm of the series e:q>ansion 
can give. The effect of non-isothermal wall is significant 
on hoat-transfor at the surface and the variations in heat- 
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transfer at the surface is greatly influenced by Prandtl 
number and the shape of variation of wall temperature 
distribution as given by values of p* 

2 . BASIC EQUATIONS 

We consider the steady two~ dimensional boundary layer 
flow on a semi- infinite horizontal plate placed in an infinite 
expanse of Boussinesq fluid. The leading edge of the plate 
is taken as the origin of coordinates / x is along the plate 
and y normal to the plate* The governing continuity, 
momentum and energy equations are 

+ §1 = 0 » (5.3) 

3x ay 

u^ + v^ = -- ^ + V , (5.4) 

3x ay ^ 3x ay^ 

i ^ (T - Tj , (5.5) 

ay 

+ = a (5.6) 

ax ay ay 

where u,v are the velocity components along x and y 
directions respectively, p is "th® dynamic pressure, ± 
is the temperature, T^ is the temperature of the ambient 
fluid, 3/ a. a and g are respectively the kinematic 
viscosity, density, thermal diffusivity, the coefficient of 
thermal e 5 q)anslon and gravitational acceleration. 
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The boundary conditions of the problem are ■ 


Y = Oi u=v = 0. 


T,(x) ^ 
w 


(5.7) 


U ^ Og p 0, T T 


whore T^( x) is the temperature of the plate. 

By introducing the following non-dimensional variables 


v = - V - 2 11 - iili VL 


P 


2 

pL 

9 * 

pv^ 


T 


V 


T„(x)-T„ ' 


y 


(5,8) 


wo transform the equations (5,3 to 5,6) in dimensionless form 


&LE + 0 / 

ax ay 


' (5.9) 


uSii + v-SH + ilii , 

ax ay - ax 


5LE s G T , 

ay w 


uiS + v^ + ~ — uT=i ^ 
ax + ^ ay + dx a 


ay 


where G is a variable Grashof number defined by 
w 


G 


(\(x) - Tj 


w 


V 


2 


(5.10) 

( 5 .1 1) 

(5.12) 


L denotes a characteristic length and a is the Prandtl 
number of the fluid* 
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The boundary conditions reduce to 
y=0: u = v = 0# T = l} 

y-*.oo; U-*0, p-^O/T-^O 


( 5 •! 3) 


3 , transformation of boundary layer Equations 

The non- isothermal condition in terms of the function 

G still remains in the boundary layer ec^uations (5*9 - 5 •12)- 

w 

In order to achieve a description in which the function 
no longer appears explicitly, we introduce in the equations 
the following transf03anations 


t = S G,, dx # 


and 


G 


3/5 


w 








G' 




w 




T = 


(5.14) 


* is the stream function and satisfies the continuity 

equation ( 5 *9 ) • 

The houndary layer eolations (5.10 - 5.12) nov. 


become 
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■.m 


(n+ 1 ) 


[(3-20) ff " (1+0) f'^+ (2-30) g*-2(l+0)g] 


= (n+l)g [f' , ( 5,155 

^'= ® ' (5.16) 

i 0 " + i~il(3-20)fe'-(n+l)0ef'= (n+l)g[f' ||- -e' ||] , 

dg 3 § 

(5.17) 


In the foregoing equations, the primes denote partial 
differentiation with respect to V and 0(g) is the 
GORTLER'S Principal Function of the problem which 
characterizes the non-isothermal surface condition and is 
given by 


3(0 


i- 

g2 

w 


dG 
w 

cix 


(5.18) 


'Ihc boundary conditions for the transformed equations are: 


f(O0) = f'(g, 0 ) = 0 , e(g,o) = 1 . 
f'(§/«>) a 9 ( 0 °®) =* 0(0“®) - 0 • 


(5.19) 


In the case of power-law distribution of wall-temperature, 
the equations (5*15 — 5.17) reduce essentially to those of 
the similarity equations dealt with by Rotera and Claassen [ 4 ] . 
To solve the pairtial differential equations (5.15 — 5.17)# 
we employ first G&’rtler series expansion technique and then 
local non-similarity method. 
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4, GORTLER SERIES EXPANSION METHOD 

To solve the flow problem by Goirtler series method we 
first invert the transformation 


§ 


X 

S G , dx . 
o 


(5,20) 


with G given by (5*1)< by using Lagrange's theorem for 
reversion of power series* This inversion process gives x 
as a power series in § which we substitute in (5*18) to 
get the following series for 3: 

+ ••• < 5 . 21 ) 


whore are to be fixed by Imowlng the coefficients 

G^, Gj, G^-.-ln (S.1) for a given surface temperature 
distribution. Ibe procedure to obtain 3(5) as a power 
series in 5 follows the same way as shown in the appendix 


of Chapter IIA* 

To solve the eguations (5.15 - 5^7) we assume the 
following series for the functions f(6.’l). gCS.’l) 

ecg,^); 


oo 

fCS.h) = E V”' 

6*’’' ' 

Tc=0 


OO 

g(5.n) = 9^'”’ 


eC6,h) - ^ 

1c«0 

5*’’^ - 


(5.22) 



Substituting equations (5.21) and (5.22) in equations 
(5*15 to 5.17) and collecting the coefficients of like 
powers of , we obtain the following set of ordinary 
differential equations in f^(T))* 


f^" + . iR+ll f'^+ ilral rj gj- |(2n+l)g^= 0/ 


1. p " 4. (bt? ) 

a ®o ^ 


= ®o ' 


= — fe'-nf^e =0^ 
5 o o o o 


(5.23) 


and 


(n^ fo^ - [ 2 + pWn+D] +pkCn+l>]f;\ 


./ . r (n+3) 

4 . iazSlr) a' - a 

+ —5 V g]^ - 5 % 

= (2f^f- + f'^+ ZT) gj + 2g^) 


k-1 


[ {pj(n+l)+ f5£]^.j-CpJ(n+l)+ 






Sk = ®lc' 


1 + ("+31 f e'- [n+(n+i)pic] f;e^+[^+ ("+i>pdeSV"8o*i; 

O' ^ ^ 


Pv(2^o8i =*o®o’ 


k-1 
+ S 
j=l 


[£(n+l)p+n!f5e^.j-{^+ ("+1’P> 

for k = l#2#3, ••• 


(5 .24;) 
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The boundary conditions are*. 

f^(o) = f'(o) = 0 , e^co) = 1, 

f;(~) - 0, g^(») - 0, e^(») - 0 (5.25) 

and 

f^^co) = f'.(o) = e^(0) = 0, 

- 0, g^(cc) ^ 0, eyX<^) - O (5.26) 

for ]c = 1,2/3/ ... 

It may be noted that the zero-order ecjuations (5*23) are 
the similarity e(guations for varying as where the 

constant exponent n is to be determined by the actual 
behaviour of as x approaches 0. In physically 

realistic situations the plate temperature sremains finite 
and hence the leading edge behaviour of G^ can be very 
adequately described by x^ with n > 0^ The rest of the 
equations aire linear ^d form a recursive set of equations# 
These linear equations depend upon the prescribed surface 
thermal conditions through the appearance of pj. To 
eliminate this dependence, we write functions f^^, 
k > 1 as linear combination of certain other functions: 

% = ^ * ®1 “ 

^ 2 -= '^ 2 ^ 22 ^’^^ ^ ^ 2 * ^ 2 ^ 22 ^^^ ^ 1 ^ 21 ^^^ ^ 

©2 * ^ 2 ^ 22 ^^^ ^ 


4 ? 
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^3 “ ^ 1^31 ^ 2 ^ 32 ^ 3 ^ 33 ^^^ 

+ + &i^2^^2^Ti)*, e^ s ^3^32^'n'i -* 

+ %'^2^32^^^ (5*27) 


and so on* The functions y,Z and W are universal functions 
of V with parameters n,p and a* These universal functions 
satisfy tho following sot of equations: 


(n+1 ) 


T 


(2f C + 3-V g' + 2g^+ f' ) # 


o o 


Hi =^ii' 


o o 


(5*28) 


1^ *22' 22^ 22-' 


Z' 

22 


(2fQf''+ f'^t 3jr? g' + 2gQ ) / 


W. 


22 ' 


(5*29) 


1-2 (^22# 232 '^22^ « ^2fQe' + 5 f ^ e ^) t 

.p(n+l)]Y,,Yri 


+ |(n+l)(Yiif'+£oYl"l)+ |(n+l)Zii+ |(n+l)zii 


2^1 


^2 ^^21 '^21' ^21^ 


= w. 


21 ' 


= [p(n+l)- 

+ |(n+l)yiie'+(n+l)y{ie„+ f(n+l)f^e£i 



175 


L, - i-iL+l) 


l'^33'S3'”33' = ~T~ '2^"+ f' + 3D g' + 2g^) , 


^33 = «33 ' 


S,<''33'^33'”33> = ^ '2fo®o + ’ 


(5.31) 


= -[p(n+l)+ S~] Y^^Y^'+[3p(rn-l)+ |( 2n+l )] Y'^Y'j_~ 
-[2p(n+l)+ i£^]Y^'Y 2 ,+ li3^ 


(n+1 ) 




7' 


''31' 


Lj,(Y 33_,Z3^,W33,)=-[-^^2+31 + p(n+i)]Y^Z''+ [n+ 2 p(n+l)] - 


7r2l 


«[i2±ll + 2p(n+l)] Y2^Z^'+[n+p(n+l)]Y'J_z'J_-i- 

+ |(n+l)Y^^'^*+(n+l)Y£zJ^+ |(n+l)f^z'{ + 


+ (n+l)fjz^l+ |(n+l)g''Y2^+(n+l)gSY2l > 


(5-32) 


^^32*^32*^32^ 


.[p(n+l)+ +3p(iM-l)]^'Y- 

•[2p(n+i)+ ■'■ 

. |(n+l)Z22+ f(n+l><fi'^22+^J''22«o’'5P ' 


Z32 ="32' 
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"^3?/ ^^32 '^^32^ = + p(n+l)]Y^^^'+[n+2p{n+l)]Yj^'2 - 

4in±3l ^ 2p(n+l)]Y22Z{{+[n+p(n+l)]Y2'2 2'j_ 
+ ~(n+l) 

+ ( n+1 ) PqY' 2 + f ^ ^ ^0^11 ^ ^0^11 '*' 

+ |(n+l) Y^^^'+ (n+l)gjY'i (5.3; 


v;here 


'■l^^ki-'^kr'^Kr' 




+ [li!^ + kp(n+l)]f‘'V 

^o\r ■*■ 


[in+^ + (n+l)pk] n 0 ^ Y'^ ( 5 . 34 ) 


for k > 1 and r « l,2/3t...#k* 
The boundary conditions ares 


and 


Y^^(0)= Xyo) = V'°’ “ ° 

- 0, W^(“) - 0 » V‘”* * ° 


for k > 1 


and r « l|2/3^ V^k. 


( 5 - 35 ) 
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The above system of equations for determining universal 
functions can now be easily solved for the various set of 
values of parameters n/p and ct« Hence the calculation of 
flow and heat transfer characteristics of the free convection 
boundary layer over a no n-' isothermal horizontal plate becomes 
a simple matter* 

A standard shooting technique as described by Nachtshelm 
and Swigert has been employed to solve the zero-order 

equations for n = 0 and a = 0.72. The solutions 
f'^(O), g (O), 6_(0) are calculated to an accuracy upto 

w V/ 

six decimal places. The linear equations (5.28 to 5 ..33) 
are solved recursively upto third-order terms in the series 
expansion (5.22) for p = 1.0,0.5#0.25 with n = O and 
a =: 0.72. The method of solution adopted here is the same 
as given in Chapterll except for the fact that here we have 
to consider cjvaluations of complementary functions thrice 
with throe sets of initial conditions at the surface instead 
of two for the problem in Chapterll. The equations v/ore 
integrated by Runge-Kutta integration scheme with step-size 
H « 1/64 and the edge of the boundary layer fixed at = 15. 

This yields 


Throughout the computation/ 


calculations were carried out 
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upto eight decimal places « For the Prandtl number# a = 0.72 

and n = 0 the universal functions are determined upto 

and including the third order terms for values of p = 1.0, 0,5,0 -.2'S 

The missing wall values for the functions f . g ; 0 and 

o o o 

the universal functions Y, Z,VJ are shown in Table 1, 

The rate of heat transfer from the plate to the fluid 
is the most important characteristic of the free convection 
boundary layer problem^ which is calculated by using 
Fourier law: 


q 


3y y=0 


0 


Introducing the dimensionless variables^ the expression 
for q becomes 


q * 


^ (T (x) 
L w 


,3/5 


T ) 


w 

[Cn+Dg] 


30(5,0) 

an 


I cy x) - T„) 




[ S 

k=0 


an ^ 

(5.36) 


NOW introducing the local heat transfer coefficient, 
local Nusselt-number, and local Grashof number as 

.r. 

h = -S , Nu_ = p , Sr- = ~ ' 
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th© dirriGnslonloss rspjrsssntcition of the local heat tiransfer 
becomes 


Nu 


X 


Gr 
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X 




(5.37) 


5. LOCAL NON~SIMIIARITY MEIHOD 

We apply this method of solution for the wall-tcamperature 
distribution given by (5.2). Before proceeding with the 
formulation of the system of equations for local non-similarity 
mothod^ WG shall obtain firsts, 3(S) iii terms of p and 
s, whore 

s-l+e^. ^5.38) 


using cejuations (5.2) and (5.38) in the definitions of g and 
m) given respectively in (5.14) and (5.18), we obtain 


§ 


(s-1) 


1/p (p-fs) 

TiSl 



(5*39) 


and 


^($) 


p(p+s)(S"'ll 

(p+1) 


Also it can be shown that 


(5.40) 


g ^ = s^(s) 


(5.41) 
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In order to describe our results for heat transfer coefficient 

G 

as a function of ^ (=s)/ we transform the non-similar form 

"'o 

of the boundar/ layer equations (5.15 to 5 .1?) in terms of s 
and V instead of [; and n using (5 ••41). Then equations 
(5*15 to 5*17) become 


f'" 4- i [(3-2l3}ff^' + (2-38)7) g'- 2(l+'3)g] 

= s8(s) [f'l'- f"f + g ] , 
g' = © # 

1. J. (3-28)f0'-8©f''= s3(s) [e f'- f 0'] 

0 5 


(5.42) 

(5.43) 

(5.44) 


where 




# g 


, 0 =* 
3s' 


3 ^ 

3S 


(5.45) 


The boundary conditions (5.19) now take the form 

f(s,o) = f'(s,o) = 0, e(s,o) = 1 - 
£'(s.“>) = g(s/“) = e(s,“) =0- 


(5.46) 


Now we proceed with the formulation of the wo 
Eqfuatlon model in the local non-simxlarity me 
solution. Tho derivation of this model follows to the 
same manner as descrlhed In ChapterllX for free convection 
houndar. layer flow over a non-lsothermal plate Inclined at 
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(1 stTiQll an^le to the horizontal and hence the details of 
the method are omitted here* 

The system of equations for the two-equation model 
can then be summarized as: 

(a) Equations (5,42 - 5.44) 

(b) the truncated equations for f|g and 0 


fl". [s II + f'+ [s If + I 


I- |(2-3P)»I s’- [ s a| + I (2+7B)]g = 


g s e 


i.||(2ff"*+f' +37)9' t2g) , 

(5.47) 

(5.48) 


1 g'% I (3-23) f [s H + |(1+P)3 e' f - [s II + 23]f^C-3f 0 


^(2f0'+ £'0) • 

5 ds 


(5.49) 


(c) Tlio boundary conditions 

f(s.o) =. =?(s,0) =l'(s. 0 ) =e(s,o) =0 , 

e(s,o) « 1’/ 

g(s,») = ° 


(5.50) 



182 


It c.in bo scon that the equations (5.42 ~ 5*44) and 
(5.47 “ 5*49) arc coupled and they rmist be solved for the 
six unknown functions f ,g,8#f,g,and 8. With s regarded 
as known^ 8 and ^ can then be determined and therefore 
those equations may bo treated as a system of coupled 
ordinary dif ferential equations of the similarity typo. 

The system of equations for the second level local 
non-similax'ity model was solved by emi^loying the initial 
value technique as modified by Nachtsheian andSwigert. 

A solution was considered to be convergent when the magnitudes 
of the boundary conditions at the edge of the boundary layer 
g( 3 ,oo)i e(s,co) and £'(s,«o), g(s,«), e(s,«>) for 
a fixed s became less than 5 x 10 ^ and 5 x 10 
respectively. In table 2^ we give the values or 
which dosoribes tho dtaiGnslonlESS characteristic of surface 
hoat-flux for a = 0.72, 0.1, 5.0 and p = 1, 0.5, 0.25. 


6 • 


results and discussion 


Wo shall consider the application of Gortlor series^ 
method of solution to the class of wall temperature variati 
= G^(l+ s as already given in (5.2). The exponent p, 

which gives tho shape of the variation, is required to be 

i 1 nr^n—intGcrx‘9.1 • Also wc noto 
a positive number/ integral or 

, ■e+.q'r* 1 enath is taken as the length 

that when the characteristic leng 

p renresents the maximum deviation of 
of the plate l, then e represen 
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from in the region x = 0 to x = ; • Here we 

have considered the wall temperature relative to the ambient 
to be .1 finite non-zero value at x = 0 and this increases 
steadily trom x = 0 to x = i!,. 

Now comparing the class of surface temperature variation 
( S * ?. ) with (5*1) we get 

n = 0 / Oji^ = e ^2 ~ ^3 “ ••• = 0 * 

Using those values and the procedure as followed in Chapter II 
v/c obtain the coefficients 0 ^ of the Principal function 
) as followsi 


n 


n ^ ** 

O nHhl 


0 I 


^ gP 

o 

02 « - 1^2p * ^ 


« [ (245p) ^ + (l+4p)] Xp^ q3p ' 

P 

- -[(3+7p)(^)^ + |(4+wp) pS: 


i'r«n tho definition of S as given in agnation (5^4), 

after simplification^ we get 


O 

(5.51) 


(5.52) 
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Using (5.52) and (5.51)^ the dimensionless representation 
of the local heat flux given by ecjuation (5.37) can be 
written in terms of s for a given p as follov/s; 


Nu^ 

X 

X 


p-hs j 


a-gjg/O) , 

371 


(5.f3) 


= - [ [0o(O)+^3^W'^(o)§P+{32w/^(O)4.p2W'2(O)}|2P 

+ {^^W'^{0)+33_32W32(0)+%W'3(0)}g2P-h...] . 

(5.54) 

The Nusselt-Grashof relation as obtained frcm (5*5 3) 
and (5.54) for the local non-sinilarity method and Gdrtler 
series expansion respectively is shown in figure 1 over the 
range 1 £ s < 1.4 for a = 0.72P and p =1 .0/0.5 #0.25. 

A comparison of local heat transfer at the surface shows 
that the agreement between the two methods of solution for 
a fixed s depend both on the number of teims used in the 
series solution and values of p. For a fixed p we find 
the region of agreement between results for heat— transfer 
from the two methods of solution Improve if 4- term sum 
instead of S-tejsn sum is used in Gdirtler series solution. 

In the region R ^diere the results from the two methods 
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deviate there is a significant trend observed for all values 
of p in using a 4- term sum and 3- term sum from the series 
solution for comparison purposes* The trend of results from 
3^term s\m for local heat transfer in R is to underestimate 
it while, the 4-term sum overestimates even as the region R 
shrink now a little. Also this alternating behaviour of 
Gortler series in relation to inclusion of additional terms 
for obtaining solution in downstream flow is to continue, 
judging by the trend of the results for universal functions 
and the value of coefficients. However, the evaluation 

of further universal functions is hindered by considerable 


difficulties. These are due not only to the fact that for 
every additional term in the series the number of differential 
oguations to be solved increases but also and even more 
forcibly the difficulties are due to the need to evaluate 
the functions for lower power term with ever increasing 
accuracy, if the functions for higher power terms are to be 
sufficiently accurate. Also the chances of carrying out an 
extens-ion of this method is limited by the fact that a 
corresponding demand on the need to do increasing algebraic 
manipulations for finding Sj anerges. On, the other hand, 
in the tw -Equation model of local non-slmilarlty method 
solving of six coupled equations, simultaneously by initial 
value techniques may show divergent behaviour unless good 
estimates of the missing wall values are Icnown, Fortunately. 
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In the cases treated here such couplings have mostly 
insignificant influence on the convergence of the solution.. 

We mention here that quasilinearization technique can be 
of value to predict good estimates of the missing conditions 
at the wall as it has already been demonstrated in Chapter III., 

The results from local non-similarity method for the 
Nusselt-Grashof relation as obtained from (5.53) are 
presented in figure 2(a), (b) as a function of s for 
O’ — 0.1, 5.0 and for p = 1.0, 0*5, 0*25 • For a fixed o 
wo see from figures (1) and (2a, b) that the effect of increase 
in value of p at a fixed s isto Increase the local Nusselt 
ntamber.. 

The effect of a non-isothermal wall on heat transfer at 
the surface is computed from the foUowipg equation 


" f^) Gr^^] (Nu yorl/^) 

X _ X X 

[q/(\- f^) (NuyGr^^y 

” X iso X X iso 

where the subscript, iso, denotes a quantity corresponding 
to an isothermal wall. For our special case (5.2) we can 


q 


(p+s)s 2/5 

[.£_] 

(p+i) 


r dV 

e;(o) 


] 


write 
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where 0^(0) is the characteristic dimensionless surface 
heat flux for isothermal vTall* Itie values of against s 
for <y = 0*72'/ 5»0 as obtained by using the two-equation 

model of tluo local non-similarity method are shown in 
figures (3 a),b)^c)) for p = 1.0, 0,5, 0,25 respectively.. 
For a fixed p, the percentage increase in heat transfer 
coefficient at the surface over an isothermal wall for a 
fixed s decreases as the Prandtl number increases,. . 
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Table 1 

Missing wall values for universal functions for 
a = 0,72, n = 0 and p = 1,0.5,0*25 



0.25 


f'' (O) 

a'(o) 

Yi'(o) 

W'^{0) 

Y''(0) 

'’’22 

W'^CO) 

Y''(0) 

7^1(0) 

W'lCO) 

Y'|(0) 

Wj3(0) 

Y”(0) 

Zj^CO) 

WS 2 ^°> 

;w|i(o) 


0.978404 
-1 .734925 
-0,357409 
0,084434 
0.093690 
-0.181933 
-0.074718 
0.121104 
-0,154230 
0.013412 
-0.001766 
I 0.025947 
0^34195 
0.045798 
-0.141510 
0.009101. 
-0..011420 
0.X337064 
-0.002411 
0.002999 
-0.004009 


0.978404 
-1 ,734925 
-0,357409 
0.052887 
0.127810 
-0,201790 
-0.105 5 69 
0.174429 
-0.184974 
0.014560 
-0.014304 
0,042405 
0,103918 
0.073588 
-0.167873 
0.010961 
-0,028079 
0.064116 
-0.001616 
0*005201 
-0.007265 


0.978404 

-1.734925 

-0.357409 

0.027968 

0.154982 

-0.215486 

-0.132497 

0,222'341 

-0.209491 

0.0l7082r 

-0.029020 

0.057273 

0.070809 

04O829I 

-0.190900 

0.014678 

-0.049559 

0,091688 

-0.001687 

0.009876 

-0.011573. 
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Fig. 3 Effect of non-rsothermal wall 
condition on heat transfer 
(local non-similarity method). 
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